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A comprehensive  survey  of  improper  ferroelectrics  was  carried 
out  to  identify  materials  with  potential  to  enhance  the  perform- 
ance of  the  pyroelectric  vidicon.  Representatives  of  four 
families  of  improper  ferroelectrics  were  selected  for  detailed 
experimental  and  theoretical  analysis:  rare-earth  molybdates, 
the  sulfate  family,  propionates,  and  boracites.  A thorough 
experimental  and  theoretical  characterization  of  the  rare-earth 
molybdates  confirmed  the  expectation  that  the  dielectric  anomaly 
could  be  nearly  absent  in  improper  ferroelectrics  and  thus  lead 
to  a high  value  of  the  f igure-of-merit  p/{cc.  Experimental 

confirmation  of  the  theoretical  models  allows  one  to  obtain  an 
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optimized  f igure-of-merit  of  10  x 10  cm  /C  for  materials  of 
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this  family  as  compared  to  3 x 10  cm  /C  for  TGS  at  40°C). 

From  the  sulfate  family,  lithium  ammonium  sulfate  and  its 

deuterated  analog  exhibited  a room-temperature  value  of  the 

f igure-of-merit  about  equal  to  that  of  TGS  at  40°C.  The  other 

analog  sulfates  examined  exhibited  a characteristically  different 

behavior  and  smaller  values  for  the  f igure-of-merit Initial 

measurements  on  the  propionates  were  analyzed  and  lead  to  a 

projected  figure  of  merit  which  is  unbounded  above  and  was 
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estimated  to  be  7.3  x 10  cm  /C  for  dicalcium  lead  propionate. 
Initial  measurements  on  boracites  were  consistent  with  data 
measured  elsewhere  which  give  the  expectation  of  a figure  of 
merit  for  the  boracites  greatly  in  excess  of  that  of  TGS.  Thus, 
the  expectation  that  improper  ferroelectrics  can  have  large 
figures  of  merit  is  confirmed  in  detail;  the  boracite  and 
propionate  families  were  identified  as  having  the  greatest 
promise  for  enhanced  device  performance. 
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PREFACE 


A program  to  develop  an  optimized  pyroelectric  vidicon  thermal 
imager  was  issued  by  the  U.  S.  Army  Mobility  Equipment  Research 
and  Development  Center,  Fort  Belvoir,  Virginia,  and  was  initiat- 
ed under  Contract  DAAG53-76-C-0053 . Modification  No.  P00002  to 
the  contract  added  the  requirement  for  research  of  improper 
ferroelectric  materials;  Dr.  Lynn  Garn  was  the  Contracting 
Officer's  Representative  for  the  Night  Vision  & Electro-Optics 
Laboratories,  USAECOM. 

This  work  was  performed  by  Philips  Laboratories,  a Division 
of  North  American  Philips  Corporation,  Briarcliff  Manor,  New 
York.  The  materials  research  task  was  under  the  overall  super- 
vision of  Dr.  Stewart  Kurtz,  Director,  Exploratory  Research 
Group.  Dr.  Wallace  Arden  Smith,  Theoretical  Physicist,  was  the 
Program  Leader.  Dr.  Gertrude  Neumark,  Physicist,  participated 
in  the  theoretical  analysis.  Measurements  at  Philips  Laboratories 
were  made  by  Dr.  Joseph  P.  Dougherty,  Material  Scientist; 

Mr.  M.  Delfino,  Chemist;  and  Mr.  G.  M.  Loiacono,  Senior  Program 
Leader,  Crystal  Growth  & Materials  Technology.  Mr.  Loiacono  was 
also  responsible  for  materials  growth.  Dr.  Edward  Stupp,  Senior 
Program  Leader  for  Solid  State  Components,  was  available  in  a 
consulting  and  advisory  capacity  on  device  implications. 

A subcontract  was  let  to  Yeshiva  University  to  perform  pyro- 
electric and  dielectric  measurements  and  to  collaborate  on  the 
analysis.  Prof.  Michael  I.  Bell,  Principal  Investigator, 

Dr.  Avner  Shaulov,  Visiting  Scientist,  and  Mr.  Y.  H.  Tsuo, 

Research  Assistant,  conducted  this  research. 

The  work  described  in  this  Final  Report  covers  the  period  from 
September  1,  1976  to  February  28,  1978. 
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1. 


INTRODUCTION,  SUMMARY  AND  CONCLUSIONS 


1 . 1 Introduction 

This  study  of  improper  f erroelectr ics  was  motivated  by  a desire 
to  identify  materials  which  would  enhance  the  performance  of  the 
pyroelectric  vidicon  (PEV).  The  pyroelectric  materials  commonly 
used  in  this  device  are  in  the  TGS-DTGFB  family  and  are  all 
proper  f erroelectrics . The  principle  limitation  with  proper 
f erroelectr ics  is  that  the  increase  in  the  pyroelectric  coef- 
ficient obtained  by  operating  closer  in  temperature  to  the 
f erroelectr ic-paraelectric  phase  transition  is  more  than  over- 
come by  a concomitant  increase  in  the  dielectric  constant. 

Thus,  the  material's  figure  of  merit,  p/ec,  where  p = pyro- 
electric coefficient,  e = dielectric  constant,  and  c = volume 
specific  heat,  decreases  as  the  ferroelectric  Curie  temperature 
is  approached.  In  the  case  of  improper  ferroelectr ics , there  is 
no  dielectric  divergence  upon  approaching  the  ferroelectric-para- 
electric  phase  transition,  and  therefore  one  obtains  a figure  of 
merit  that  increases  dramatically  up  to  the  Curie  point. 

The  investigations  at  Philips  Laboratories,  in  conjunction  with 
Yeshiva  University,  began  with  a comprehensive  literature  survey 
and  selection  of  representative  materials  from  a number  of  known 
classes  of  improper  ferroelectr ics . Samples  of  each  material 
were  obtained  or  grown,  and  a comprehensive  set  of  thermal  and 
electrical  measurements  was  performed.  The  measurements 
program  consisted  of  independent  determinations  of  p,  e,  c,  p/c 
and  p/ec.  With  five  independent  measurements  of  three  quan- 
tities, we  were  thus  able  to  eliminate  the  discrepancies  and 
artifacts  that  obscured  the  published  literature.  Combining 
these  precision  measurements  with  an  extended  thermodynamic 
theory,  we  were  able  to  confirm  our  initial  expectation  that  the 
improper  f erroelectrics  would  have  high  figures  of  merit.  In 
every  class  of  improper  ferroelectr ics  investigated,  we  found 
the  figure  of  merit  was  comparable  or  distinctly  superior  to 
that  of  TGS. 
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1.2 


Summary  of  Results 


( 1 ) Molybdates 

A detailed  analysis  of  the  dielectric,  pyroelectric  and  heat 
capacity  measurements  on  gadolinium,  terbium,  and  gadolinium- 
dysprosium  molybdate  revealed  an  excellent  agreement  with  the 
theoretically  expected  behavior.  Moreover,  these  separate 
analyses  of  the  dielectric,  pyroelectric  and  thermal  properties 
can  be  combined  to  predict  the  behavior  of  the  figure  of  merit. 
Independent  measurements  of  the  figure  of  merit  are  in  excellent 
agreement  with  the  theoretical  projections. 

These  measurements  and  their  analysis  confirm  the  theoretical 
expectation  that  improper  f erroelectrics  can  overcome  the 
limitations  in  the  usual  proper  ferroelectrics  brought  about  by 
the  divergence  of  the  dielectric  constant  near  the  transition 
temperature.  The  improper  ferroelectrics  have  well-behaved 
dielectric  constants  and  thus  have  a figure  of  merit  that  rises 
as  the  transition  temperature  is  approached. 

Using  the  theoretical  expression  for  the  figure  of  merit,  we  are 

able  to  optimize  this  expression  with  respect  to  variations  in 

the  di f f icult-to-assess  free  energy  parameters.  The  optimum 

figure  of  merit  is  then  given  in  terms  of  readily  measured  or 

estimated  quantities.  Using  values  characteristic  of  the 

molybdate  family,  we  estimate  an  optimum  figure  of  merit  of 
3 2 

about  10  x 10  cm  /C  (for  comparison,  we  find  the 
figure  of  merit  of  TGS  to  be  3 x 10^  cm^/C  at  40°C) . 


( 2 ) Sulfate  Family 

Five  materials  from  the  family  of  ferroelectric  sulfates  were 

studied:  LiNH.SO.,  LiND . SO . , LiCsSO, , LiKSO . and  LiHo(Se0o),. 

44  44  4 4 3 32 

Lithium  ammonium  sulfate  and  its  deuterated  analog  exhibited 
similar  properties.  The  high  room- temperature  value  for  the 
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figure  of  merit  (about  equal  to  that  of  TGS  at  40°C)  suggested 
further  study  of  related  materials.  Both  lithium  cesium  sulfate 
and  lithium  potassium  sulfate  were  consequently  investigated; 
these  materials,  however,  exhibited  distinctly  different  behaviors 
which  eliminate  them  for  practical  pyroelectric  applications. 

Lithium  trihydrogen  selenite  was  investigated  because  of  its 

2 

high  value  of  spontaneous  polarization  (~  15  yC/cm  ),  and  the 
fact  that  the  related  material  RbH3(SeC>3)2  has  been  identified 
as  an  improper  ferroelectric.  We  find  no  evidence  that 
LiH3(SeC>3)2  is  an  improper  ferroelectric.  The  measured  value  of 
that  ratio  p/e  is  about  one-fourth  that  of  TGS.  While  this 
eliminates  LiH3(SeC>3)2  for  practical  applications,  the  large 
value  of  the  spontaneous  polarization  suggests  that  related 
materials  might  have  potentially  useful  figures  of  merit. 

( 3 ) Propionates 

Initial  measurements  of  the  pyroelectric  and  dielectric  proper- 
ties of  dicalcium  lead  propionate  (DLP)  reveal  the  presence  of 
an  internal  bias  field.  Annealing  samples  of  DLP  did  not 
eliminate  the  effects  of  internal  bias  fields  as  reported  in  the 
literature;  indeed  the  properties  of  the  samples  were  severely 
degraded.  Initial  measurements  on  the  specific  heat  of  DLP  and 
dicalcium  strontium  propionate  (DSP)  showed,  as  expected,  no 
large  anomaly  at  the  transition  temperature. 

A theoretical  analysis  based  on  a free  energy  appropriate  to  a 
second-order  pseudo-proper  phase  transition  was  carried  out;  the 
predictions  agree  with  our  measurements  as  well  as  certain 
published  data.  In  this  analysis,  the  figure  of  merit  has  a 
maximum  as  a function  of  temperature  which  can  be  expressed  in 
terms  of  the  free-energy  parameters.  Using  data  from  the 
literature  (which  is  compatible  with  our  measurements),  we  can 
estimate  these  free-energy  parameters  and  project  a figure  of 
merit  of  4.2  x 10^  cm^/C  for  DSP  and  7.3  x 10^  cm^/C  for  DLP. 


(4) 


Boracites 


Initial  dielectric  and  pyroelectric  measurements  were  performed 
on  samples  of  boracites  obtained  from  Plessey  Company.  Major 
difficulties  in  fully  poling  these  samples  were  encounter- 
ed. The  measurements  are  consistent  with  the  data  reported  by 
Plessey  Company  and  with  the  resulting  expectations  of  a high 
figure  of  merit  for  the  boracites.  Specific  heat  measurements 
were  made  on  samples  of  Zn-Br,  Cr-Cl,  Fe-I,  Ni-Br  and  Cu-Cl 
boracite.  In  addition,  infrared  vibrational  spectra,  electronic 
spectra,  and  second-harmonic  measurements  were  performed  on  some 
of  the  samples  in  order  to  better  characterize  the  material.  A 
number  of  initial  observations  can  be  made: 


• The  cubic-orthorhombic  phase  transition  appears  to  be 
strongly  first-order.  This  is  indicated  by  the  large 
latent  heat  of  the  transition,  a detectable  thermal 
hysteresis,  and  the  discontinuous  jump  in  the  second 
harmonic  intensity  at  the  transition. 

• The  different  members  of  the  boracite  family  may  be 
structurally  dissimilar . This  is  indicated  by  large 
differences  in  the  molar  heat  capacity  of  Fe-I  and 
Cr-Cl  boracite  in  the  cubic  phase.  Large  differences 
in  the  molar  latent  heats  also  support  this  con- 
clusion. 


1 . 3 Conclusions 

The  conclusions  of  this  study  are: 

• Confirmation  that  improper  f erroelectr ics  have  the 
expected  potential  for  greatly  improved  pyroelectric 
figure  of  merit. 

• Identification  of  the  boracite  and  propionate  mater- 
ials as  the  families  of  improper  f erroelectrics  with 
greatest  potential. 
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Interest  focuses  on  the  boracites  for  a number  of  reasons. 

First , measurements  (Ref.  1-1)  of  the  Fe-I  boracite  show  a figure 
of  merit,  p/ec,  which  is  one-fourth  that  of  TGS  at  room  temper- 
ature, and  which  subsequently  rises  to  exceed  the  figure  of  merit  of 
DTGFB  as  the  temperature  approaches  the  Curie  temperature  at 
72°C . This  behavior  is  characteristic  of  improper  f erroelectr ics . 
The  room- temperature  figure  of  merit  of  the  Ni-Br  boracite  was 
measured  (Ref.  1-1)  and  found  to  be  larger  than  that  of  TGS; 
one  would  expect  it  to  rise  dramatically  with  increasing  temper- 
ature to  result  in  a figure  of  merit  distinctly  superior  to  that 
of  DTGFB  (estimate:  a factor  of  3 greater).  Second,  the 
measured  (Ref.  1-2)  piezoelectric  constants  of  the  boracites  are 
one-tenth  those  of  TGS.  This  will  substantially  aid  in  tube 
ruggedization  by  decreasing  the  sensitivity  to  vibration. 

Third , pedestal-noise  suppression,  when  fully  implemented  in  the 

vidicon  tube,  will  change  the  f igure-of-merit  from  p/ec  to 
3/2 

p/e  c.  Since  the  principle  virtue  of  inproper  ferro- 
electrics  is  their  low  dielectric  constant,  this  will  further 
enhance  the  superiority  of  the  boracites  by  another  factor  of 
1.7.  In  addition,  the  measured  (Ref.  1-1)  resistance  of  the 
boracites  are  in  the  range  of  lO"^  - 1012  ohm-cm,  which 
could  possibly  eliminate  the  need  for  pedestal  generation. 
Elimination  of  this  noise  source  should  improve  the  overall 
tube  performance  significantly.  Finally,  the  boracites  are 
chemically  stable  and  mechanically  hard.  This  should  facilitate 
target  fabrication  and  simplify  tube  operation  in  a hard  vacuum. 

The  propionates  are  of  interest  for  three  major  reasons.  First, 

the  measurements  and  their  analysis  for  dicalcium  lead  propionate 

3 2 

( DLP ) allow  us  to  project  a figure  of  merit  about  7.3  x 10  cm  /C . 
Our  analysis  also  shows  that  DLP  belongs  to  a sub-class  of  im- 
proper ferroelectr ics  which  has  no  upper  bounds  on  the  potential 
figure  of  merit.  Thus,  we  expect  that  other  members  of  this 


Ref.  1-1.  R.W . Whatmore,  Annual  Research  Report,  Plessey  C.  Ltd., 
Nov.  1976-Oct . 1977,  Ref.  CD. 6500088. 

Ref.  1-2.  A.S.  Sonin  and  I.S.  Zheludev,  Sov.  Phys.  Cryst.  8, 

217  (1963). 
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family,  or  suitably  substituted  DLP , might  have  distinctly 
higher  figures  of  merit.  Finally,  the  growth  of  large  crystals 
of  materials  in  this  family,  using  standard  solution-growth 
techniques,  is  possible. 


2. 


MOLYBDATES 


A detailed  analysis  of  the  dielectric,  pyroelectric  and  heat 
capacity  measurements  on  gadolinium,  terbium,  and  gadolinium- 
dysprosium  molybdate  reveals  an  excellent  agreement  with  the 
theoretically  expected  behavior.  Moreover,  these  separate 
analyses  of  the  dielectric,  pyroelectric  and  thermal  properties 
can  be  combined  to  predict  the  behavior  of  the  figure  of  merit. 
Independent  measurements  of  the  figure  of  merit  are  in  excellent 
agreement  with  the  theoretical  projections. 

These  measurements  and  their  analysis  confirm  the  theoretical 
expectation  that  improper  f erroelectrics  can  overcome  the 
limitations  in  the  usual  proper  f erroelectrics  brought  about  by 
the  divergence  of  the  dielectric  constant  near  the  transition 
temperature.  The  improper  ferroelectrics  have  well-behaved 
dielectric  constants  and  thus  have  a figure  of  merit  that  rises 
as  the  transition  temperature  is  approached. 

Using  the  theoretical  expression  for  the  figure  of  merit,  we  are 

able  to  optimize  this  expression  with  respect  to  variations  in 

the  dif f icult-to-assess  free  energy  parameters.  The  optimum 

figure  of  merit  is  then  given  in  terms  of  readily  measured  or 

estimated  quantities.  Using  values  characteristic  of  the 

molybdate  family  we  estimate  an  optimum  figure  of  merit  about 
3 2 

10  x 10  cm  /C  (for  comparison  we  find  the  figure  of  merit 
of  TGS  to  be  3 x 103  cm2/C  at  40°C). 

2 . 1 Review 

This  section  reviews  the  previous  theoretical  and  experi- 
mental work  on  the  rare-earth  molybdates. 

2.1.1  Review  of  Theory 

The  theoretical  interpretation  of  gadolinium  molybdate  as  an 
improper  ferroelectric  was  advanced  independently  and  simul- 
taneously by  four  different  sources:  Dvorak  (Refs.  2-1  and 
2-2),  Levanyuk  and  Sannikov  (Ref.  2-3),  Aizu  (Ref.  2-4)  and 
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Pytte  (Ref.  2-5).  A subsequent  interesting  discussion  of 
the  theory  for  the  molybdates  is  presented  in  a paper  by  Dorner, 
Axe  and  Shirane  (Ref.  2-6)  as  part  of  a description  of  their 
neutron  scattering  experiment  on  terbium  molybdate. 

The  most  complete  analysis  is  contained  in  the  works  of  Dvorak 
(Refs.  2-1  and  2-2)  which  are  based  on  a knowledge  of  the 
precise  space  groups  of  the  high  and  low  symmetry  phases.  This 
evaluation  is  not  intended  to  detract  from  the  contributions  of 
the  other  workers  (Refs.  2-3,  2-4,  2-5)  who  inferred  the  es- 
sential features  in  the  behavior  of  the  molybdates  based  on  less 
complete  information. 

( 1 ) Dvorak  Analysis 

Dvorak  presents  his  analyses  of  the  behavior  of  gadolinium 
molybdate  in  two  papers  (Refs.  2-1,  2-2).  In  the  first  paper 
the  form  of  the  thermodynamic  potential  applicable  to  the 
molybdates  is  determined,  based  on  the  knowledge  of  the  precise 
space  groups  in  the  ferroelectric  and  paraelectric  phases  and 
the  observation  of  cell  doubling  at  the  transition.  In  the 
second  paper  (Ref.  2-2),  the  dielectric  and  elastic  behavior  of 
gadolinium  molybdate  is  deduced. 

Dvorak's  analysis  of  gadolinium  molybdate  is  a direct  application 
of  Landau's  theory  (Ref.  2-7)  of  continuous  phase  transitions. 
Three  experimental  facts  are  used  in  this  analysis:  the  space 
group  of  the  high  symmetry  phase  D^  (PT2^m);  the  space  group 
of  the  low  symmetry  phase  C®v  (Pba2);  the  volume  of  the  unit 
cell  in  the  low  symmetry  phase  is  double  that  of  the  high 
symmetry  phase.  In  a detailed  argument,  Dvorak  shows  that  only 
one  physically  irreducible  representation  of  the  high  symmetry 
space  group  can  induce  precisely  this  symmetry  reduction  at  the 
phase  transition.  This  is  a two-dimensional  representation 
belonging  to  the  M point  of  the  brillouin  zone.  The  form  of  the 
free-energy  (in  his  notation)  deduced  from  this  analysis  is. 
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F = Fq  + ^a(T)  (qj  + q*)  + (q}  + q*)  + |e2  qlq2  + B3qiq2  (ql  ' q2)  + 

12  231-171  2 2 1 2 

+ ?<qI  + q2V  + 5 *33  P22  ♦ 2 CU  (uxx  + V)  + 2 c33u„  + 

+ I c66uxy  * c12uxxV  + c13<"xx  + V “zz  + a36Pzuxy  + 

+ Y1  Pzqlq2  + r2Pz<ql  “ q2>  + 41  “zz(ql  + q2J  + 

+ 62|uxx  + V (qS!  + q2>  + 53uxyqlq2  + 54  Vql  ' q2> 


where  (<Ji»q2)  t*ne  soft  mode  amplitude,  Pz  is  the  z component 
of  the  polarization  and  are  the  components  of  the  strain 

tensor.  The  essential  difference  between  this  free  energy  and 
the  form  advanced  by  other  workers  is  the  presence  of  the  terms 
with  coefficient  g^,  Y2  and  6^.  It  is  precisely  these  terms 
that  cause  the  minimum  of  the  free  energy  in  the  ferroelectric 
phase  to  lie  off  the  lines:  q^  = + q2?  q^  = 0,  q2  £ 0;  and 
ql  ^ q2  = O’  That  the  physically  realized  states  in  the  low 
temperature  phase  occur  off  these  lines  is  essential  to  get  the 
correct  correspondence  with  the  observed  space  group  of  the 
ferroelectric  phase.  These  terms  can  lead  to  a qualitatively 
different  behavior  for  certain  of  the  elastic  and  dielectric 
properties  at  the  phase  transition,  namely,  such  terms  in  the 
free  energy  would  produce  a small  discontinuous  jump  at  the 
transition  temperatures  superimposed  on  the  temperature  depend- 
ence expected  from  a free  energy  without  such  terms.  Existing 
experimental  data  on  dielectric  and  elastic  properties  do  not 
reveal  such  effects  with  reasonable  magnitude. 

Nevertheless,  these  terms  must  be  present  in  the  free-energy 
however  small  their  magnitude  or  their  effect  on  dielectric  and 
elastic  properties  - to  properly  account  for  the  reduction  in 
symmetry  at  the  phase  transition. 
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(2) 


Levanyuk  and  Sannikov  Analysis 


The  first  interpretation  of  the  phase  transition  in  gadolinium 
molybdate  as  driven  by  a zone  boundary  phonon  was  offered  by 
Levanyuk  and  Sannikov  (Ref.  2-3).  The  authors  recognized  that 
the  observed  cell  doubling  at  the  transition  temperature  indicat- 
ed that  the  order  parameter  describing  the  transition  must  be 
other  than  the  polarization.  They  were  then  able  to  apply  their 
earlier  more  generalized  analysis  (Ref.  2-8)  of  ferroelectric 
phase  transitions  where  the  polarization  was  not  the  order 
parameter.  Without  carrying  out  a detailed  group  theoretical 
analysis  of  the  transition,  Levanyuk  and  Sannikov  hypothesized 
that  the  transition  was  driven  by  a two-component  order  para- 
meter and  wrote  down  a thermodynamic  potential  to  describe  the 
phase  transition  (in  their  notation): 

$ = | (n2  + s2)  + (n2  + c2)2  + ~y  (ns)2 

+ l (n2  + s2)3  + * p2  + | u2  + 

XPu  + anSP  + bncu  - PE  - au 

where  (n,  S)  = two  component  order  parameter,  P = z component  of 
the  polarization  and  u = the  x-y  component  of  the  strain.  Using 
this  hypothesized  potential,  Levanyuk  and  Sannikov  calculated  the 
expected  temperature  dependence  of  the  elastic  and  dielectric 
constants  and  found  them  to  be  in  agreement  with  the  available 
experimental  results. 

Levanyuk  and  Sannikov  did  note  that  there  are  many  phase  transi- 
tions between  the  known  point  groups  of  the  high  and  low 
symmetry  phases  which  are  exactly  described  by  just  the  thermo- 
dynamic potential  they  used.  However,  without  knowledge  of  the 
precise  space  groups  of  the  high  and  low  symmetry  phases  they 
were  unable  to  determine  a priori  the  precise  form  of  the 
potential. 
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The  form  of  the  potential  analyzed  by  Levanyuk  and  Sannikov  has 
the  substantial  advantage  that  analytic  solutions  are  readily 

I calculated.  Moreover,  as  we  shall  argue  in  Par.  2.2,  it  is 

plausible  to  use  just  such  a thermodynamic  potential  as  an 
approximation  to  the  complete  thermodynamic  potential  derived  by 
Dvorak . 

» 

( 3 ) Aizu  Analysis 

Aizu's  analysis  (Ref.  2-4)  of  the  phase  transition  in  gadolinium 
molybdate  is  based  on  knowledge  of  the  high  and  low  symmetry 
space  groups,  as  well  as  the  observed  cell  doubling  at  the  phase 
transition.  Aizu  correctly  concludes  that  the  phase  transition 
is  brought  about  by  an  instability  in  a doubly-degenerate 
zone-boundary  phonon.  The  method  used  to  deduce  the  form  of  the 
thermodynamic  potential  is  not  the  usual  calculation  involving 
the  irreducible  representations  of  the  high  symmetry  space 
group.  Rather  Aizu  enumerates  the  symmetry  transformations 
which  must  leave  the  potential  invariant  and  then  writes  down 
the  simplest  form  of  the  thermodynamic  potential  consistent  with 
these  symmetries.  In  order  to  simplify  the  computation  Aizu 
introduces  an  additional  hypothetical  symmetry,  while  carefully 
emphasizing:  "This  assumption  serves  to  simplify  the  form  of 

the  free-energy  function.  However,  it  is  not  essentially 
requisite  to  the  theory.  If  it  were  found  to  make  the  results 
of  the  theory  too  roughly  coincident  with  the  observations,  it 
would  be  withdrawn."  The  resulting  thermodynamic  potential  is 
equivalent  to  the  one  analyzed  by  Levanyuk  and  Sannikov.  Aizu 
demonstrates  that  the  predictions  of  his  analysis  are  in  good 
agreement  with  the  extant  elastic,  dielectric  and  optical 
measurements  made  on  gadolinium  molybdate. 

( 4 ) Dorner,  Axe  and  Shirane  Discussion 

In  reporting  their  neutron  scattering  measurements  on  terbium 
molybdate,  Dorner,  Axe  and  Shirane  (Ref.  2-6)  discuss  the 
theoretical  analysis  of  their  result.  The  major  contribution  of 
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this  work  is,  of  course,  the  direct  observation  of  the  zone- 
boundary soft  mode  which  causes  the  phase  transition.  However, 
the  authors  make  two  useful  observations  which  help  illuminate  the 
theory.  First,  they  include  the  complete  anistropic  sixth  order 
term  in  the  potential  for  the  soft  mode;  previous  workers  had  in- 
cluded only  an  isotropic  term  for  simplicity.  Second,  they  note 
that  the  essential  difficulty  in  carrying  out  an  analytic 
solution  of  the  complete  thermodynamic  potential  advanced  by 
Dvorak  is  that  there  is  a temperature  dependence  in  the  phase 
angle  4>  = arctan  (q,/q^).  To  circumvent  this  difficulty  the 
authors  append  the  hypothesis  that  this  angle  does  not  vary  very 
much  and  analyze  their  data  according  to  a thermodynamic  potential 
essentially  equivalent  to  that  of  Levanyuk  and  Sannikov  and 
Aizu.  The  experiment  is  found  to  be  in  excellent  agreement  with 
the  theoretical  expressions.  This  provides  post-hoc  justification 
for  their  hypothesis  of  constant  <J>.  (In  our  analysis  in  Par. 

2.2  the  hypothesis  of  constant  <J>  is  shown  to  be  quite  plausible.) 

( 5 ) Pytte  Analysis 

An  analysis  of  the  phase  transition  in  gadolinium  molybdate  is 
given  by  Pytte  (Ref.  2-5)  from  the  soft  mode  perspective.  Using 
the  known  space  group  of  the  high  symmetry  phase  and  the  cell 
doubling  observed  at  the  transition,  Pytte  identifies  the  mode 
whose  instability  causes  the  transition  as  one  of  three  zone 
boundary  modes.  Introducing  a phenomological  Hamiltonian  which 
describes  both  the  phase  transition  in  this  mode  and  the  soft 
mode  coupling  to  a polar  and  an  elastic  mode,  Pytte  is  able  to 
account  for  the  extant  elastic  and  ferroelectric  measurements 
on  gadolinium  molybdate. 

While  Pytte  observes  that  this  soft  mode  is  doubly  degenerate, 
the  Hamiltonian  involves  only  the  amplitude  for  this  soft  mode; 
no  angular  coordinate  is  introduced  which  would  allow  for  a 
temperature  variation  in  the  particular  linear  combination  of 
modes  which  constitutes  the  soft  mode.  Thus,  while  it  captures 
the  essential  points,  this  analysis  is  not  sufficiently  detailed 
to  give  an  analysis  comparable  to  that  of  Dvorak. 


2.1.2  Review  of  Measurements 


This  section  reviews  the  previous  pyroelectric,  dielectric  and 
thermal  measurements  made  on  rare-earth  molybdates. 

( 1 ) Pyroelectric  Coefficient  and  Spontaneous  Polarization 

Measurements 

Measurement  Techniques  and  Boundary  Conditions.  There  are  two 
essentially  different  approaches  to  determining  the  pyroelectric 
coefficient  of  a ferroelectric  material.  One  is  to  directly 
measure  the  pyroelectric  current  while  varying  the  temperature 
either  in  a dc  fashion  (Byer-Roundy  technique)  or  in  an  ac 
fashion  (Chynoweth  technique).  The  other  approach  is  to  measure 
the  spontaneous  polarization  as  a function  of  temperature  either 
by  analog  integration  of  the  pyroelectric  current  during  dc 
heating,  or  from  a ferroelectric  hysteresis  loop  at  a sequence 
of  temperatures  (Sawyer-Tower  bridge  technique).  The  resulting 
curve  can  then  be  numerically  differentiated.  All  these  tech- 
niques have  inherent  dangers.  At  low  dc  heating  rates  the 
pyroelectric  current  is  small  and  hard  to  measure  accurately; 
at  higher  heating  rates  it  is  difficult  to  maintain  thermal 
equilibrium  in  the  sample.  The  ac  heating  technique  requires  a 
normalization  and  an  independent  measurement  of  the  heat  capacity 
of  the  sample.  Tracing  ferroelectric  hysteresis  loops  introduc- 
es many  complexities  involved  in  the  domain  switching;  numerical- 
ly differentiating  experimental  data  requires  even  more  courage. 
It  is  not  surprising  then  that  there  are  few  precision  measure- 
ments of  the  pyroelectric  coefficient  of  the  rare-earth  moly- 
bdates, even  though  there  are  fairly  many  measurements  of 
spontaneous  polarization. 

The  mechanical  boundary  conditions  for  measurements  of  the  pyro- 
electric current  and  spontaneous  polarization  are  generally  good 
approximations  to  constant  (zero)  stress.  The  electrical 
boundary  conditions  are  generally  constant  applied  field  of  zero 
magnitude,  but  some  pyroelectric  current  measurement  use  a 
finite  field  to  maintain  the  sample  in  a poled  state  especially 
during  a cooling  part  of  the  cycle. 
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Direct  Measurements  of  Pyroelectric  Coefficient.  The  literature 
contains  only  three  efforts  to  directly  measure  pyroelectric 
currents.  The  two  efforts  to  measure  the  pyroelectric  current 
from  GMO  were  plagued  by  an  experimental  artifact  that  we  shall 
call  multiple  peaking.  The  essential  characteristic  of  this 
artifact  is  that  in  the  vicinity  of  the  transition  temperature 
the  pyroelectric  current  exhibits  a number  of  spikes  which  are, 
in  a given  sample,  reproducible  to  a fair  degree.  These  spikes 
correspond  to  different  portions  of  the  sample  going  through  the 
f erroelectric-paraelectric  phase  transition  over  a finite  range 
of  temperatures. 

Rabinovich  (Ref.  2-9)  observed  the  multiple  peaking  in  gadolinium 
molybdate  and  attributed  it  to  the  motion  of  domain  walls  caused 
by  inhomogenities  ("inclusions,  residual  stress,  etc.")  in  the 
single  crystal  sample.  The  only  firm  conclusion  Rabinovich 
draws  from  his  measurements  is  that  the  change  in  the  polari- 
zation in  the  region  of  the  phase  transition  amounts  to 
4 x 10  ® C/cm^. 

A second  effort  to  directly  measure  the  pyroelectric  current 

from  GMO  is  reported  in  three  papers:  Ullman,  Ganguly  and  Hardy 

(Ref.  2-10),  Ullman,  Ganguly  and  Zeidler  (Ref.  2-11)  and  Ullman, 

Cheung,  Rakes  and  Ganguly  (Ref.  2-12).  These  works  ascribe  the 

multiple  peaking  to  a thermal  gradient  near  the  surface  of  the 

sample  which  leads,  through  differential  thermal  expansion,  to  a 

stress-induced  polarization  switching.  The  data  taken  far  from 

the  transition  temperature  should  be  free  of  this  effect;  values 

-9  2 

for  the  pyroelectric  coefficient  of  1.4  x 10  C/cm  K at 
143  C and  2.5  x 10  ^ C/cm^  K at  151°C  are  given.  Their 
value  at  151  C is  just  9%  below  the  value  which  can  be  deduced 
from  the  data  of  Sawaguchi  and  Cross  (Ref.  2-13)  described 
below . 

The  other  direct  measurement  of  the  pyroelectric  current  was 
performed  on  TMO  by  Keve,  Abrahams,  Nassau  and  Glass  (Ref. 

2-14).  In  this  case  inhomogeneities  in  the  sample  probably  so 
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broadened  the  phase  transition  region  that  it  was  not  possible 
to  resolve  any  discontinuities  at  the  transition  temperature 
l (Ref.  2-15).  In  the  region  where  the  pyroelectric  coefficient 

does  not  vary  rapidly  with  temperature,  this  data  should  be 
quite  reliable.  Unfortunately  no  analysis  of  the  data  is 
presented . 

Temperature  Dependence  of  Spontaneous  Polarization,  GMO.  Of 
the  several  papers  which  describe  the  temperature  dependence  of 
the  spontaneous  polarization  of  rare-earth  molybdates,  the  most 
comprehensive  study  was  performed  by  Sawaguchi  and  Cross  (Ref. 
2-13)  on  GMO  from  4.2  K to  473  K with  particular  care  near  the 
phase  transition.  Low  frequency  (0.01  Hz)  hysteresis  loops  were 
used  to  obtain  the  spontaneous  polarization  down  to  75  K;  below 
this  temperature,  the  charge  integration  technique  was  used 
because  the  coercive  field  is  too  large  to  obtain  saturated 
loops.  The  resulting  data  were  subjected  to  a detailed  analysis 
which  included  a fit  to  the  form  expected  from  Aizu's  (Ref.  2-16) 
analysis  of  improper  ferroelectrics , namely, 

P(T)  = [1.78  (433  - T)1/2  + 2.59]  x 10_6  C/cm2,  T < 432.3  K. 

The  fit  is  good  only  within  5 K of  the  transition.  A significant- 
ly better  fit  is  obtained  to  the  expression, 

P ( T)  = 3.68  (433.6  - T)1/3  x 10-8  C/cm2,  T < 432.3  K. 

A unique  contribution  of  this  work  is  the  extrapolated  value  of 

the  spontaneous  polarization  at  absolute  zero:  2.86  + 0.1  x 

-7  2 -7 

10  C/cm  . The  value  measured  at  20°C,  2.01  + 0.15  x 10 

2 

C/cm  , and  the  discontinuity  at  the  transition  temperature  of 

— 8 2 

432.3  K,  4.0  + 0.2  x 10  C/cm  , are  also  useful  for  com- 
parisons with  other  works. 

Only  one  (Ref.  2-17)  of  the  other  measurements  using  the  loop 
tracing  technique  on  GMO  was  able  to  distinguish  the  discontinu- 
ous jump  to  zero  from  a smooth  decrease  to  zero  of  the  spontaneous 
polarization  at  the  transition  temperature.  Evidence  for  a 
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first  order  transition  is  reported  by  Kumada  (Refs.  2-18,  2-19) 
who  observed  double  hysteresis  loops  above  the  transition  temper- 
ature Sawaguchi  and  Cross  (Ref.  2-20)  expected  to  observe  such 
double  loops  but  were  unable  to  find  the  phenomenon. 

Cummins  (Ref.  2-17)  also  finds  a discontinuous  jump  of  about  4.2  x 
_o  2 

10  C/cm  at  the  159°C  transition  temperature.  His  measure- 
ments on  GMO  were  made  at  "low  frequency  (hand  traced)"  and  span 

from  25°C  through  the  transition  temperature.  Cummins'  value 

-7  2 

for  the  20°C  polarization,  2.0  x 10  C/cm  , also  agrees  with 
the  measurements  of  Sawaguchi  and  Cross  (Ref.  2-13). 

In  the  first  report  of  ferroelectricity  in  the  rare-earth 

molybdates,  Borchardt  and  Bierstedt  (Ref.  2-21)  measured  the 

ferroelectric  hysteresis  loops  of  GMO  at  60  Hz  from  25°C  through 

the  phase  transition  at  159°C.  The  data  has  a scatter  of  about 

10%  around  a smooth  trend  and  no  conclusion  about  the  order  of 

the  transition  can  be  drawn.  The  room  temperature  value  of  the 

-7  2 

spontaneous  polarization  is  given  as  1.7  x 10  C/cm  . This 
is  the  first  report  of  ferroelectricity  in  rare-earth  molybdates. 

Other  early  measurements  on  the  ferroelectric  hysteresis  of  GMO 

are  reported  in  a short  note  by  Aizu,  Kumada,  Yumoto  and  Ashida 

(Ref.  2-22)  and  further  elaborated  on  by  Kumada,  Yumoto  and 

Ashida  (Ref.  2-23).  The  data  is  taken  at  1 MHz  from  25°C 

through  the  phase  transition.  There  is  little  scatter  about  a 

smooth  trend  which  goes  continuously  to  zero  at  the  transition 

-7  2 

temperature.  The  value  of  1.8  x 10  C/cm  is  reported  for 
the  polarization  at  20°C. 

Temperature  Dependence  of  Spontaneous  Polarization,  TMO  and 
Other  Molybdates.  Terbium  molybdate  was  examined  by  Keve, 
Abrahams,  Nassau  and  Glass  (Ref.  2-14)  using  polarization  re- 
versal, charge  integration  and  pyroelectric  measurements.  A value 

-7  2 

of  1.90  + 0.05  10  C/cm  is  obtained  for  the  polarization  at 
25°C,  this  agrees  with  the  earlier  measurement  by  Borchardt  and 
Bierstedt  (Ref.  2-24)  noted  below.  The  report  that  the  polari- 
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zation  decreases  continuously  to  zero  at  the  phase  transition  is 
probably  due  to  inhomogeneities  broadening  the  transition  region 
(Ref.  2-15). 

In  a survey  of  twelve  ferroelectric  rare-earth  molybdates, 

Borchardt  and  Bierstedt  (Ref.  2-24)  give  spontaneous  polarization 

at  two  or  three  selected  temperatures  for  each  material.  These 

data  were  probably  taken  from  60  Hz  hysteresis  loops  as  was 

described  in  the  authors'  earlier  study  (Ref.  2-21)  of  GMO . 

Values  for  the  polarization  at  25°C  vary  from  1.4  to  2.4  x 

10  7 C/cm2;  the  value  reported  for  GMO,  1.7  x 107  C/cm2,  is 

precisely  the  same  as  their  earlier  measurement.  The  25°C 

-7  2 

spontaneous  polarization  of  TMO  is  given  as  1.8  x 10  C/cm  . 

( 2 ) Dielectric  Permittivity 

Thermal  and  Mechanical  Boundary  Conditions.  Conceptually,  four 
distinct  dielectric  permittivities  can  be  distinguished  depend- 
ing on  the  thermal  (adiabatic  or  isothermal)  and  mechanical 
(constant  stress  or  constant  strain)  boundary  conditions.  In 
practice,  dielectric  measurements  are  made  with  a small  ac  probe 
field  with  no  external  mechanical  stress.  Which  dielectric 
constant  is  measured  in  a given  experiment  depends  on  the 
relationship  of  the  frequency  of  the  probe  field  to  the  thermal 
and  mechanical  time  constants  of  the  sample  under  study.  These 
thermal  and  mechanical  time  constants  depend  not  just  upon  the 
properties  of  the  materials  but  also  upon  the  dimensions  of  the 
sample  and  the  experimental  environment. 

Since  typical  thermal  time  constants  are  on  the  order  of  magni- 
tude of  seconds,  one  may  assume  that  the  typical  low-frequency 
experiment  which  uses  a probe  field  with  frequencies  near  a 
kilohertz  measures  an  adiabatic  dielectric  permittivity.  The 
relationship  between  the  adiabatic  and  isothermal  dielectric 
susceptibilities  (e  = 1 + X)  is  given  by: 


xT/xs 


1 = Tp2/cEXS  = Tp2/cPXT, 


33 


m 4 


where  p is  the  pyroelectric  coefficient,  T is  the  absolute 

E P 

temperature  and  c and  c are  the  specific  heats  at  constant 
field  and  polarization  respectively  (the  appropriate  mechanical 
boundary  condition  being  the  same  for  every  quantity) . Thus  the 
difference  between  them  becomes  significant  only  just  below  the 
phase  transition  where  the  pyroelectric  coefficient  becomes 
large . 

With  regard  to  mechanical  boundary  conditions,  the  value  obtain- 
ed at  low  frequencies  (typically  kHz)  corresponds  to  the 
adiabatic  dielectric  permittivity  at  constant  (zero)  stress.  It 
is  only  at  very  high  frequencies  (typically  many  MHz)  that 
the  crystal  cannot  equilibrate  mechanically  as  fast  as  the  field 
changes  so  that  the  resulting  value  measured  is  the  adiabatic 
dielectric  permittivity  at  constant  (zero)  strain.  The  permit- 
tivity at  constant  stress  is  higher  than  the  value  at  constant 
strain,  because,  at  constant  stress,  the  applied  field  can 
electrostrictively  cause  a distortion  in  the  lattice  which 
enhances  piezoelectr ically  the  polarization  already  caused  by 
the  direct  effect  of  the  applied  field. 

Adiabatic,  Mechanically  Clamped  Dielectric 'Constant  Measurements. 
For  the  rare-earth  molybdates,  two  measurements  of  the  adiabatic 
dielectric  permittivity  at  constant  (zero)  strain  appear  in  the 
literature.  Both  were  performed  at  19  MHz  and  give  values  which 
are  essentially  unchanged  even  as  the  temperature  is  scanned 
through  the  paraelectric-f erroelectric  phase  transition.  Cross, 
Fouskova  and  Cummins  (Ref.  2-25)  report  a value  of  about  9.5  for 
GMO  in  the  range  from  25°C  to  180°C.  Keve,  Abrahams,  Nassau 
and  Glass  (Ref.  2-14)  obtain  a value  of  9.2  + 0.5  for  TMO  in  the 
range  from  25°C  to  170°C.>  These  observations  on  the  mechanically 
clampdd  dielectric  constant  were  the  first  indication  that  the 
rare-earth  molybdates  are  truly  special  ferroelectrics . The 
correct  understanding  of  this  point  is  essential  to  their 
description. 
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Adiabatic,  Mechanically  Free  Dielectric  Constant  Measurements. 

All  the  reported  low-frequency  measurements  on  the  rare-earth 
molybdates  give  essentially  the  same  qualitative  behavior  for 
the  temperature  dependence  of  the  adiabatic  zero-stress  dielec- 
tric constant.  Namely,  from  a magnitude  of  about  10  at  room 
temperature  the  permittivity  increases  in  the  vicinity  of 
transition  temperature  to  a peak  value  about  15-20%  larger  and 
then  abruptly,  at  the  transition  temperature,  drops  back  to 
essentially  its  room  temperature  value  and  remains  constant  in 
the  paraelectric  phase.  Most  measurements  exhibit  a small 
(1-2%)  tailing  which  extends  a few  degrees  above  the  transition; 
none  of  the  papers  offer  an  opinion  as  to  whether  this  tail  is 
an  intrinsic  property  of  the  material  or  an  experimental  artifact. 

The  first  measurement  of  the  low  frequency  dielectric  permit- 
tivity of  gadolinium  molybdate  was  given  by  Cross,  Fouskova  and 
Cummins  (Ref.  2-25).  They  used  a probe  field  of  20  V/cm  at  1 kHz 
to  scan  from  25°C  to  180°C.  The  background  value  of  the  dielec- 
tric constant  is  reported  as  10.0  and  the  maximum  value  of  11.6 
occurs  at  the  transition  temperature,  159  C. 

These  values  are  confirmed  by  Fousek  and  Konak  (Ref.  2-26)  who 
focus  attention  on  an  essential  experimental  point  for  all 
low-frequency  dielectric  measurements,  namely  the  requirement 
that  single  domain  samples  be  used  so  that  domain  wall  motion 
does  not  spuriously  enhance  the  measured  value  in  the  ferro- 
electric phase.  This  is  a particularly  important  consideration 
for  measurements  made  near  the  transition  temperature.  Fousek 
and  Konak  report  measurements  on  single-domain  GMO  (poled  by 
application  of  a 20  kV/cm  field  prior  to  measurement)  in  the 
temperature  range  100°C  to  190°C  at  a frequency  of  1 kHz.  They 
obtain  a background  value  of  10.0  for  the  dielectric  permit- 
tivity which  rises  to  a peak  of  11.6  at  the  159°C  transition 
temperature.  However,  the  same  measurements  on  a multi-domain 
sample  give  a peak  value  of  12.9;  moreover,  the  permittivity  is 
consistently  higher  than  the  single  domain  sample  all  the  way 
down  to  100°C.  No  differences  between  the  two  samples  appear  in 
the  paraelectric  phase,  as  expected. 
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Cummins  (Ref.  2-17)  summarizes  verbally  the  result  of  his 
measurements  on  GMO  at  1 kHz;  the  data  is  not  displayed  graph- 
ically. From  a value  of  10  at  25°C,  he  reports,  the  permit- 
tivity increases  by  about  10%  through  the  transition  at  159°C. 

At  variance  with  the  above  results  is  the  measurements  of  Smith 
and  Burns  (Ref.  2-27).  GMO  was  studied  from  25°C  to  170°C  at  a 
frequency  of  1 kHz.  A background  value  of  8.0  and  a peak  value 
of  10.7  were  reported.  These  values  were  probably  not  determined 
with  high  precision  as  the  authors  only  wanted  the  results  to 
argue  that  the  dielectric  permittivity  "is  essentially  temper- 
ature-independent". The  prinicpal  goal  of  their  work  was  the 
analysis  of  the  electro-optic  coefficients,  so  we  do  not  regard 
this  discrepancy  as  significant. 

The  only  complete  temperature  dependence  for  the  low  frequency 
permittivity  of  terbium  molybdate  is  reported  by  Keve,  Abrahams, 
Nassau  and  Glass  (Ref.  2-14)  at  400  Hz  from  25°C  to  170°C.  They 
reported  that  the  background  value  of  10.0  +0.5  rises  to  a peak 
of  about  15  at  the  transition  temperature  of  163  + 1°C. 

Dielectric  Constant  at  Room  Temperature.  Besides  the  cited 
measurements,  there  have  been  some  measurements  simply  of  the 
room  temperature  value  of  the  low-frequency  dielectric  permit- 
tivity. Aizu,  Kumada,  Yumoto  and  Ashida  (Ref.  2-22)  obtain  at 
400  Hz  a value  of  10.5  for  GMO  at  25°C.  Borchardt  and  Bierstedt 
(Ref.  2-24)  give  dielectric  constants  for  twelve  ferroelectric 
rare-earth  molybdates.  They  do  not  describe  the  details  of 
their  measurements  but  we  may  plausibly  expect  that  they  are 
room  temperature  values  obtained  at  a low  frequency,  since  in  an 
earlier  work  (Ref.  2-21)  the  same  authors  obtain  a value  of 
about  10  for  the  dielectric  constant  of  GMO  using  a 60  Hz 
ferroelectric  loop  tracer.  The  reported  values  range  from  a low 
of  9.5  for  europium  molybdate  to  a high  of  13  for  europium- 
gadolinium  and  europium-terbium  mixed  crystals.  The  values  of 
10  for  GMO  and  11  for  TMO  are  consistent  with  the  results  of 
other  workers. 
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Pressure  Dependence  of  Dielectric  Constant.  Shirokov,  Mylov, 
Baranov  and  Prokhortseva  (Ref.  2-28)  performed  an  interest- 
I ing  study  of  the  effect  of  hydrostatic  pressure  on  the  temper- 

ature dependence  of  the  low-frequency  dielectric  permittivity  of 
GMO.  They  measure  at  10  kHz  for  hydrostatic  pressure  up  to  2.6 
kbar  and  obtain  a background  value  of  10  and  a peak  value  of  11 
for  the  permittivity.  The  shape  of  the  temperature  dependence 
does  not  change  with  increasing  pressure  except  for  an  overall 
shift  of  the  entire  curve  to  higher  temperatures.  This  shift 
varies  linearly  with  hydrostatic  pressure  and  nas  a magnitude  of 
29.5  + 0.3  K/kbar.  A similar  study  on  Sm2(MoC>4)0  by 
Midorikowa,  Ishibashi  nd  Takagi  (Ref.  2-29)  finds  a temperature 
shift  of  29.4  K/kbar. 

In  summary  then,  all  of  the  measurements  on  the  rare-earth 
molybdates  are  in  reasonable  agreement  with  each  other.  More- 
over they  correspond  qualitatively  to  what  is  expected  from  the 
theory.  However,  the  literature  contains  no  detailed  quanti- 
tative analysis  of  the  temperature  dependence  of  the  dielectric 
permittivity  of  rare-earth  molybdates. 

( 3 ) Specific  Heat  Measurements 

Three  measurements  of  the  specific  heat  of  rare-earth  molybdates 
have  been  reported.  The  measurements  of  Fouskova  (Ref.  2-30)  and 
Cheung  and  Ullman  (Ref.  2-31)  were  made  on  gadolinium  molybdate, 
while  those  of  Keve,  Abrahams,  Nassau  and  Glass  (Ref.  2-14)  were 
made  on  terbium  molybdate.  Only  Fouskova  reports  absolute 
measurements;  the  other  papers  calibrate  their  relative  measure- 
ments by  her  absolute  determination. 

Fouskova  (Ref.  2-30)  reports  measurements  on  gadolinium  molybdate 
made  with  a Perkin-Elmer  differential  scanning  calorimeter  in 
the  temperature  100°C  to  200°C.  From  a value  of  about  0.104 
cal/gK  at  100°C,  the  specific  heat  gradually  rise  about  15%  as 
the  phase  transition  is  approached  and  then  shows  a sharp  peak 
at  the  phase  transition  at  about  160°C.  Above  the  transition 
the  specific  heat  drops  abruptly  to  0.110  cal/gK  and  remains 
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essentially  constant  out  to  200°C.  Sample-to-sample  variations 
of  1 K in  the  transition  temperature  are  reported.  The  transi- 
tion entropy  is  estimated  to  be  0.22  + 0.08  cal/mol  K and  the 


transition  heat  (in  a 0.5  K neighborhood  of  the  transition)  is 
estimated  to  be  0.048  cal/g  + 10%.  No  detailed  analysis  of  the 
temperature  dependence  of  the  specific  heat  is  reported.  The 
author  estimates  a precision  of  + 5%  in  the  absolute  determina- 
tion; the  data  exhibit  a scatter  of  less  than  1%  about  a 


smooth  trend.  The  samples  were  silver  plated  on  all  sides  to 
insure  the  electrical  boundary  condition  of  E = 0;  mechanical 
boundary  conditions  are  not  mentioned. 


Cheung  and  Ullman  (Ref.  2-31)  performed  measurements  on  gadolin- 
ium molybdate  in  the  temperature  range  from  50°C  to  240°C  using 
an  ac  heating  technique.  This  technique  produces  a relative 
measurement  whose  scale  is  calibrated  by  setting  the  value  of 
the  specific  heat  just  above  the  transition  to  the  value  report- 
ed by  Fouskova  (Ref.  2-30).  The  measurements  are  reproducible 
to  within  0.1  K of  the  transition  and  no  hysteresis  effects 
appear  between  heating  and  cooling.  The  qualitative  features  of 
the  specific  heat  temperature  dependence  are  the  same  as  report- 
ed by  Fouskova;  however,  Cheung  and  Ullman  carry  out  a much  more 
detailed  analysis  of  the  data.  The  data  fit  well  to  linear 
temperature  dependence  both  far  below; 

cp  = 2.43  x 10"4  T - 6.28  x 10~4  cal/gK,  50°C  < T < 110°C 
and  above; 

c = 3.65  x 10-5  T + 0.10  cal/gK,  159°C  < T < 240°C 
P 

the  transition  temperature,  although  the  two  linear  regions  do 
not  join  together  smoothly.  The  data  just  below  the  transition 
were  fitted  to  a number  of  different  expressions.  A power  law 
divergence  at  the  transition  temperature  gave  an  exponent  of  2.1 
but  fit  the  data  very  poorly  within  0.5  K of  the  transition. 

Fits  to  a logarithmic  divergence  at  the  transition  were  perform- 
ed in  two  temperature  ranges  below  the  transition  with  the 
following  results; 


Cp  = 0.09  - 0.0059  In [ ( 4 32  - T)/432]  cal/gK,110°C  < T < 157. 4°C 


cp  = 0.0787  - 0.0077  ln[(432  - T)/432]  cal/gK, 1 56 . 5° C < T < 158. 9°C 

The  heat  and  entropy  of  transition  are  given  as  0.108  cal/mol  and 
0.25  cal/mol  K respectively.  No  quantitative  estimate  of  the 
accuracy  of  the  fits  is  presented,  but  the  graphical  display  of 
data  points  versus  fitted  curve  is  impressive.  Estimates  of  the 
accuracy  of  the  experimental  data  is  not  made  explicitly. 

However  the  authors  convey  the  impression  that  they  consider  the 
data  accurate  to  only  a few  percent,  since  they  regard  a 2% 
difference  between  cp  and  cv  as  negligible  and  employ  an 
approximation  good  to  1%  to  calculate  the  transition  entropy. 

The  authors  state  that  care  has  been  exercised  to  eliminate 
thermal  gradients  and  mechanical  stresses,  but  do  not  comment  on 
electrical  boundary  conditions. 

Keve,  Abrahams,  Nassau  and  Glass  (Ref.  2-14)  report  their 
measurements  of  the  specific  heat  of  terbium  molybdate  in 
arbitrary  units.  Quantitative  estimates  of  the  transition  heat 
38  cal/mol)  and  transition  entropy  (0.23  cal/mol  K)  are  obtained 
by  setting  the  magnitude  of  the  specific  heat  of  TMO  at  100°C  to 
the  same  value  that  Fouskova  (Ref.  2-30)  reported  for  GMO  at 
that  temperature.  The  qualitative  temperature  dependence  of  the 
specific  heat  differs  from  the  reports  of  Fouskova  (Ref.  2-30) 
and  Cheung  and  Ullman  (Ref.  2-31)  in  that  the  sharp  peak  in  the 
specific  heat  at  the  transition  temperature  has  an  appreciable 
tail  extending  some  5 to  10  K into  the  high  temperature  phase. 

We  regard  the  quantitative  estimates  as  plausible  and  ascribe 
the  qualitative  discrepancy  to  an  experimental  artifact.  Such 
an  assessement  merits  explanation.  First,  the  high  temperature 
tail  on  the  specific  heat  is  not  due  to  any  gross  inhomogeneities 
in  the  sample  since  there  is  no  similar  tailing  in  the  pyro- 
electric coefficient  (nor  the  dielectric  constant  and  birefring- 
ence). There  probably  are  inhomogeneities  rounding  out  the 
measured  values  in  the  vicinity  of  the  phase  transition  but 
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these  probably  extend  over  a range  of  1 or  2 K.  Second,  just 
such  a tailing  effect  can  be  produced  by  the  inertia  of  the 
measuring  apparatus  during  a scan  with  increasing  temperature 
(measuring  specific  heat  while  heating  is  the  standard  experi- 
mental procedure).  Thus  this  entire  high  temperature  tail 
- should  be  associated  with  the  peak  at  the  phase  transition. 

From  this  interpretation  it  is  also  clear  that  the  heat  and 
entropy  of  transition  can  give  reasonable  values  even  though 
some  of  the  contribution  to  them  is  shifted  a few  degrees 
towards  higher  temperatures.  No  estimate  of  the  accuracy  of  the 
measurements  is  given.  Neither  the  elastic  nor  the  electrical 
boundary  conditions  ar  mentioned. 

For  the  specific  heats  of  rare-earth  molybdates  then,  the 
reported  measurements  are  in  general  agreement  and  have  qualita- 
tively the  expected  behavior.  No  attempt  has  been  made  to  fit 
the  measurements  to  the  temperature  dependence  expected  for  an 
improper  ferroelectric. 


2.2 


Theory 


In  this  section  the  complete  thermodynamic  potential  which 
describes  the  phase  transition  in  the  rare-earth  molybdates  is 
examined  in  detail.  After  reducing  out  the  mechanical  strains 
and  introducing  some  reasonable  approximations,  expressions  for 
the  temperature  dependence  of  the  pyroelectric  coefficient, 
dielectric  constant  and  specific  heat  are  derived.  These 
expressions  are  subsequently  compared  with  experimental  results 
and  found  to  be  in  excellent  agreement. 

2.2.1  Complete  Thermodynamic  Potential 

The  experimental  determination  of  the  space  groups  of  both  the 
paraelectric  phase  D2d^P^lm)  and  the  ferroelectric  phase 
C2v(Pba2)  in  gadolinium  molydbate  allow  the  unequivocal  de- 
termination of  the  form  of  the  thermodynamic  potential  associat- 
ed with  the  ferroelectric  phase  transition  in  the  molybdates: 


6(qi,q2,P,u)  = *o(T)  + ^±1  (q2  + q2j  + (2qiq2)2  + (q2  - q2) 2 

+ T <2qlq2>  (ql  " q2>  + T (2qlq2)2(qi  +q2> 

"^2  2 222  2 ^3  2 22  2 

+ 6 (<31  " V (ql  + q2)  + 6 (2qlq2)  (ql  ’ q2)  (ql  + q2) 

, 1 -1  _2  1 „ 2 1 2 2, 

+ -r-  Y P + — C.  -U  + - C,.  u + u ) 

2 z 2 66  xy  2 11  xx  yy 

12  , 

+ 2 C33uzz  + C12UxxUyy  + C13  Uxx  + uyy  uzz 

+ SlPz(2qlq2)  + SPz  {ql  ~ q2}  + dPzUxy 

+ blUxy(2qiq2)  + b2Uxy (qi  “ q2) 

+ b4Uzz(qi  + q2)  + b3(Uxx  + uyy’ + q2) 


-PE-au  -au  -au  -au. 
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where,  (qlfq2)  = two  component  order  parameter  describing  the 
soft  mode,  P = component  of  the  polarization  along  the  ferro- 
electric z axis,  Ez  = component  of  the  electric  field  along  the 
z axis,  = components  of  the  strain  tensor,  and  = compon- 

ents of  the  stress  tensor.  The  coefficients  8^,  82'  63,  Yj 
Y-,,  Y3,  X d,  a^,  a2/  b^,  b2»  b^,  and  b^  are  all  taken 


1' 


'1'  p2 ' 

2 r I 31  A 1 u , ^ j , 01 3 f a2  ' u I ' ^ ^ ^ 

to  be  temperature  independent,  while  the  coefficient  a = 
ao(T-TQ)  is  taken  to  have  a linear  temperature  dependence. 

The  thermodynamic  potential  is  written  in  this  form  to  emphasize 
the  fact  that  the  polarization  and  strain  couple  to  two  independ- 
ent quantities  which  are  quadratic  in  the  order  parameter,  to 
2 2 

wit  2q^2  and  q^  - q2* 


( 1 ) Physical  Interpretation  of  Terms  in  Thermodynamic  Potential 

The  physical  meaning  of  the  various  terms  becomes  apparent  when 
one  considers  some  special  cases.  The  terms  with  coefficients 
a(T),  B^,  B2f  §3,  y3,  y2*  Y3/  describe  the  temperature 
dependence  of  the  soft  mode  if  the  other  displacements  (strain 
u^j  and  polarization  P)  were  held  constant.  These  terms  would 
then  describe  a phase  transition  for  the  soft  mode  which  is 
either  first  or  second  order  depending  on  the  sign  of  the  B's. 

The  term  with  coefficient  x~^  describes  the  behavior  of  the 
polarization  if  the  other  displacements  (strain  u^  and  soft 
mode  q^,q2)  were  held  constant.  In  such  a case  the  system 
would  have  a temperature  independent  electric  susceptibility  of  x* 
Similarly  the  terms  with  coefficients  describe  the  behavior 

of  the  various  strain  components  if  the  other  displacements 
(polarization,  P,  and  soft  mode  amplitude,  q3,q2)  were  held 
constant.  In  such  a case  these  terms  would  yield  temperature 
independent  elastic  constants. 

In  most  realizable  experiments,  it  is  not  the  displacements  but 
rather  the  conjugate  forces  which  are  held  constant.  The  force 
conjugate  to  the  polarization  is  just  the  applied  electric  field 
and  the  forces  conjugate  to  the  strain  components  are  just  the 
corresponding  stresses.  However,  the  force  conjugate  to  the 


soft  mode  amplitude  usually  goes  unnamed  since  there  is  no  way 
to  realize  experimentally  any  value  other  than  zero  for  this 
force.  In  describing  the  usual  experiments  performed  with 
constant  forces,  one  must  take  into  account  the  couplings  amongst 
the  displacements.  The  terms  with  coefficient  a^  and  a2 
describe  the  direct  interaction  of  the  soft  mode  with  the 
polarization,  while  the  terms  with  coefficients  b^,  b2  b^  and 
b4  describe  the  coupling  of  the  soft  mode  and  the  various 
strain  components.  These  terms  are  the  lowest  order  interactions 
allowed  by  crystallographic  symmetry;  thev  are  quadratic  in  the 
soft  mode  amplitude  and  linear  in  the  polarization  and  strains. 

As  the  soft  mode  goes  through  its  phase  transition,  these  terms 
give  rise  to  a spontaneous  polarization  and  a spontaneous  elastic 
distortion  as  well  as  introducing  a temperature  dependence  in 
the  dielectric  susceptibility  and  elastic  constants.  Finally 
the  usual  bilinear  piezoelectric  coupling  between  the  polarization 
P and  the  strain  u is  given  by  the  term  with  coefficient  d. 

Since  the  phase  transition  in  the  molybdates  is  caused  by  an 
instability  in  the  soft  mode,  it  is  usual  to  regard  the  inter- 
action terms  as  describing  how  the  phase  transition  in  the  soft 
mode  transmits  its  effect  to  the  polarization  and  strains  thereby 
causing  a spontaneous  polarization  and  elastic  distortion. 

However,  these  coupling  terms  also  mediate  the  effect  of  the 
polarization  and  strain  (or  their  conjugate  forces:  electric 
field  and  stress)  back  on  the  soft  mode  itself. 

( 2 ) Estimates  of  Magnitudes 

It  is  difficult  to  make  general  statements  about  the  relative 
magnitudes  of  the  various  terms,  since  the  usual  experiments  per- 
formed with  constant  forces  give  information  only  about  various 
combinations  of  the  terms.  However,  there  is  one  experiment,  viz., 
the  high  frequency  dielectric  permittivity,  which  corresponds  to 
a constant  displacement  (strain).  This  enables  one  to  directly 
assess  the  magnitude  of  the  coefficients  a^,  and  a2*  The 
observation  that  the  high  frequency  dielectric  constant  for  the 
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rare-earth  molybdates  does  not  vary  with  temperature  means  that 

when  the  strain  is  held  constant  the  polarization  is  essentially 

decoupled  from  the  soft  mode,  so  that  ^ = &2  = 0.  Since  the 

polarization  has  only  very  weak  direct  coupling  to  the  soft 

mode,  the  observed  spontaneous  polarization  and  temperature 

dependent  dielectric  permittivity  must  arise  indirectly  from  the 

soft  mode  coupling  to  the  strain  u which  in  turn  couples 

xy 

piezoelectrically  to  the  polarization. 

One  aspect  of  the  interaction  of  the  soft  mode  and  the  strain 
can  be  learned  from  the  observed  large  presure  dependence  of  the 
transition  temperature  (Ref.  2-28).  This  result  implies  that  the 
isotropic  strains,  uxx,  uyy  and  uzz,  are  fairly  strongly  coupled 
to  the  soft  modes. 

2.2.2  Reduced  Thermodynamic  Potential 

Since  our  experiments  are  performed  under  stress-free  conditions, 
we  need  only  a reduced  thermodynamic  potential  describing  the 
polarization  and  the  soft  mode.  Using  the  zero  stress  conditions, 
= 0 while  = 0,  allows  one  to  express  the  strains 

in  terms  of  the  polarization  and  soft  mode  amplitude: 

uxy  = " C66  fbi(2c3i'32)  + b2  (ql  ~ q2)  3 


uxx  = - 

[<sn+s 

12 } b3  + 

S13b4] 

(qj  + 

q2 ) 

uyy  = 

[<su  + 

S12)  b3  + 

S13b4 

] (q i 

+ q2> 

uzz  = " 

[2s13b3 

+ S33b4] 

(qj  + 

q^) 

where,  s^  is  the  inverse  of  the  matrix  c^. 
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Inserting  these  into  the  complete  thermodynamic  potential  yields 
the  reduced  thermodynamic  potential  corresponding  to  conditions 
of  zero  stress. 


<J,(q1fq2'p>T) 


l 4.  a (T)  / 2 2.  1 ? ^9  9 9 9 

°(T)  + 1“  (ql  + q2)  + ~4  (2qlq2)  + T (qi  ~ q2> 

+ (2qiq2) (qi  - q2}  + t (2qiq2)2(qi  + q2) 

+ T (qi2  - q2)2  <qi  + q2>  + T (q2  - q22)  (ql2  + q2> 
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+ 2 P + alP(2qlq2)  + a2P(qi  " V " PE 


where , 
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It  is  important  to  remember  that  the  interaction  terms  between 
the  polarization  and  the  soft  mode  amplitude  give  the  effective 
coupling  under  stress-free  conditions.  Since  we  know  that  the 
direct  interaction  between  the  polarization  and  the  soft  mode  is 
essentially  negligible,  the  physical  origin  of  this  coupling 
lies  in  the  piezoelectric  coupling  of  the  polarization  to  the 
strain  which  itself  couples  directly  to  the  soft  mode  amplitude. 


I 
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2.2.3  Approximation  to  Reduced  Thermodynamic  Potential 


Even  the  reduced  thermodynamic  potential  does  not  yield  simple 
analytic  expressions  to  compare  with  the  measured  temperature 
dependence  of  the  polarization,  dielectric  permittivity  and 
specific  heat.  However,  as  noted  in  Par.  2.1  a variety  of 
apparently  different  approximations  do  yield  analytic  expressions 
which  compare  well  with  experiment.  The  common  thread  in  all 
these  approximations  is  that  the  phase  angle  <(>  = arctan  (q^/q2) 
is  temperature  independent.  In  this  section  we  show  why  such  an 
approximation  is  plausible. 


As  a first  step  it  is  convenient  to  express  the  thermodynamic 
potential  in  polar  coordinates  (q^  = ncos<!>,  q2  = nsin<f>). 


<J>  (n,<i>,P,T)  = 4>Q(T)  + n2  + [(3-l  + S2)  + (P2  - 8^  cos4(J) 

4 6 

+ 63sin4({i]  + [ (yl  + y2)  + (Y2  “ Yj.)  cos4<ti  + Y3sin4ij>] 


+ P2  + 

2 


P(a^sin24>  + a2cos2<(>)n  - PE 


The  polarization  can  be  reduced  out  of  this  thermodynamic 
potential  using  the  condition,  3$/3p  = 0 to  give  the  polari- 
zation in  terms  of  the  soft  mode, 

P = x{E  - n2(a1sin2<()  + a2cos2i}>)  } 

and  yield  the  reduced  thermodynamic  potential  involving  only  the 
sof  t mode , 

2 

$(h,<J>»T)  + <t>o  + ~ [ot(T)  + 2xE(a1sin<j>  + a^osiji)] 

4 

+ — (f  + gcos4(}>  + hsin4<()) 

6 

+ L2  C(Y1  + Y2*  + (Y2  “ cos4<J’  + Y3sin44>] 
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where 


(1/2)  (Bj.  + B2)  - (a2  + a2)X 
(1/2) (B2  - B1)  + (a2  + a2)X 


(1/2) - 2a1a2x 


To  determine  the  behavior  of  the  soft  mode  amplitude  one  need 
only  solve  the  equations  34>/94>  = 0 and  3$/3n  = 0 for  n(T) 
and  <j>  (T ) . These  equations  take  the  form: 


2 i 4 

n < n 


[ <Y2  - Y1)sin4<ji  - Y3cos4<f>]/3 


+ n [gsin4<(>  - hcos4<j>]  + 


+ 2xE  [a^cos2<j>  + a2sin2i}i] 


( 4 

n < n [(y1 


+ y2)  + (v2  ~ X2)cos4(t>  + Y3sin4(j)]/2 


+ n [f  + gcos4<(>  + hsin4ij)] 


+ a(T)  + 2xE  [a1sin2<|)  + a2cos2({i]|  = 0 

For  the  solution  to  these  equations  to  be  a minimum  of  $ , one 

2 2 2 

must  also  have  that  iL£  > o and  > 0 since  = 0 at  all 

3*2  3n2  3*9n 

extrema. 


While  detailed  solutions  to  these  equations  require  numerical 
computation,  some  general  characteristics  of  the  solutions  can 
be  discerned  immediately.  The  first  equation  does  not  involve 
the  temperature  dependent  coefficient  a(T);  the  solution  to 
this  equation  then  gives  the  locus  of  points  in  the  n~<J>  plane 
which  are  minima  of  $.  The  temperature  dependence  of  the 
actual  minimum  is  determined  by  the  solution  to  the  second 
equation  which  involves  a(T)  explicitly.  Thus  the  minimum  of 
$ moves  along  the  curve  determined  by  the  first  equation 
at  a rate  determined  by  the  second  equation.  To  illustrate 
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the  possible  behaviors  we  have  plotted  in  Figure  la  and  lb 
the  solution  to  the  first  temperature-independent  equation  for  a 
variety  of  values  of  the  parameters  = Y2  and  Y3  = 0 
in  all  plots  shown).  These  curves  illustrate  the  point  that, 
except  for  small  values  of  n,  there  is  essentially  no  variation 
in  the  angle  <P . It  is  thus  plausible  to  neglect  any  explicit 
angular  dependence  at  the  very  onset. 


2.2.4 


Expressions  to  Compare  with  Experiments 


If  the  angular  dependence  is  completely  eliminated  from  the 
reduced  themodynamic  potential  one  obtains: 

$(n,P,T)  = * (T)  + 2i2l  n2  + 7 n4  + \ n6 
0 2 4 6 


+ aPn2  + P 2 


PE 


2 2 2 

where,  n = q.  + q„  is  the  square  of  the  order  parameter 

1 z -1 
and  the  coefficients  ct  = a0(T-T0),  0 , y,  a,  and  x are  re~ 

lated  to  the  coefficients  appearing  in  the  complete  thermodynamic 

potential.  The  precise  relationship  of  the  coefficients  in  this 

approximate  thermodynamic  potential  to  the  coefficients  in  the 

complete  thermodynamic  potential,  depends,  of  course,  on  which 

approximation  is  adopted.  We  fit  all  of  our  experimental  data 

to  this  form  for  the  thermodynamic  potential. 

The  first  step  in  the  analysis  is  to  reduce  out  the  polarization 
using  the  condition  3<I>/8P  = 0 to  express  the  polarization  in 
terms  of  the  soft  mode  amplitude, 

P ( T , E ) = X [E  - ah2 (T,E ) ] 


and  to  obtain 
mode  alone: 


the  reduced  thermodynamic  potential  for 


o (n)  = 


o (T) 

O 


a (T) 


. 6 4 . Y 

+ 4 n + 6 


the  soft 


where  a(T)  = aQ(T  - TQ)  + 2axE  and  8=6-  2xa2. 

The  temperature  dependence  of  h is  determined  by  the  condition 
34>/3n  = 0,  namely. 


n[a(T)  + 8n2  + xn4]  = 0. 
We  thus  obtain  explicitly  the  solution. 
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n (t,e)  = 


(f)V2  [Ti(E) 


- T] 


1/2 


T < T (E) 
c 

T > T (E) 
c 


where 


TX(E) 


T 

o 


_^E  + j*i_ 

a 4a  Y 

o o 


and 


VE) 


T 

o 


- ^E  + 
a 

o 


16  a y 
o 


We  have  assumed  in  this  solution  that  8 < 0 so  that  a first 
order  transition  occurs  at  Tc,  which  is  always  less  than 


The  polarization  is  then, 


P(T,E)  = 


P2  + K(T  - T) 


1/2 


\ XE 


where 


T < T 

c 


T > T 

c 


Pj_  (E) 
K 


(XE 


_ Ml) 

2Y 


- ax 


1/2 


Evidently  the  effect  of  an  external  electric  field  is  merely  to 
shift  the  transition  temperature;  it  does  not  cause  any  smoothing 
of  the  discontinuities  as  occurs  in  a proper  ferroelectric. 
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The  pyroelectric  coefficient  is  obtained  by  simply  differentiat- 
ing with  the  respect  to  T,  to  give 


P(T,E)  = 


I (T1  ' T>~1/2 


T < T 


T > T 


A singularity  in  the  measurements  at  the  transition  temperature 
is  caused  by  the  discontinuous  jump  in  the  polarization. 

Similarly  the  dielectric  constant  is  obtained  by  differentiating 
the  polarization  with  respect  to  E , to  give 


e(T,E)  = 


1 + x + ^(^  - T) 


1 + X 


T < T 


where 


„ 2 . 1/21/2 

Ki  = a x/a0  y 


Since  Tc  is  always  less  than  T^  there  is  no  diveregence  in 
the  dielectric  constant,  only  a step  discontinuity  at  the 
transition. 

To  find  the  specific  heat,  one  must  first  derive  the  entropy 
from  the  thermodynamic  potential, 

9$  2 n2 

S(n,P,T)  =--[n.P,Tl 


and  then  differentiate  to  obtain  the  specific  heat, 
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3 s 

c(E.T)  = T — (T,E) 


h s 

T [n(T,E) ,P (T,E) , T] 

9T 


c (T, E) 


c (T)  + K-T  [T. 

O 2 1 


co(T) 


T < T 

c 

T > T 

c 


where 

K = a 3/2/47  = K2/4K 

2 O -L 


. > * 


[> 


C0(T)  = T 


32»o 

3T2 


Here  c (T)  is  the  other  contribution  to  the  specific  heat  and 
is  assumed  to  vary  only  slowly  with  temperature. 
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2.3 


Dielectric,  Pyroelectric  and  Thermal  Measurements 


This  section  describes  measurements  of  the  dielectric  constant  e, 
pyroelectric  coefficient  p,  and  specific  heat  c as  functions 
of  temperature  performed  on  gadolinium  molybdate  ( GMO ) , terbium 
molybdate  (TMO ) and  gadolinium  dysprosium  molybdate  (GDMO ) . The 
effect  of  an  applied  electric  field  on  the  dielectric  properties 
is  also  investigated. 

2.3.1  Sample  Preparation 

Thin  slices  of  GMO,  TMO,  and  GDMO  were  cut  from  single  crystals 
normal  to  the  pyroelectric  axis.  Samples  were  polished  to 
thicknesses  from  0.03  to  0.1  mm  and  electroded  on  both  faces 
with  evaporated  antimony.  In  order  to  minimize  edge  effects,  a 
special  technique  was  used  to  obtain  fully  electroded  samples. 
The  technique  is  based  on  the  observations  that  in  GMO  the 
ferroelastic  domain  boundaries  consist  of  (110)  and  (110) 
crystallographic  planes  which  are  also  cleavage  planes  (Ref. 
2-22).  (We  have  confirmed  that  these  observations  apply  to  TMO 
and  GDMO  as  well. ) Samples  were  prepared  using  the  following 
procedure:  (1)  identify  (110)  and  (110)  planes  by  visual 

observation  of  the  domain  pattern  through  crossed  polarizers, 

(2)  evaporate  electrodes,  and  (3)  remove  unelectroded  regions  by 
cleaving  along  (110)  and  (110)  planes.  After  preparation, 
samples  were  mounted  freely  on  3 mil  diameter  silver  leads. 

Before  each  measurement  samples  were  poled  with  a steady  field 
of  ~ 5 kV/cm  applied  during  slow  cooling  through  the  Curie 
point  and  down  to  room  temperature. 
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2.3.2 


Dielectric  Constant 


^ ( 1 ) Experimental  Apparatus 

The  dielectric  constant  was  measured  by  determining  the  ratio  of 
the  capacitance  of  a sample  to  that  of  a reference  capacitor 
placed  in  series  with  it.  A small  ac  signal  (at  most  10  V/cm  at 
f 8 kHz)  was  applied  to  the  sample  and  reference,  and  synchronous 

detection  (i.e.,  a lock-in  amplifier)  was  used  to  measure  the 
voltage  appearing  across  the  reference  capacitor.  Measurements 
were  made  in  the  oven  of  a Deltatherm  III  differential  thermal 
♦ analyzer  (DTA)  which  had  been  modified  to  provide  additional 

feed-throughs  and  connections  for  the  sample  and  thermocouple. 

To  achieve  good  temperature  uniformity,  the  sample  was  enclosed 
in  a closed  copper  cylinder  filled  with  inert  gas  (helium).  A 
chromel-alumel  thermocouple  was  used,  and  the  thermocouple  tip 
placed  with  1 mm  of  the  sample  surface  to  ensure  accurate 
temperature  measurements.  By  connecting  the  outputs  of  the 
lock-in  amplifier  and  thermocouple  to  an  x-y  recorder,  we 
obtained  continuous  recordings  of  the  dielectric  constant  as  a 
function  of  temperature. 

The  capacitance  of  the  reference  was  measured  with  an  accuracy 
of  2%,  and  the  sample  dimensions  determined  to  within  5%,  so 
that  the  dielectric  constant  was  obtained  with  an  accuracy  of 
better  than  10%. 

A second,  completely  independent,  exerimental  apparatus  was  used 
to  cross-check  the  measurements.  This  second  aparatus  employed 
a Boonton  capacitance  meter  at  1 MHz  and  a Hewlett  Packard  4270A 
capacitance  bridge  at  1 kHz,  10  kHz,  100  kHz  and  1 MHz.  The 
test  signals  for  these  instruments  were  less  than  100  mV.  In 
all  cases  the  experimental  data  for  both  apparatuses  agreed 
within  better  than  10%. 
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(2) 


Results  of  Measure mentis 

The  dielectric  constants  measured  on  samples  prepared  as  describ- 
ed above  are  shown  in  Figure  2.  The  results  tor  GMO  and  TMO 
are  in  general  agreement  with  those  reported  by  other  investi- 
gators, although  we  obtain  values  consistently  somewhat  smaller 
than  those  reported  elsewhere.  The  behavior  of  GDMO  is  very 
similar  to  that  of  the  other  two  materials.  No  effect  of 
applied  electric  field  on  the  dielectric  constant  could  be 
observed  in  any  of  these  materials  for  fields  up  to  10  kV/cm. 

2.3.3  Pyroelectric  Coefficient 

( 1 ) Byer-Roundy  Apparatus 

The  pyroelectric  coefficient  was  measured  by  the  direct  tech- 
nique of  Byer  and  Roundy  (Ref.  2-32),  in  which  the  crystal  was 
connected  to  a microvolt  ammeter  (Keithley,  Model  150B)  and 
heated  slowly  at  a constant  rate,  typically  2°C/min.  The 
pyroelectric  current  produced  by  the  sample  was  recorded  direct- 
ly as  a function  of  temperature  using  an  x-y  recorder  and  the 
sample  temperature  was  recorded  separately  as  a function  of 
time.  The  same  samples  and  oven  were  used  as  for  the  dielectric 
constant  measurements. 

The  major  source  of  error  in  these  measurements  is  nonlinearity 
of  the  sample  temperature  as  a function  of  time.  The  measure- 
ments were  repeated  several  times,  and  we  were  able  in  general 
to  achieve  a heating  rate  which  was  constant  to  better  than  5%. 

On  the  few  occasions  when  larger  variations  occurred  they  were 
taken  into  account  by  correcting  the  experimental  data. 

A second,  completely  independent,  Byer-Roundy  apparatus  was  used 
to  cross-check  these  measurements.  This  second  apparatus 
employed  a Keithley  616  Digital  Electrometer.  In  all  cases  the 
results  of  both  apparatus  agreed  well  within  the  experimental 
errors . 
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(2) 


Results  of  Byer-Roundy  Measurements 


Figure  3 shows  the  measured  pyroelectric  coefficients  of  fully- 
electroded  samples  of  GMO,  TMO,  and  GDMO.  The  behavior  of 
the  pyroelectric  coefficient  is  essentially  the  same  in  each. 

No  influence  of  an  applied  electric  field  on  p was  observed  for 
fields  up  to  10  kV/cm. 

Our  data  yield  larger  pyroelectric  coefficients  for  GMO  than 
reported  by  Ullman  et  al.  (Ref.  2-10).  We  obtain  2.0  x 10-9 
C/cm2K  at  1 43 °C  and  2.8  x 10-9  C/cm2K  at  151°C,  which  are, 
respectively,  50%  and  20%  larger  than  the  results  of  Ref.  2-10. 
For  TMO,  we  find  1.8  x 10  9 C/cm2K  at  140°C,  which  agrees 
well  with  the  results  of  Keve  et  al.  (Ref.  2-14)  and  with  the 
value  inferred  from  the  dynamic  measurements  described  below. 

The  importance  of  using  fully-electroded  samples  when  measuring 
the  pyroelectric  coefficient  is  demonstrated  by  the  fact  that 
partially-electroded  samples  exhibit  multiple  peaks  in  the 
pyroelectric  current  similar  to  those  reported  in  Ref.  2-10. 

This  structure  is  reproducible  from  run  to  run  on  a given 
sample,  but  varies  from  one  sample  to  another.  Thicker  samples 
(100-200  nm)  frequently  exhibit  three  or  more  peaks,  while  very 
thin  samples  (25-50  ym)  sometimes  show  a single  main  peak  and 
one  or  more  smaller,  secondary  ones.  As  noted  above,  these 
multiple  peaks  were  attributed  by  Ullman  et  al.  (Ref.  2-12)  to 
the  presence  of  temperature  gradients  in  the  samples.  We  were 
unable  to  eliminate  this  phenomenon  by  using  very  thin  samples 
freely  suspended  in  a strictly  isothermal  enclosure,  as  suggested 
in  Ref.  2-12.  This  led  us  to  consider  the  possibility  that  the 
multiple  peaks  are  due  to  electrical  or  mechanical  boundary 
conditions  imposed  by  partial  electrodes,  and  we  have  found  that 
the  use  of  fully-electroded  samples  prepared  as  described  in 
Par.  2.3.1  greatly  reduces  or  eliminates  this  effect.  When 
full-surface  electrodes  are  used,  even  samples  several  hundred 
microns  in  thickness  show  only  a single  current  peak  on  both 
heating  and  cooling. 
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PYROELECTRIC  COEFFICIENT  (10  C/cm  K) 


(3) 


Chynoweth  Apparatus 


A second,  independent  measurement  of  the  pyroelectric  coefficient 
was  obtained  by  determining  the  ratio  p/c  using  the  dynamic 
response  technique  introduced  by  Chynoweth  (Ref.  2-33).  The 
crystal  was  illuminated  by  a train  of  rectangular  infrared 
pulses  of  short  duration  (typically  15  ms)  so  that  the  increase 
of  its  temperature  did  not  exceed  0.01  K.  The  resulting  pyro- 
electric voltage  response  was  measured  across  a load  resistor 
7 8 

(typically  10  - 10  ft)  and  FET  source  follower  circuit 

connected  in  parallel  to  the  sample  (Ref.  2-34).  The  pyroelectric 
response  voltages  obtained  were  displayed  on  an  oscilloscope, 
and  the  rms  value  of  the  voltage  response  was  measured  with  a 
lock-in  amplifier  (Keithly  Autolock  840).  The  outputs  of  the 
lock-in  and  thermocouple  were  recorded  simultaneously  on  an  x-y 
recorder . 

This  method  yields  relative,  not  absolute,  values  of  p/c  unless 
the  radiation  absorbed  by  the  sample  is  known.  It  is  not 
practical  to  measure  the  incident  flux  or  sample  absorption 
directly,  so  the  system  was  calibrated  by  taking  the  values 
p = 1.8  x 10-9  C/cm2  K and  c = 2.0  J/cm3K  for  TMO  at  140°C. 

( 4 ) Results  of  Chynoweth  Measurements 

The  ratio  p/c  measured  on  f ully-electroded  samples  of  GMO,  TMO, 
and  GDMO  is  shown  in  Figure  4.  It  is  readily  seen  that  both 
the  absolute  magnitude  and  the  temperature  dependence  of  this 
ratio  are  very  similar  in  the  three  materials.  These  results, 
taken  together  with  the  values  of  c reported  in  Par.  2.3.4, 
yield  good  agreement  with  the  direct  measurement  of  p reported 
above . 

Effects  due  to  the  use  of  partially-electroded  samples  are  also 
observable  in  the  dynamic  response  measurements.  Instead  of  the 
sharp  maximum  at  the  transition  temperature  reported  above  for 
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f ully-electroded  samples  and  predicted  by  the  phenomenological 
theory  (Section  4),  we  find  that  the  response  has  a smaller  peak 
than  in  the  fully  electroded  case  and  exhibits  a shoulder  on  the 
high-temperature  side.  The  origin  of  this  behavior  and  its 
relationship  to  the  multiple  peaks  observed  in  the  pyroelectric 
current  are  not  yet  understood.  We  suspect  that  the  solution  to 
this  problem  lies  in  the  coupled  f erroelectric-f erroelastic 
behavior  of  these  materials  which  causes  the  electrical  and 
elastic  boundary  conditions  to  interact.  An  explanation  of 
these  effects  may  also  have  implications  for  many  other  materials 
which  are  simultaneously  ferroelectric  and  ferroelastic. 

2.3.4  Specific  Heat  Measurements 

( 1 ) Description  of  Apparatus 

Specific  heat  measurements  were  made  using  a dc  heating  technique. 
The  instrument  used  was  a Perkin-Elmer  differential  scanning 
colorimeter,  Model  DSC-2,  interfaced  to  a Tektronix  programmable 
calculator.  Model  31.  The  system  is  equipped  with  an  Autoscan 
Zero  mini-computer  for  automatic  baseline  corrections.  Single 

crystal  samples,  typically  between  20  and  80  mg,  were  weighed 

-4 

with  a precision  of  1 (10  ) on  a Perkm-Elmer  autobalance, 

Model  AD-2. 

Open  aluminum  sample  pans  were  used  as  crucibles;  no  external 
mechanical  stresses  were  applied  to  the  samples.  Samples  were 
flat  plates  about  1 mm  thick  and  3 mm  on  a side.  In  all  cases, 
the  flat  side  of  the  sample  made  good  thermal  contact  with  the 
crucible . 

Typical  data  was  taken  using  heating  rates  of  1.25  K/min  in  the 
vicinity  of  the  phase  transition.  Lowering  the  heating  rate  by 
half  did  not  effect  the  observed  data.  In  the  vicinity  of  the 
phase  transition,  data  were  recorded  at  intervals  of  0.1  K.  All 
reported  data  were  taken  during  the  heating  portion  of  a cycle; 
the  observations  made  during  the  cooling  portion  were  in  all 
cases  super imposable  on  the  heating  curves  to  within  1%. 


( 2 ) Accuracy  of  Measurements 

The  calorimeter  was  calibrated  relative  to  the  melting  point  and 
the  enthalpy  of  fusion  of  5N  indium  metal;  Tu  = 429.87  K and 
AHj  = 3266.8  J/mol,  and  of  6N  lead  metal;  TM  = 600.65  K 
and  AHf  = 4768.1  J/mol.  Both  the  melting  points  and  the 
enthalpies  of  fusion  were  found  to  be  independent  of  heating  and 
cooling  rates  from  0.625  to  20  K/min  and  reproducible  to  0.02% 
and  0.75%,  respectively. 

The  raw  date  were  further  refined  by  normalization  to  the  NBS 
values  for  sapphire  (Ref.  2-35),  corrected  on  a sapphire 
standard,  run  under  the  identical  operating  parameters. 

The  overall  accuracy  of  the  absolute  value  of  the  specific  heats 
was  better  than  1%.  Relative  measurements  on  a single  material 
are  accurate  to  0.25%. 

( 3 ) Results 

The  isobaric  molar  heat  capacities  in  the  vicinity  of  the 
ferroelectric  phase  transition  of  the  three  rare-earth  molyb- 
dates are  shown  in  Figures  5 to  7.  The  corresponding  thermo- 
dynamic values  are  given  in  Table  1.  The  transition  temperature 
is  taken  as  the  peak  in  the  heat  capacity.  The  heat  of  transi- 
tion AH  is  obtained  by  graphical  integration  of  the  excess  heat 
capacity . 

The  integration  is  performed  by  extrapolating  a line  from  the 
heat  capacity  curve  in  the  low  temperature  phase  to  a line  in 
the  high  temperature  phase.  This  is  done  over  a temperature 
interval  as  small  as  8 K in  TMO  and  as  large  as  15  K in  GDMO. 

The  entropy  of  the  transition  is  calculated  by  dividing  AH  by 
Tc>  This  assumes  that  the  transition  occurs  over  a small  enough 
temperature  interval  so  that  the  error  in  this  approximation  is 
minimal.  Thus  values  of  AH  and  AS  are  upper  bounds  on  the 
true  latent  heat  and  entropy  discontinuity  at  Tc  as  they 
contain  contributions  also  from  the  rise  in  the  specific  heat 
just  below  the  transition  temperature. 
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TABLE  1.  THERMODYNAMIC  VALUES  FOR  THE  FERROELECTRIC  PHASE 
TRANSITION  IN  THE  MOLYBDATES. 


T 

c 

(K) 

H ( J/mol  1) 

S (J/mol  1K~1) 

MW 

TMO 

431.7 

+ 0.2 

255  + 

26 

0.59 

+ 0.06 

797.66 

GMO 

430.8 

+ 0.2 

304  + 

30 

0.71 

+ 0.07 

794.31 

GDMO 

429.8 

+ 0.2 

288  + 

29 

0.67 

+ 0.07 

799.56 

Two  conclusions 

may  be 

drawn  from  this 

data. 

First,  the 

ferro- 

electric  phase  transition  in  the  molybdates  is  strongly  first 
order  as  evidenced  by  the  large  latent  heat  of  transition  and 
the  sharp  discontinuity  in  the  heat  capacity.  Second,  the 
position  of  the  transition  temperatures,  the  very  similar  heat 
capacities  above  and  below  the  transition,  and  the  similar  heats 
of  transition  suggest  that  the  effect  of  varying  the  rare-earth 
metal  upon  the  crystal  thermodynamics  is  quite  small. 

( 4 ) Multiple  Peaks  in  the  Specific  Heat  of  Rare-Earth 

Molybdates 

In  the  measurements  on  the  specific  heats  of  rare-earth  molyb- 
dates a number  of  samples  exhibited  multiple  peaks  in  the 
specific  heat  near  the  phase  transition  (Fig.  8).  Moreover, 
this  phenomenon  of  multiple  peaking  can  be  induced  in  a sample  by 
subjecting  it  to  mechanical  stresses  (Fig.  9).  The  occur- 
rence of  multiple  peaks  is  then  interpreted  to  be  caused  by 
residual  mechanical  stresses  in  the  samples.  This  interpreta- 
tion is  confirmed  by  annealing  a sample  which  exhibited  multiple 
peaks  and  noting  the  disappearance  of  the  multiple  peaks  in  the 
specific  heat. 

Three  single  crystal  plates  of  TMO  which  exhibited  multiple 
peaking  of  the  heat  capacity  were  wrapped  in  platinum  foil  and 
heated  to  960°C  in  alumina  crucibles  under  a flowing  molecular 
oxygen  atmosphere  of  1.0  1/min.  The  temperature  was  maintained 
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Figure  8.  Specific  heat  of  gadolinium  dysprosium  molyb- 
date vs.  temperature.  This  sample  exhibits  the  phenomenon 
of  multiple  peaks  at  the  transition. 
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Figure  9.  Specific  heat  of  terbium  molybdate  vs.  temperature  follow- 
ing various  mechanical  processings:  (a)  initial  scan,  (b)  after 
applying  finger  pressure  for  5 sec. , (c)  after  applying  finger  pressure 
for  5 sec  more,  (d)  after  pressing  between  two  glass  slides  for  10  sec. 


at  960°C  + 10°C  for  24  hours  and  cooled  down  to  room  temperature 
at  a linear  rate  of  12°C/hr.  Samples  which  were  heated  1080  + 
20°C  melted,  even  though  the  reported  melting  point  for  TMO  is 
1170°C  (Ref.  2-36). 

The  measured  molar  heat  capacity  of  one  sample  of  TMO,  before  and 
after  annealing,  is  shown  in  Figures  10  and  11,  respectively. 

It  is  cleav'  that  proper  temperature  annealing  does  eliminate  the 
multiple  peaking,  although  one  sample  of  TMO  still  showed  some 
evidence  of  residual  multiple  peaking.  The  magnitude  of  cp 
for  the  annealed  samples  is  in  excellent  agreement  with  our 
previous  measurements  on  unannealed  TMO  samples  which  exhibited 
no  multiple  peaking.  It  is  also  worth  noting  that  the  measured 
heat  of  transition  is  the  same  in  the  multiple  and  single  peaked 
samples,  as  one  would  expect  since  the  multiple  peaking  is  a 
strain  induced  effect. 
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Lgure  10.  Specific  heat  vs.  temperature  for  TMO  sample 
.u  i i j.  i 1 ^ noakc  n a r transition  temperature. 
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2.4 


Measurements 

This  section  describes  the  fit  of  the  dielectric  constant, 
pyroelectric  coefficient,  and  specific  heat  measurements  to  the 
expressions  derived  from  the  reduced  thermodynamic  potential 
given  in  Paragraph  2.2.  Excellent  agreement  is  obtained  in  all 
cases.  The  Clausius-Clapeyron  equation,  giving  the  electric 
field  dependence  of  the  transition  temperature,  is  also  checked 
and  found  to  be  in  agreement  with  the  data. 

2.4.1  Analysis  of  Dielectric  Constant 

For  the  temperature  dependence  of  the  dielectric  constant  along 
the  polar  axis,  one  expects  (see  Par.  2.2.4),  in  the  paraelectric 
phase : 

e = ep  = 1 + X 

assumed  to  be  independent  of  T near  Tc;  and  in  the  ferroelectric 
phase: 

e = £p  + K1  (T1  - T ) -1/2  (1) 

where  is  a constant,  and  T^  is  related  to  the  Curie-Weiss 
temperature  Tq  and  the  transition  point  Tc  by: 

T - T - (T  - T )/3.  (2) 

1 c c c 

Figure  12  shows  our  previously  reported  results  for  e in  GMO, 

TMO , and  GDMO , together  with  the  behavior  calculated  from  Equation 
(1)  using  e , K^,  and  T^  chosen  to  give  the  best  fit.  The 
agreement  is  remarkably  good  below  the  transition  temperature, 
and  the  error  is  less  than  5%  above  it.  The  values  of  e , 

, and  T^  obtained  in  this  way  are  listed  in  Table  2. 


DIELECTRIC  CONSTANT 


TEMPERATURE  CC) 

(a)  TMO 


TEMPERATURE  CC) 

(b)  GMO 


++++*  + + 


TEMPERATURE  CC) 

(c)  GDMO 


Figure  12.  Dielectric  constant 
vs.  temperature  for  (a)  terbium 
molybdate,  (b)  gadolinium 
molybdate,  and  (c)  gadolinium- 
dysprosium  molybdate.  The 
crosses  are  the  experimental 
data  and  the  solid  line  rep- 
resents the  theoretical  ex- 
pected temperature  dependence 
fit  to  only  the  data  points 
which  lie  below  the  transition 
tempe  rature . 
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TABLE  2.  PARAMETERS  e , K, , AND  T OBTAINED 
BY  FITTING  THE  OBSERVES  DIELECTRIC  CONSTANT . 


£P 

( K1/2 ) 

T1(°C) 

« 

TMO 

10.6 

1.25 

161.3 

GMO 

9.5 

1.87 

160.2 

t 

GDMO 

9.8 

1.57 

154.7 

2.4.2  Analysis  of  Pyroelectric  Coefficient 

The  expected  (see  Par.  2.2.4)  temperature  dependence  of  the 
spontaneous  polarization  P in  the  ferroelectric  phase  is, 

S 

Ps  = K (T1  - T)1/2  + Px,  (3) 

where,  K = a constant,  = temperature  above  which  no  ferro- 
electric phase  exists,  and  P^  = hypothetical  polarization  at 
T^.  The  temperature  T^  is  defined  above. 

From  Equation  (3)  the  pyroelectric  coefficient  p = - dP  /dT  is 

s 

given  by: 

p = ( K/2 ) (T1  - T)_1/2  (4) 

In  order  to  find  the  numerical  values  of  K and  T. , we  plotted 

2 -*■ 

1/p  as  a function  of  the  temperature  T,  for  TMO , GMO,  and  GDMO 

(Fig.  13).  It  can  be  seen  that  the  experimental  points  lie  along 

straight  lines  as  expected  from  Equation  (4).  The  agreement  is 

excellent,  and  the  values  obtained  for  K and  T^  are  given  in 

Table  3. 
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TABLE  3.  PARAMETERS  K AND  T,  OBTAINED  BY  FITTING 
DYNAMIC  MEASUREMENTS  OF  THE  PYROELECTRIC  COEFFICIENT. 


K 

*1 

(10"8C/cm2K1/2 ) 

(°C) 

TMO 

1.77 

161.3 

GMO 

1.72 

159.6 

GDMO 

1.66 

154.9 
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2.4.3  Transition  Temperatures 


Defining  Ae  = e - £p,  we  can  derive  from  Equations  (1)  and  (4) 
a useful  relationship  for  determining  Tq: 

Ae(Tc)/Ae(TQ)  = 2 (5) 


Using  Equation  (5)  and  the  experimental  data  shown  in  Figure 
we  have  determined  Tq  for  each  of  the  materials  studied.  The 
results  are  given  in  Table  4r  together  with  the  observed  transit- 
ion temperatures  Tc  and  the  values  obtained  for  T^  above. 

Slightly  different  values  of  T^  are  obtained  by  fitting  the 
pyroelectric  coefficient  and  the  dielectric  constant,  but  the 
differences  are  not  significant  in  view  of  the  accuracy  with 
which  we  can  determine  the  sample  temperature  and  the  uncertainty 
in  Tc  caused  by  the  experimental  width  of  the  transition. 


TABLE  4.  THE  CURIE-WEISS  TEMPERATURE  T AND  TRANS- 
ITION TEMPERATURE  T AS  DETERMINED  FROM°EQ  . (15)  AND 
THE  MAXIMUM  OF  £,  RESPECTIVELY.  ALSO  GIVEN  FOR  COM- 
PARISON IS  T1  AS  FOUND  FROM  FITTING  THE  PYROELECTRIC 
COEFFICIENT  p or  THE  DIELECTRIC  CONSTANT  e.  ALL 
TEMPERATURES  ARE  IN  °C . 


To 

T 

C 

T^  (from  p) 

T^  (from  e ) 

TMO 

159.7 

161.0 

161.3 

161.3 

GMO 

157.0 

159.5 

159.6 

160.2 

GDMO 

152.6 

154.2 

154.9 

154.7 

t 


2.4.4  Influence  of  an  Applied  Electric  Field 

The  expected  value  of  the  rate  of  change  of  the  transition 
temperature  with  external  field  can  be  given  in  terms  of  the 
constants  K and  K^,  namely, 

dTc/dE  = 2K.,/K  . (6) 

Inserting  the  experimental  values  for  GMO  leads  to  dT  /dE  = 

-5-1  c 

1.9  x 10  K/V  cm  . This  value  is  quite  small  compared 

to  the  result  for  the  proper  ferroelectric  BaTiO,,  where 

-3-1  J 

dTc/dE  = 1.6  x 10  K/V  cm  , and  explains  why  we  have 

observed  no  change  in  Tc  in  the  p or  e measurements  at 

fields  up  to  10  kV/cm. 

A general  thermodynamic  argument  based  on  the  Clausius-Clapeyron 
equation  shows  that, 

dTc/dE  = AP/AS,  (7) 

where,  AP  and  AS  are  the  discontinuous  changes  in  the  polar- 
ization and  entropy,  respectively,  which  occur  at  the  phase 
transition.  Inserting  the  experimental  values  (Refs.  2-13, 

2-19,  2-30)  for  GMO  (Ap  = 4 x 10-8  C/cm2  and  As  = 2.2  x 
10-3  J/cm3K ) , yields  dTc/dE  = 1.8  x 10-5  K/V  cm-1,  in 
very  good  agreement  with  the  result  obtained  above  from 
Equation  (6). 
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2.4.5  Analysis  of  Thermal  Measurements 

The  specific  heat  measurements  can  be  compared  to  the  theoretical 
expressions  for: 

(a)  Discontinuity  of  the  entropy  at  the  transition  temperature, 


AS  = (3K2/2K1)  (T1-Tc)1/2 


(b)  Latent  heat  at  the  Curie  point. 


AQ  = TcAS  = (3K2Tc/2K1)  (T^T^172 


(c)  Temperature  dependence  of  the  specific  heat, 


c = 


cQ  + (K2/4K1)  T (T1-T)"1/2  for  T < Tc 

c for  T > T 

o c 


The  values  of  K,  , T^  and  Tc  are  determined  above  from  the 
analysis  of  the  pyroelectric  and  dielectric  properties.  Table  5 
summarizes  these  values,  using  the  T^  value  determined  in  the 
dielectric  measurements. 


A comparison  of  the  predicted  values  and  measured  upper  bounds 
for  AS  and  AQ  is  given  in  Table  6. 


TABLE  5.  PARAMETERS  DETERMINED  PREVIOUSLY  FROM 
PYROELECTRIC  AND  DIELECTRIC  MEASUREMENTS. 


Tc 

T1 

K 

K1 

(K) 

(K) 

(10_8C/cm2K1/2) 

(K172 

TMO 

434.2 

434.5 

1.77 

1.25 

GMO 

432.7 

433.4 

1.72 

1.87 

GDMO 

427.4 

427.9 

1.66 

1.57 
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TABLE  6.  COMPARISON  OF  CALCULATED  VALUES  AND  MEASURED  UPPER 


BOUNDS  FOR  ENTROPY 

DISCONTINUITY 

AND  THE  LATENT 

HEAT  . 

AS ( J/mol  K) 

AQ( J/mol) 

Calculated 

Measured 

Calculated 

Measured 

Upper  Bound 

Upper  Bound 

TMO 

0.40 

0.59  + 0.06 

173 

255  + 26 

GMO 

0.39 

0.71  + 0.07 

169 

304  + 30 

GDMO 

0.36 

0.67  + 0.07 

157 

288  + 29 

Figure  14  shows  experimental  results  for  c in  TMO  (Isomet  #1)  as 
measured  (the  temperature  scale  of  the  experimental  data  was 
shifted  by  2.5  K to  agree  with  T£  = 161°C  as  determined  from 
our  dielectric  and  pyroelectric  measurements)  together  with 
the  behavior  calculated,  using  the  data  of  Table  5 and  choosing 
c to  give  the  best  fit.  For  c we  found  a constant  value 
of  2.07  J/cnrK  over  the  temperature  range  150-165°C.  The 
agreement  between  the  calculated  curve  and  the  experimental 
points  is  good  except  in  the  immediate  vicinity  of  the  transi- 
tion point.  In  this  region  an  agreement  between  the  theory  and 
the  experiment  is  not  expected  because  the  anomaly  in  the 
specific  heat  at  the  transition  is  always  spread  over  a small 
temperature  range  due  to  instrumental  resolution. 


SPECIFIC  HEAT  (J/crn  K) 


TEMPERATURE  CK) 

-j 

Figure  14.  Specific  heat  of  TMO  (Isomet  #1). 

Crosses  are  experimental  values  previously 
reported;  solid  line  is  theoretically  cal- 
culated behavior- 
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2 . 5 Figures  of  Merit 

The  figure  of  merit  of  a pyroelectric  material  used  as  a vidicon 

1 1/2 

target  is  M(l)  = p/ec,  while  M^)  = p/e  ' c is  the  figure  of 
merit  for  a small-area  detector  made  of  that  material.  We  have 
developed  an  experimental  technique  which  permits  direct  measure- 
ment of  these  figures  of  merit,  reducing  or  completely  eliminat- 
ing the  need  to  calculate  products  or  ratios  of  independent 
measurements.  Results  obtained  by  this  method  are  reported 
below  for  the  rare-earth  molybdates.  To  permit  comparison  with 
the  behavior  of  a typical  proper  ferroelectric,  we  also  present 
data  for  TGS. 

2.5.1  Measurement  Method 

The  measurement  method  is  a modification  (Ref.  2-34)  of  the 
dynamic  technique  of  Chynoweth  (Ref.  2-31),  in  which  we  measure 
voltage  rather  than  current  response.  A brief  analysis  of  the 
technique  is  given  below. 

Consider  a thin  plate  of  pyroelectric  material  electroded  on  its 
major  faces  which  are  perpendicular  to  the  polar  axis.  In 
parallel  to  the  sample  are  connected  a load  resistor  R and  a 

Li 

preamplifier  of  input  capacitance  and  input  resistance  R^. 

A chopped  beam  of  radiation  with  a 50%  duty  factor  (Fig.  15a) 
is  incident  upon  the  front  electrode  and  heats  the  sample.  The 
chopping  period  t is  assumed  to  be  much  shorter  than  the 
thermal  time  constant  of  the  sample.  It  can  be  shown  that  in 
the  steady  state,  the  mean  temperature  T of  the  sample  rises 
linearly  during  a radiation  pulse  and  falls  linearly  between 
pulses  (Fig.  15b).  The  slopes  of  the  rise  and  fall  are  equal 
(except  for  sign)  and  given  by  + Fq/2cL,  where  Fq  = radiation 
power  absorbed  per  unit  area  of  the  sample  during  a single 
pulse,  and  L = thickness  of  the  sample.  As  a result  of  the 
temperature  variations,  the  spontaneous  polarization  of  the 
material  is  changed,  causing  a charge  flow  through  the  sample 
and  the  load. 


» 
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The  pyroelectric  current  (Ref.  2-37)  pA(dT/dt)  consists  of 
positive  and  negative  square  pulses  (Fig.  15c)  of  height 
pAFQ/2cL,  where  A = electrode  area.  For  the  voltage  V(t) 
developed  across  the  sample  and  the  load,  one  obtains  the 
equation : 


d(CV) 

dt 


dT 

dt  ' 


(1) 


where  C represents  the  sum  of  and  the  sample  capacitance 
C , and  R the  parallel  combination  of  RT , R.  , and  the  sample 
resistance  R . Normally  R.  and  R are  several  orders  of 

S 1 o 

magnitude  larger  than  Rr , so  that  to  a good  approximation  R = R 
Solving  Eq.  (1)  one  finds  that  various  response  shapes  can  be  ob 
tained  depending  on  the  ratio  between  the  electronic  time  con- 
stant of  the  system  Tg  = RC,  and  the  period  t of  the  radiation 
pulses.  When  r is  much  greater  than  Tg,  a rectangular  re- 
sponse is  obtained  (Fig.  15d).  In  this  case  the  peak-to-peak 
value  V of  the  response  is  independent  of  t and  directly 
proportional  to  the  ratio  (p/c): 

VR  = (FqARl/L) (p/c).  (2) 

When  t is  much  less  than  t , a triangular  response  is 
obtained  (Fig.  15e)  for  which  the  peak-to-peak  value  VT  is 
proportional  to  the  chopping  period  t, 

VT  = (FqAt/4L) (p/cC)  . (3) 

When  the  capacitance  Cg  of  the  sample  dominates,  C - Cc  = eQeA/L 
and  VT  is  directly  proportional  to  M^.  If  is  not  negligible, 
can  be  expressed  in  terms  of  VT  and  a correction  factor  g, 
which  can  be  calculated  from  the  ratio  (VT/VR) : 

M1  = (4g/For)VT  , (4) 

where 

g = [1  " 4RLCiT_1(VT/VR)]'1. 
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Using  Eqs.  (2)  and  (4),  M1//2  can  exPressed  in  the  form: 

Ml/2  = (2/F0)(gL/ARLx)1/2  (V^)172  (5) 

Equations  (4)  and  (5)  form  the  basis  for  our  direct  experimental 
determinations  of  the  figures  of  merit. 

2.5.2  Experimental  Techniques 

Thin  slices  of  TGS,  GMO,  TMO,  and  GDMO  were  cut  from  single 
crystals  normal  to  the  pyroelectric  axis.  Samples  were  polished 
to  thicknesses  from  0.C3  to  0.1  mm  and  electroded  on  both  faces 
with  evaporated  antimony.  Attached  to  the  electrodes  were  3 mil 
diameter  silver  leads  by  which  the  samples  were  suspended.  A 
small  sample  holder  with  a LiF  window  was  used,  and  the  sample 
temperature  was  measured  to  ~ 0.2  C with  a calibrated  thermo- 
couple attached  very  close  to  the  rear  electrode.  The  sample 
mount  was  attached  to  a heated  copper  block  and  placed  in  a 
vacuum-tight  chamber  with  a LiF  window.  Before  each  measure- 
ment, samples  were  poled  with  a steady  field  of  2-5  kV/cm 
applied  during  slow  cooling  through  the  Curie  point  and  down  to 
room  temperature.  The  radiation  source  was  a Globar  with  a 
variable  speed  chopper  (Ithaco  382A).  To  provide  rectangular 
pulses,  the  width  of  the  chopper  teeth,  as  well  as  the  distance 
between  them,  was  much  larger  than  the  diameter  of  the  radiation 
beam. 

5 9 

The  measuring  circuit  consisted  of  a load  resistor  (10  - 10  ft ) 

and  a FET  source  follower  circuit  connected  in  parallel  to  the 
sample  (Ref.  2-38).  The  pyroelectric  response  voltages  obtained 
were  displaced  on  an  oscilloscope,  and  the  rms  value  of  the 
voltage  response  was  measured  with  a lock-in  amplifier  (Keithley 
Autolock  840).  The  outputs  of  the  lock-in  and  thermocouple  were 
recorded  simultaneously  on  an  x-y  recorder. 


I 
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2.5.3  Experimental  Results 

Figure  16  shows  the  temperature  dependence  of  p/c  and  p/eeoc 
in  TGS  as  measured  by  the  dynamic  method.  The  calibration  in 
this  measurement  was  obtained  by  assuming  for  the  pyroelectric 
coefficient  (Ref.  2-39)  p = 3 x 10  ® C/cm2  K and  for  the  specific 
heat  (Ref.  2-40)  c = 2.2  J/cm^  K at  25°C.  Comparing  the  be- 
havior of  p/c  and  p/eeQc  shows  that  e plays  a dominant  role 
in  determining  the  temperature  dependence  of  ; the  increase 
in  the  dielectric  constant  causes  a steady  drop  in  as  the 
transition  point  is  approached. 

Figure  17  shows  the  temperature  dependence  of  p/eeQc  in  GMO, 

TMO , and  GDMO.  The  results  for  p/c  in  these  materials  were 
presented  in  Par.  3.3.4.  In  the  molybdates,  in  sharp  contrast 
to  the  behavior  exhibited  by  TGS,  p/eecc  increases  in  the 
same  manner  as  p/c  as  the  transition  point  is  approached.  This 
result  is  expected,  since  at  zero  stress  these  materials  exhibit 
only  a small  increase  in  dielectric  constant  near  the  transition 
temperature  (Par.  3.2.2). 


2.5.4  Analysis  of  Results 

The  results  of  Section  2.2  lead  to  predictions  for  both  figures 
of  merit  in  the  molybdates: 

M ( 1 ) = K/{2(Ep  + K^Tj-T)"172]  [cq(T1-T)1/2  + K2T/4K1]}, 

M (i)  = K/{ 2 [ep  + K1(T1-T)_1/2]1/2  [co(T1-T)1/2  + K2T/4K1]}. 

In  Figures  18  and  19  the  calculated  values  of  M(l)  and 
M(y)  are  compared  with  our  experimental  results  for  TMO. 

The  crosses  in  Figure  18  represent  the  experimental  values  of 
p/ec  as  measured  directly  by  the  dynamic  method.  In  Figure 
19  the  crosses  represent  experimental  values  of  p/e*7^c 
obtained  by  taking  the  square  root  of  the  product  of  the  measured 
values  of  p/ec  and  p/c.  The  solid  lines  were  calculated  using 
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for  K,  , T^,  and  Tc  the  values  given  in  Table  5.  For  ep 
we  used  a value  of  10  which  is  ~ 5%  less  than  the  value  of 
10.6  obtained  from  independent  measurement  of  e.  For  cq  we 
chose  the  value  1.97  J/cm^  K which  gives  a better  fit  to  the 
specific  heat  data  of  TMO  #2.  Clearly,  there  is  good  agreement 
between  theory  and  experiment  for  both  figures  of  merit. 

2.5.5  Optimization  of  M(l) 

It  is  readily  seen  from  the  expression  for  M(l)  obtained  in 

previous  section  that  this  figure  of  merit  attains  a maximum  at 

T unless  unrealistically  small  values  of  e and  c are 
c p o 

assumed.  (Even  in  such  a case  the  value  of  at  its  maximum 
differs  only  insignificantly  from  its  value  at  Tc> ) One  can 
evaluate  at  Tc  and  obtain  the  peak  values  of  the 
f igure-of-merit.  This  peak  value  is  a function  of  both  the 
readily  measured  or  estimated  quantities  ep,  co,  and  Tc 
and  the  more  difficult  to  assess  quantities  (T^  - Tc),  K and 
Kj  (which  are  expressible  in  terms  of  the  free  energy  parameters). 
This  peak  value  of  the  figure  of  merit  can  now  be  optimized  with 
respect  to  variations  in  the  free  energy  parameters  to  give  an 
optimum  f igure-of-merit, 

M °Pt  = i (e  c T (6) 

1 4 'epco  c' 

This  optimum  is  attained  for  values  of  the  free  energy  parameters 
such  that, 


Tx  - Tc  = K2Tc/4cQep,  (7a) 

(K/K1 ) = (4cQ/epTc)1/2.  (7b) 

For  the  three  molybdates  examined,  the  free  energy  parameters 
are  such  that  this  optimum  value  is  very  nearly  attained.  Thus, 
this  expression  for  provides  a useful  guide  to  how 

large  a figure  of  merit  might  be  realized  in  the  molybdate 
family . 
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For  materials  in  the  molybdate  family,  cq  is  about  2 j/cm3 
K,  and  Cp  is  about  10;  therefore,  the  principal  variation  in 
stems  from  variations  in  T^.  As  Tc  varies  from  160°C 
down  to  40°C , only  varies  from  9.0  x 103  cm2/C  up  to 

10.6  x 10  cm  /C.  For  compariosn  one  should  recall  that  TGS 
has  a figure  of  merit  of  - 3.0  x 103  cm2/C  at  40°C.  (Fig.  16) 
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3. 


SULFATE  FAMILY 


Five  materials  from  the  family  of  ferroelectric  sulfates  were 

studied:  LiNH4SC>4,  LiND4SC>4,  LiCsSC>4,  LiKS04  and  LiH3(SeC>3)2. 

Lithium  ammonium  sulfate  and  its  deuterated  analogue  exhibited 

similar  properties.  The  high  room-temperature  value  for  the 

figure  of  merit  (about  equal  to  that  of  TGS  at  40°C)  suggested 

further  study  of  related  materials.  Both  lithium  cesium  sulfate 

and  lithium  potassium  sulfate  were  consequently  investigated; 

these  materials,  however,  exhibited  distinctly  different  behaviors 

which  eliminate  them  for  practical  pyroelectric  applications. 

Lithium  trihydrogen  selenite  was  investigated  because  of  its 

2 

high  value  of  spontaneous  polarization  (~  15  yC/cm  ) and  the 
fact  that  the  related  material  RbH3(SeC>3)2  has  been  identified 
as  an  improper  ferroelectric.  We  find  no  evidence  that 
LiH3(Se03)2  is  an  improper  ferroelectric.  The  measured  value  of 
the  ratio  p/e  is  about  one-fourth  that  of  TGS.  While  this 
eliminates  LiH3(SeC>3)2  for  practical  applications,  the 
large  value  of  the  spontaneous  polarization  suggests  that 
related  materials  might  have  potentially  useful  figures  of 
merit. 


3 . 1 Lithium  Ammonium  Sulfate  and  Deuterated  Lithium 

Ammonium  Sulfate 

The  room  temperature  value  of  the  f igure-of-merit  for  these 

3 2 

materials  is  about  3 x 10  cm  /C . On  the  basis  of  this  high 
value  we  conclude  that  further  investigation  of  related  mater- 
ials would  be  worthwhile. 

Pyroelectric  and  second  harmonic  generation  measurements  show 
the  low  temperature  (below  273  K)  phase  of  these  materials  to  be 
polar,  thus  clarifying  an  ambiguity  in  the  literature. 


3.1.1  Growth  and  Characterization  of  Materials 

Large  (3x3x3  cm)  colorless  single  crystals  of  LiNH4S04  were 
grown  from  aqueous  solution  containing  equimolar  concentrations 
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of  (NH^^SO^  and  Lj^SO^  H2O.  The  deuterated  crystals  were 

grown  from  D^O  (99.9%  D)  solutions  saturated  with  previously 

recrystallized  LiND^SO^.  The  final  growth  solution  had  a 

deuteration  level  of  97.7  + 0.03%.  The  evaporation  method  of 

3 

growth  was  employed  using  a nitrogen  gas  flow  of  2 ft  /hr.  The 

3 

rate  of  solvent  removal  was  controlled  at  2 cm  /day.  The 
saturated  LiNH^SO^  and  LiND^SO^  solutions  were  at  pH  3.88 
and  3.93,  respectively.  The  growth  solution  was  controlled  at 
318.00  + 0.05  K and  stirred  at  25  rpm.  Single  crystals  from 
spontaneously  nucleated  seeds  were  always  twinned  and  had  a 
pseudohexagonal  morphology  (Fig.  20).  Untwinned  sections  of 
these  crystals  were  selected  and  only  c-axis  plates  were  used  as 
seeds  for  subsequent  growth. 


Figure  20:  Pseudo-hexagonal  Morphology  of* 
LiNH.SO.  (Crossed  Polarizers). 


There  exists  some  confusion  in  the  literature  on  the  notation  of 
the  crystallographic  axes  for  LiNH^SO^.  We  used  the  conven- 
tion c < a < b for  the  room  temperature  orthorhombic  lattice. 

O O 

Thus  the  lattice  parameters  are  a = 8.78  A,  b = 9.12  A,  and 
c = 5.26  A. 

3.1.2  Specific  Heat  Measurements 

Measurements  of  the  temperature  dependence  of  the  specific  heat 
were  performed  with  a Perkin-Elmer  differential  scanning  calori- 
meter Model  DSC-2.  Sample  weights  were  in  the  10  to  40  mg  range. 
Data  were  recorded  at  a scan  rate  of  10°  min  and  a sensitivity 
of  2 meal  sec  The  apparatus  was  calibrated  with  high 
purity  indium  (Tu  = 429.87  K and  AH  = 6.80  cal  g-3 ) . The 
error  in  H for  indium  was  less  than  0.1%  while  the  error  in  T„ 

M 

was  0.3  K.  The  estimated  absolute  accuracy  of  these  measurements 
is  + 5%. 

The  temperature  dependence  of  the  specific  heat  for  LiNH4SC>4 
and  LiND^SO^  is  shown  in  Figure  21.  Upon  deuteration  of 
LiNH^SO^,  the  transition  temperature  shifts  upward  8.2  K; 
this  is  at  variance  with  the  results  of  our  dielectric  and 
pyroelectric  measurements  where  a shift  of  only  about  2 K is 
observed.  The  experimental  data  for  both  crystals  are  summarized 
in  Table  7.  A least  squares  fit  of  these  data  over  the  range 
415  K to  445  K yields,  for  LiNH4SC>4 

cp  (cal.g-1°K_1 ) = - 0.118  + 1.17  (10_3)T, 

and  for  LiND.SO. 

4 4 

Cp  ( cal . g-1 °K_1 ) = 0.054  + 1.01  (10_3)T. 

Figure  22  illustrates  a high  precision  scan  of  LiNH4s04  (dT/dt  = 
0.6°  min  and  1 mcal/sec  sensitivity)  which  revealed  multiple 
peaking  behavior  near  the  ferroelectric  transition  temperature. 


95 


• I D3) 


° Li N H4SO4  Tc  = 46I-5°K±0.5* 
• Li  NO4SO4  Tc  = 469.7°  K ± 0.5° 
[9  7 % D] 


TEMPERATURE (K) 

Figure  21.  Temperature  behavior  of  c in  LiNH  SO  and 

LiND  SO  . P 

4 4 
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SPECIFIC  HEAT  (ARBITRARY  UNITS) 


^ 1 T 


Similar  behavior  has  been  observed  in  Gd2(Mo04)^,  Tb2(MoC>4)2, 
Cr^B^O^Cl  and  Ni^B^O-^Br*  and  can  be  attributed  to  bulk 
strains  in  the  crystalline  samples  (see  Par.  2.3.4). 

The  specific  heat  measurements  also  show  that  both  LiNH4S04  and 
LiND4SC>4  are  stable  in  air  up  to  at  least  548  K.  This  contra- 
dicts an  earlier  report  that  LiNH4SC>4  decomposes  at  465  K 
(Ref.  3-1). 

3.1.3  Dielectric  Measurements 

Samples  for  electrical  measurements  were  fabricated  and  their 
surfaces  polished  with  a water,  ethanol  and  glycerine  mixture.  A 
domain  structure  in  the  form  of  narrow  stripes  were  clearly 
observed  with  domain  boundaries  parallel  to  the  (100)  planes. 
Electrodes  were  formed  by  deposition  of  antimony. 

( 1 ) High  Temperature  Phase  Transition 

In  contrast  to  previous  reports,  (Refs.  3-2,  3-3)  but  in  agree- 
ment with  Mitsui  et  al.,  (Ref.  3-4)  we  observed  that  the  polarity 
of  LiNH4S04  can  be  reversed  at  room  temperature  by  application 
of  a dc  electric  field  along  the  c-axis.  However,  due  to  the 
large  room  temperature  coercive  field  (~  50  kV/cm)  (Ref.  3-4), 
we  were  not  able  to  generate  good  60  Hz  hysteresis  loops  except 
at  elevated  temperatures.  Figure  23  shows  a typical  loop 
obtained  at  437  K.  The  spontaneous  polarization  and  coercive 
field  estimated  from  this  loop  LiNH4S04  are  0.17  uC/cm  and 
14  kV/cm. 

The  temperature  dependence  of  the  dielectric  constant  of 
LiNH4SC>4  was  measured  under  a helium  atmosphere  at  8 kHz.  At 
293  K,  a value  of  e = 11.2  is  obtained,  and  then  rises  slow- 
ly with  increasing  temperature  to  13  and  then  jumps  to 
about  22  at  Tc  = 459  k (Fig.  24).  The  dielectric  constant 
follows  a Curie-Weiss  law  in  a narrow  temperature  region  above 
Tc  as  previously  reported  by  Mitsui  et  al.  (Ref.  3-4). 

Between  Tc  and  463  K,  the  temperature  dependence  of  e is 
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given  by  e = ec  + C/(T-To),  with  ec  = 13.2,  Tq  = 

and  C = 1.8  K.  This  is  illustrated  in  Figure  25.  We  observed 

no  appreciable  effects  from  dc  bias  fields  up  to  2 kV  cm'^  on 

ec  and  Tc.  The  electrical  properties  of  LiND^SO^  are  almost 

identical  to  the  results  obtained  for  LiNH.SO..  The  ferro- 

4 4 

electric  transition  temperature  is  461  + 1 K.  From  these 
results  we  conclude  that  the  hydrogen  bonding  in  these  crystals 
is  only  weakly  coupled  to  the  ferroelectric  transition. 

( 2 ) Low  Temperature  Phase  Transition 

The  dielectric  constant  of  fully  electroded  samples  of  LiND4SC>4 
was  measured  in  the  vincinity  of  the  low  temperature  transition 
from  80  Hz  to  8 kHz.  Within  1°C  of  the  low  temperature  transition 
the  measurement  becomes  unstable.  A step  change  of  the  dielec- 
tric constant  was  consistently  observed  between  the  upper  and 
lower  phases  (Fig.  26). 

3.1.4  Pyroelectric  Measurements 
( 1 ) High  Temperature  Phase  Transition 

Using  the  dynamic  pyroelectric  response  method,  the  temperature 
dependence  of  the  f igure-of-merit  (M  = p/eeQc)  in  a crystal 
of  LiNH^SO^  over  a range  20°C-200°C  was  measured  (Fig.  27). 

At  room  temperature,  M ~ 3000  cm^/C , which  is  about  that 
obtained  for  TGS  at  40°C.  Upon  heating  the  sample  from  room 
temperature  up  to  - 170°C,  M decreases  by  a factor  of  2.5  and 
then,  as  the  temperature  is  increased,  M increases  rapidly 
passing  through  a sharp  peak  at  the  Curie  point  (Tc  = 185  + 

1°C).  Similar  behavior  is  obtained  in  the  measurement  of  the 
ratio  of  pyroelectric  current  to  specific  heat  (p/c)  as  shown  in 
Figure  28.  Comparison  of  Figures  27  and  28  shows  that  the 
influence  of  the  dielectric  constant  in  determining  the  temper- 
ature dependence  of  M is  only  appreciable  in  the  phase  transition 
region.  From  the  ratio  p/c  and  p/eeQc,  we  obtained  values 
of  10.5  and  17  for  e at  room  temperature  and  at  T , re- 
spectively. 
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The  application  of  a bias  field  (5  kV/cm)  along  the  c-axis  was 
observed  to  enhance  the  peak  value  of  p/c  by  ~ 35%  and  induce  a 
small  pyroelectric  response  above  Tc>  Within  our  maximum 
temperature  resolution  (~  0.2°C),  we  observed  no  shift  in  the 
transition  temperature. 

Byer-Roundy  measurements  on  c-cut  LiNH^SO^  samples  were  also 
performed  in  vacuum  and  helium.  In  both  cases  large  background 
currents  are  observed  near  Tc»  The  pyroelectric  coefficient 
at  50°C  determined  by  this  method  is  about  6 x 10  ^ C cm  ^°C 
Measurements  on  "a"  and  "b"  cut  samples  also  exhibit  the  large 
background  current.  This  spurious  current  is  a bulk  property 
and  could  possibly  be  induced  through  some  defect  (i.e.,  impurity) 
mechanism  or  some  contact  potential. 

(2)  Low  Temperature  Phase  Transition 

The  dynamic  measurement  of  the  temperature  behavior  of  p/c  is 
shown  in  Figure  29  for  LiNH^SO^ . Upon  cooling,  p/c  increases 
slowly  and  becomes  nearly  constant  near  the  transition  temperature 
(Tfc  = 5.5°C),  while  the  pyroelectric  response  drops  sharply 
and  becomes  unstable  so  that  a value  for  p/c  cannot  be  determined. 
Below  4°C,  the  pyroelectric  response  becomes  stable  with  a value 
of  p/c  which  is  about  20%  of  the  value  obtained  above  5.5°C.  On 
heating,  the  pyroelectric  response  is  stable  and  nearly  constant 
up  to  9°C,  then  increases  and  becomes  unstable.  A sharp  increase 
is  observed  at  11°C  after  which  the  response  again  becomes 
stable.  In  the  deuterated  crystals,  the  transition  is  shifted 
toward  lower  temperature  Tfc  = 1.5°C  (cooling)  and  a temperature 
hysteresis  of  about  6°C  is  observed.  The  lower  transition  was 
also  examined  in  both  LiNH^SO^  and  LiND^SO^  by  measuring 
the  temperature  behavior  of  the  pyroelectric  current  during 
heating  and  cooling.  The  pyroelectric  current  of  a partially 
electroded  sample  of  LiND4SC>4  near  Tfc  exhibited  the  same 
multiple  peaking  behavior  seen  in  partially  electroded  samples 
of  Gd^  (MoC>4  ) 2 • A fully  electroded  sample  of  LiND4SC>4  re- 
sult*. - in  a marked  reduction  of  the  number  of  multiple  peaks 
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Figure  29.  Temperature 
dependence  of  p/c,  in 
LiNH^SO  upon  heating 
and  cooling  through  the 
low  temperature  phase 
transition . 


but  did  not  completely  eliminate  this  behavior  (Figs.  30  and 
31).  We  conclude  that  the  reduction  of  clamping  achieved  by  the 
use  of  full  face  electrodes  coupled  with  the  fact  that  thermal 
annealing  eliminates  the  multiple  peaking  behavior,  is  conclusive 
proof  that  this  effect  can  be  attributed  to  strains  present  in 
the  bulk  crystal.  These  strains  are  induced  in  the  crystal 
during  growth  and/or  during  sample  fabrication.  It  is  clear 
that  all  quantitative  measurements  should  be  performed  only  on 
suitably  annealed  samples. 

From  the  sign  of  the  pyroelectric  current  released  near  Tc,  we 

deduce  that  the  spontaneous  polarization  of  LiND4S04  on  the 

low  temperature  side  of  Tfc  is  smaller  than  that  above  Tfc. 

The  jump  in  spontaneous  polarization,  measured  by  integration  of 
, _ o 

the  pyroelectric  current  released  is  0.27  pC/cm  (+  10%). 
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3.1.5  Second  Harmonic  Generation  Measurements 

The  nature  of  the  two  phase  transitions  in  LiNH4S04  and  LiND4SC>4 

was  investigated  using  standard  optical  second  harmonic  tech-  ® 

niques  (SHG)  (Ref.  3-5).  Index  matching  oils  were  used  only  for 

temperatures  in  the  range  293  to  313  K.  The  average  index  for 

LiNH4S04  is  1.496  and  for  LiND4S04  is  1.500.  Both  materials 

were  found  to  be  phase  matchable  for  a 1.06  ym  fundamental  at  ® 

room  temperature.  Table  8 summarizes  the  SHG  analysis.  Contrary 

to  data  previously  reported,  the  low  temperature  phase  in  both 

compounds  is  acentric  (Refs.  3-2,  3-4).  Therefore,  the  symmetry 

of  this  phase  must  be  either  m or  2 and  confirms  two  previous  * 

assignments  (Refs.  3-3,  3-6).  There  is  a pronounced  thermal 

hysteresis  (Fig.  32)  in  the  temperature  behavior  of  the  SHG 

intensity,  on  the  order  of  6°  for  the  high  temperature  transition 

and  10°  at  the  lower  transition.  This  hysteresis  in  the  low  " 

temperature  transition  is  somewhat  larger  than  that  observed  in 

our  pyroelectric  measurements.  For  the  high  transition,  the 

observed  hysteresis  is  at  variance  with  our  dielectric  and 

pyroelectric  measurements  in  which  no  detectable  hysteresis  is  * 

observed. 

3.1.6  Conclusions 

( 1 ) Pyroelectric  Detector  Applications 

The  room  temperature  value  for  the  f igure-of-merit  is  about  70% 
of  that  of  TGS ; moreover  it  varies  only  slowly  with  temperature. 

This  is  a respectable  value  and  we  examined  related  materials  to 
see  if  any  further  improvements  could  be  obtained.  A number  of 
related  materials  were  prepared.  Many  were  clearly  centrosym- 


metric  at  room  temperature  for  example,  (LiRbSC>4,  LiTlS04, 
CsNH4S04,  KNH4SC>4);  others  either  formed  hydrates  (LiNaS04)  or 
proved  insoluble  (LiNH4Te04),  thus  precluding  the  growth  of 
crystals.  Only  LiKS04  and  LiCsS04  were  studied  in  detail. 


Our  investigation  of  these  two  materials,  reported  below,  shows 
them  to  be  essentially  different  from  LiNH4S04. 


TABLE  8.  OPTICAL  SECOND  HARMONIC  POWDER  ANALYSIS 
DATA  (PARTICLE  SIZE:  212-300  pm). 


Signal  Level 
Relative 
a-Si02 


Crystal 


Phase 


Temperature 
Interval  (K) 


Result 


We  conclude  that  further  study  of  related  materials  is  merited 
on  the  basis  of  the  high  f igure-of-merit  obtained  for  LAS  and 
DLAS.  However,  without  even  one  related  material  it  is  difficult 
to  determine  precisely  what  direction  should  be  taken.  The 
neutron  scattering  analysis  of  these  two  materials  already  in 
progress  might  reveal  enough  about  their  microscopic  structure 
to  suggest  plausible  directions. 

( 2 ) Low  Temperature  Phase  Transition 

On  cooling,  LiNH^S04  is  reported  to  exhibit  a first  order 
phase  transition  from  the  orthorhombic  (C2v-mm2)  structure  to 
monoclinic  at  10°C.  Yuzvak  et  al . (Ref.  3-3)  reported  a step-wise 
increase  in  spontaneous  polarization  at  this  phase  transition 
and  assigned  the  low  temperature  phase  to  the  polar  class  2.  A 
recent  study  by  Poulet  and  Mathieu  (Ref.  3-6)  confirmed  this 
polar  assignment  from  measurements  of  piezoelectric  response 
induced  by  stress  along  the  two-fold  axis  below  10°C.  However, 
Mitsui  et  al.  (Ref.  3-4)  concluded  from  the  temperature  behavior 
of  optical  second  harmonic  generation  that  the  low  temperature 
phase  was  centrosymmetric  class  2/m.  Aleksandrova  et  al. 

(Ref.  3-2)  also  concluded  on  the  basis  of  SHG  measurements  that 
the  low  temperature  phase  was  centrosymmetric. 

Our  pyroelectric  and  SHG  measurements  confirm  the  data  of  Yuzvak 
et  al.  (Ref.  3-3)  that  the  low  temperature  phase  of  LiNH4S04 
is  polar.  However,  our  pyroelectric  current  measurements 
indicate  that  the  phase  transition  is  accompanied  by  a decrease 
in  spontaneous  polarization.  Our  data  also  shows  a low  temper- 
ature polar  phase  for  LiND4S04. 
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3 . 2 Lithium  Cesium  Sulfate  and  Lithium  Potassium  Sulfate 

On  the  basis  of  optical  second  harmonic  generation  and  pyroelectric 
measurements  on  LiCsSO^,  we  resolve  the  ambiguity  in  the 
literature  about  the  space  group  of  this  material  and  conclude 
that  lithium  cesium  sulfate  does  not  belong  to  a pyroelectric 
space  group. 

Separate  measurements  of  the  pyroelectric  coefficient,  dielectric 

constant  and  specific  heat  of  LiKSO.  yield  a f igure-of-mer it 
3 2 4 

of  1.25  x 10  cm  /C . This  is  less  than  one-third  the  figure- 
of-merit  of  TGS.  We  concluded  that  this  material  does  not  have 
attractive  properties  for  practical  applications. 

3.2.1  Preparation  and  Characterization  of  Materials 

The  double  sulfates  LiKS04  and  LiCsS04  were  prepared  by 
combining  aqueous  solutions  of  the  appropriate  alkali  sulfates 
(Apache)  according  to  stoichiometry.  Single-crystals  were 
grown  from  aqueous  solution  via  isothermal  evaporation  at 
313  K.  Lattice  parameters  were  obtained  from  X-ray  precession 

O 

diagrams.  For  LiKSO.,  the  hexagonal  unit  cell  has  a = 5.15  A 

O ** 

and  c = 8.63  A,  and  for  LiCsSO.,  the  orthorhombic  unit  cell 

o o ^ o 

has  a = 9.45  A,  b = 5.45  A,  and  c = 8.80  A.  These  values 
are  in  good  agreement  with  the  previously  reported  work  of 
Bradley  (Ref.  3-7)  and  Hahn  et  al.  (Ref.  3-8). 

Flame  emission  analyses  showed  300  and  500  ppm  Na  impurity  in 
the  crystals  of  LiKS04  and  LiCsS04,  respectively.  Atomic 
absorption  analyses  showed  300  ppm  Rb  in  the  LiCsS04  crystals. 

3.2.2  Optical  Second  Harmonic  Generation  Measurements 

Powders  ground  from  single  crystal  samples  of  LiKS04  and  LiCsSO^ 
were  examined  by  optical  second  harmonic  analysis  (Ref.  3-9) 
from  260  to  560  K.  LiKS04  exhibited  a temperature  invariant 
(+  10%)  second  harmonic  signal  on  the  order  of  0.5  a-Si02. 

This  is  consistent  with  the  reported  non-centrosymmetr ic  crystal 
structure,  space  group  P6-j,  (Ref.  3-7)  and  indicates  the  absence 
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of  any  phase  transitions  in  the  temperature  range.  In  crystals 

_3 

of  LiCsSO^,  no  detectable  second  harmonic  at  a level  of  10 

ot-SiC>2  was  observed  from  260  to  560  K,  indicating  the  f 

absence  of  a non-centrosymmetric  phase  in  this  temperature 
interval . 

3.2.3  Specific  Heat  Measurements 

Single  crystals  of  LiKSO^  and  LiCsSO^  were  examined  from 
280  to  580  K with  a computer-controlled  Perkin-Elmer  dif- 
ferential scanning  calorimeter,  Model  DSC-2.  The  instrument 
calibration,  operating  parameters,  and  data  refinement  were 
performed  as  previously  described  (Ref.  3-10).  No  phase 
transitions  in  either  material  were  observed  in  agreement 
with  the  second  harmonic  generation  experiments. 

Two  specific  heat  determinations  were  made  from  300  to  520  K 
on  two  different  samples  of  each  material.  The  four  measure- 
ments on  each  material  agreed  within  1%.  The  resulting  mean 
values  are  reported  as  molar  heat  capacities  and  are  given 
in  Table  9 in  intervals  of  20  K.  Using  a least  squares 
curve  fit,  the  temperature  dependence  of  the  molar  heat 
capacity  of  the  temperature  interval  300  < T(K)  520  is  for 
LiKS04, 

c (J  mol-1K-1 ) = 81.1  + 1 . 247 ( 10-1  )T  + 7 . 2 ( 10-6 )T2 

hr 

and  for  LiCsSC>4, 

cp(J  mol-1K-1  ) = 63.0  + 1.585(10-1)T  + 4 . 8 ( 10-6  )T2 . 

The  c values  are  considered  accurate  to  + 1.5%  relative 
P 

to  the  NBS  values  for  sapphire  (Ref.  3-11). 

The  similarities  of  the  molar  heat  capacities  reflect 
similar  interatomic  bonding  energies  in  the  two  materials. 

The  temperature  slope  of  the  molar  heat  capacities  are  quite 
different,  namely,  for  LiKSC>4,  1.32  (10  *)  J mol  2, 
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TABLE  9 

ISOBARIC  MOLAR 

HEAT  CAPACITIES  OF 

LiKSO , 

4 

(MW  = 142.10)  and 

LiCsSO  . 

4 

(MW  = 235.' 

Cp  ( J/mol 

~ 1 K~ 1 ) 

T ( K ) 

LiKSO  . 

4 

LiCsSO 

300 

119.0  + 1 

. 5% 

111  . 0 + 

320 

121.7 

114.0 

340 

123.8 

116.9 

360 

126.9 

120.9 

380 

128.8 

123.6 

400 

130.7 

126.5 

420 

132.9 

128.5 

440 

135.3 

131.5 

460 

138.0 

135.7 

480 

142.0 

1 3 ~ - 7 

500 

145.8 

143.2 

520 

148 . 0 

146.9 

-1  -1  -2 

and  for  LiCsSC>4,  1.60(10  J mol  K over  the  temperature 
range  measured.  This  difference  can  be  attributed  to  dissimilar 
crystal  structures,  namely,  LiKS04  is  hexagonal,  space  group 
P6^  while  LiCsS04  is  orthorhombic,  probable  space  group  Pcmn. 

3.2.4  Pyroelectric  Measurements 

Pyroelectric  measurements  in  LiCsS04  were  performed  by  a modi- 
fication of  the  Chynoweth  technique  (Ref.  3-13).  In  crystals  of 
LiCsS04,  all  three  principal  crystallographic  directions  were 
studied  from  room  temperature  to  about  473  K.  No  pyroelectric 

response  was  observed.  This  result  places  an  upper  limit  on  the 

-10  2 

pyroelectric  coefficient  of  10  C/cm  K suggesting  that  the 
structure  is  non-polar. 


For  LiKSO^ , single  crystal  samples  oriented  with  surfaces  normal 

to  [001]  were  studied  by  the  Byer-Roundy  technique  (Ref.  3-13). 

The  sample  plates  were  polished  on  both  sides  to  a thickness  of 

0.66  mm  using  3 ym  alumina  grit  and  mineral  oil.  Both  faces  of 

2 

the  sample  (area  0.32  cm  ) were  fully  electroded  with  silver 

paint.  The  measurements  were  performed  under  a helium  atmosphere 

at  a heating  rate  of  6.5  K min-'*'.  The  pyroelectric  coefficient 

-9  -2  -1 

was  found  to  be  1.6  + 10%  x 10  C cm  K and  temperature  in- 
variant from  313  to  343  K.  The  measurement  above  343  K was 

impeded  by  a large  background  current.  Early  work  by  Ackermann 

_ g 

(Ref.  3-14)  showed  a pyroelectric  coefficient  of  1.61  x 10  C 
cm  2K  -*-  at  293  K which  rises  slowly  to  1.78  x 10  9C  cm  2K  1 
at  352  K,  in  reasonable  agreement  with  our  measurements. 

The  application  of  a dc  field  of  5 kV  cm-'*  for  1 hour  and 
2 kV/cm  ■*■  for  an  additional  12  hours  did  not  produce  a change 
in  sign  of  the  polarization  in  LiKSO^.  This  suggests  that  LiKSO^ 
is  not  ferroelectric  in  agreement  with  the  observations  of  Ando 
(Ref.  3-15). 

3.2.5  Dielectric  Measurements 

The  dielectric  constant  of  LiKSO.  was  measured  at  8 kHz.  For 

4 

the  dielectric  constant,  of  LiKSO^,  we  measure  a value 

of  7.2  + 10%  at  room  temperature  which  increases  to  a value  of 

7.6  + 10%  at  473  K.  Ando  (Ref.  3-15)  measured  a value  of  about  6 

for  £33  at  100  kHz  over  the  same  temperature  interval.  This 

small  discrepancy  may  be  attributed  to  the  dispersion  of  the 

dielectric  constant  in  LiKSO.. 

4 

3.2.6  Fiqure-of-Merit  and  Conclusions 

The  absence  of  a measureable  SHG  signal  and  pyroelectric  coef- 
ficient in  LiCsS04  indicates  that  the  correct  space  group  is 
Pcmn  and  not  Pc21n,  thus  resolving  any  ambiguity  (Ref.  3-8)  in 
the  determination  of  the  structure  of  this  material.  Since  this 
is  not  a polar  space  group  we  conclude  that  this  material  has  no 
pyroelectric  effect. 


3 

Combining  a density  of  2.39  gm/cm  for  LKS  obtained  from  the 
lattice  parameter  measurements  with  the  values  of  the  pyroelec- 
tric coefficient,  dielectric  constant  and  specific  heat  describ- 

3 2 

ed  above,  one  obtains  a value  of  1.25  x 10  cm  /C  for  the 
f igure-of-merit  of  LiKSO^  at  room  temperature.  This  value 
varies  only  slightly  from  room  temperature  up  to  70°C.  Since 
this  value  is  less  than  one-third  the  f igure-of-merit  of  TGS , we 
conclude  that  lithium  potassium  sulfate  does  not  have  attractive 
properties  for  practical  applications. 
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Lithium  Trihydrogen  Selenite 

Lithium  trihydrogen  selenite  LiH-^SeO^^  is  ferroelectric  at 

room  temperature  with  a very  large  value  of  spontaneous  polari- 

2 

zation  (~15  yC/cm  ) (Ref.  3-16).  Hydrostatic  pressure 
studies  (Ref.  3-17)  have  shown  the  Curie  point  to  be  145°C,  35° 
above  the  melting  point. 

Dielectric  hysteresis  measurements  were  made  on  several  samples 
of  LiH ^ ( SeO^ ) 2 using  a 60  Hz  Sawyer-Tower  bridge.  Sample 

to  sample  variations  were  very  large.  Measured  values  of  P 

. - 2 2 s 
varied  from  7.7  UC/cm  to  14  yC/cm  . In  addition,  some 

hysteresis  loops  also  showed  an  internal  bias  of  as  much  as  1.8 

kV/cm.  A hysteresis  loop  representative  of  one  of  the  better 

samples  is  shown  in  Figure  33. 
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Figure  33.  Ferroelectric 

hysteresis  loop  for  LiH^SeOj^- 

The  vertical  axis  is  1 yC/division 

and  the  horizontal  axis  is  200/V 

division.  For  samples  of  area 

A ^ = 0.227  cm2  and  thickness  t = 

0.50  mm,  Ps  = 12.1  yc/cm2  and 

E =3.8  kV/cm. 
c 
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Measured  pyroelectric  coefficients  in  the  lithium  trihydrogen 
selenite  samples  were  found  to  be  directly  related  to  the  degree 
of  poling  which  could  be  achieved.  These  results,  shown  in 
Table  10,  demonstrate  that  the  samples  can  have  "easy"  or  "hard" 
poling  directions. 
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TABLE  10.  COMPARISON  OF  EASY  AND  HARD  POLING 
DIRECTIONS  FOR  A LiH3(Se03)2  SAMPLE  AT  24°C. 


kV/cm 

Field  Direction  yC/cm 1 @ 60  Hz  yC/cm-K 


+ 

7.7 

3.2 

0.92  x 

IQ"8 

11.3 

5.0 

1.45  x 

IQ’8 

Pyroelectric  coefficient  and  dielectric  constant  were  measured 
on  a sample  that  showed  no  bias  effects;  the  sample  was  poled  at 
high  temperature  (84°C)  and  the  pyroelectric  data  taken  on 
cooling  without  a bias  field.  The  pyroelectric  and  dielectric 
data  are  shown  in  Figure  34.  These  data  suggest  that  LiH^SeO^^ 
might  be  approaching  a proper  phase  transition.  This  would  be 

4 

in  agreement  with  a Curie  constant  4 x 10  K measured  by 
Samara  and  Anderson  (Ref.  3-17). 

The  ratio  p/e  is  smaller  than  that  found  in  TGS-type  crystals. 

The  large  value  of  P does  not  produce  a correspondingly  large 
pyroelectric  current  since  P is  relatively  constant  with 

S 

temperature  in  the  range  of  the  data,  about  100°C  below  the 

Curie  point.  Over  the  temperature  range  40-80°C,  the  change  in 

2 

P (calculated  by  integrating  p)  is  1.35  yC/cm  , less  than 

s 2 
10%  of  P . TGS  by  comparison  has  a change  of  P of  1.2  yC/cm 

S o 

for  0-40°C,  nearly  35%  of  Pg. 

These  results  suggest  that  modifications  of  lithium  trihydrogen 
selenite  might  exhibit  useful  properties  for  practical  applica- 
tions . 
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4.  PROPIONATES 

Initial  measurements  of  the  pyroelectric  and  dielectric  proper- 
ties of  dicalcium  lead  propionate  (DLP)  reveal  the  presence  of 
an  internal  bias  field.  Annealing  samples  of  dicalcium  lead 
propionate  did  not  eliminate  the  effects  of  internal  bias  fields 
as  reported  in  the  literature;  indeed  the  properties  of  the 
samples  were  severly  degraded.  Initial  measurements  on  the 
specific  heat  of  DLP  and  dicalcium  strontium  propionate  (DSP) 
show,  as  expected,  no  large  anomaly  at  the  transition  temperature, 


A theoretical  analysis  based  on  a free  energy  appropriate  to  a 
second-order  pseudoproper  phase  transition  was  carried  out;  the 
predictions  agree  with  our  measurements  as  well  as  with  certain 
published  data.  In  this  analysis,  the  f igure-of-merit  has  a 
maximum  as  a function  of  temperature  which  can  be  expressed  in 
terms  of  the  free-energy  parameters.  Using  data  from  the 
literature  (which  is  compatible  with  our  measurements),  we  can 
estimate  these  free-energy  parameters  and  project  a maximum  fig- 
ure-of-merit  of  4.2  x 10^  cm2/C  for  DSP  and  7.3  x 10^  cm2/C 
for  DLP. 


Initial  Pyroelectric  and  Dielectric  Measurements  on  DLP 


A Byer-Roundy  measurement  on  an  as-grown  DLP  sample  was  perform- 
ed in  inert  gas  (helium).  The  values  of  the  pyroelectric 

-9  2 

coefficient  thus  determined  were  1.0  x 10  C/cm  K at  30°C, 
and  8.5  x 10  ^ C/cm2  K at  the  peak  (Fig.  35).  The  long  tail 
of  the  pyroelectric  current  above  Tc  - 56°C  is  probably  due 
to  the  internal  bias  field  which  normally  exists  in  as-grown 
DLP  samples  (Ref.  4-1). 


Using  the  dynamic  method  we  measured  the  temperature  dependence 
of  the  ratio  p/c  in  an  as-grown  sample  of  DLP,  applying  various 
dc  fields  in  the  direction  opposite  to  the  internal  bias 
(Fig.  36).  The  calibration  in  this  measurement  was  obtained  by 
setting  at  40°C:  p = 2 x 10~^  C/cm2  K,  c = 1.8  J/cm^  K. 

As  the  dc  field  was  increased  from  0 to  12  kV/cm,  the  peak  in 
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Figure  35.  Pyroelectric  coefficient  of  as-grown  DLP 
sample  as  measured  by  the  Byer-Roundy  technique. 
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p/c  became  sharper  and  shifted  slightly  to  lower  temperature. 

This  behavior  is  similar  to  the  pyroelectric  behavior  of  proper 
f erroelectrics  like  TGS  (Ref.  4-2)  and  colemanite  (Ref.  4-3) 
under  a biasing  field,  and  can  be  explained  in  a similar  way. 

The  temperature  dependence  of  e was  measured  at  800  Hz  during 
slow  heating  and  cooling  (Fig.  37).  The  heating  curve  normally 
shows  a lower  peak  than  the  cooling  curve.  The  effects  of 
external  dc  fields  applied  in  the  same  direction  as  the  internal 
bias  field  (+  3 kV/cm)  and  opposite  the  internal  bias  (-  3 kV/cm) 
are  shown  in  Figure  38. 
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COOLING  NO  BIAS 
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Figure  37.  Temperature  dependence  of  dielectric 
constant  of  as-grown  DLP. 
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Figure  38.  Temperature  dependence  of  dielectric 
constant  of  as-grown  DLP  under  dc  field  of  3 kV/cm 
applied  in  the  direction  of  the  internal  bias  (+) 
and  opposite  to  the  internal  bias  (-). 
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Annealing  Experiments  on  Dicalcium  Lead  Propionate 

In  an  attempt  to  eliminate  the  effect  of  internal  bias  fields  in 
dicalcium  lead  propionate,  samples  were  annealed  in  air  for  24 
hours  at  200°C.  After  annealing,  the  sample  had  a light  tan 
discoloration  compared  with  its  water-clear  initial  appearance. 

Measurements  of  the  dielectric  permittivity  of  the  annealed 
samples  yielded  a room-temperature  value  of  6.9,  in  agreement 
with  the  unannealed  samples.  However,  the  peak  at  the  transi- 
tion temperature  was  only  7.5  as  opposed  to  a peak  value  of 
12-14  observed  in  the  as-grown  samples.  The  dielectric  permit- 
tivity did  not  change  with  externally  applied  bias  fields  as 
large  as  5 kV/cm. 

Measurements  of  the  pyroelectric  coefficient  by  the  Byer-Roundy 
technique  revealed  no  detectable  pyroelectric  current  in  the 
annealed  samples. 

Thus,  the  annealing  procedures  did  not  simply  eliminate  the 
effect  of  internal  bias  fields  as  reported  in  the  literature. 
Rather,  the  properties  of  the  samples  were  seriously  degraded. 

Two  explanations  are  possible,  namely,  either  chemical  impurities 
in  the  samples  produced  the  observed  degradation  or  the  proper 
annealing  technique  was  not  fully  described  in  the  literature. 

We  lean  towards  the  latter  hypothesis. 

4 . 3 Initial  Measurements  of  Specific  Heat  of  Dicalcium 

Strontium  Propionate  and  Dicalcium  Lead  Propionate 

Initial  specific  heat  measurements  were  made  on  samples  of 
Ca2  Sr  (CH^CH^OO ) g and  Ca2  Pb  (CH^CH^OO)  ^ which  were  obtained 
from  outside  sources;  the  samples  were  visibly  flawed. 

Measurements  made  from  270  K to  360  K with  an  accuracy  of  about 
5%  did  not  reveal  any  anomalous  behavior  at  the  f erroelectric- 
paraelectric  phase  transition  temperature  (283  K for  DSP  and 
331  K for  DLP ) . This  is  consistent  with  the  expectation 
that  the  anomaly  in  the  specific  heat  is  quite  small  (1-4%)  for 
these  materials. 
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The  molar  heat  capacity  of  dicalcium  strontium  propionate  has 
the  form: 

cp  ( J/molK ) = - 375.9  + 3.10  T 290K  < T < 360K. 

Roughly  the  same  magnitude  with  a similar  linear  temperature 
dependence  was  obtained  for  the  lead  salt.  For  DSP  at  310  K, 
the  value  for  the  specific  heat  is  21%  lower  than  that  reported 
by  Nakamura,  et  al.  (Ref.  4-4). 

4 . 4 Analysis  of  the  Pyroelectric  Behavior  of  DLP  above 

According  to  Takashige  et  al.  (Ref.  4-1),  the  dielectric  constant 
of  unbiased  DLP  crystals  obeys  a modified  Curie-Weiss  law: 


e = c + c/(t  - t ) 

c o 


(1) 


with  ec  = 4 and  C = 69  K.  Based  on  this  result,  the  pyro- 
electric behavior  of  as-grown  DLP  above  T can  be  analyzed  in 

c 

the  following  way. 


From  a few  degrees  above  Tc  and  upwards,  where  nonlinearities 
in  the  dielectric  constant  can  be  neglected,  it  can  be  assumed 
that  the  polarization  induced  by  the  biasing  field  Eb  is, 

P = eo(E  - l)Eb  = co((ec  - D + C/(T  - V]Eb'  (2) 

-14 

where  eQ  = 8.85  x 10  F/cm,  is  expressed  in  V/cm  and 
P is  in  C/cm^.  The  pyroelectric  coefficient  is  thus, 

p = -dP/dT  = eoCEb/(T  - To)2.  (3) 


,-1/2 


It  follows  that  plotting  p versus  T one  should  obtain  a 


linear  plot  with  a slope  given  by  (eQCEb)  . Figure  39  shows 
obtained  by  the  Byer-Roundy  measurement  of  p,  as  a 
function  of  T.  The  agreement  with  the  expected  behavior  is 


reasonably  good.  The  value  obtained  for  E.,  using  C = 69  K, 


is  about  15  kV/cm  in  agreement  with  the  value  estimated  by 
Takashige  et  al.  (Ref.  4-1),  from  their  dielectric  measurements. 


J 


» 
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Both  dicalcium  lead  propionate  (DLP)  and  dicalcium  strontium 
propionate  (DSP)  exhibit  two  phase  transitions,  and  the  structur- 
es of  the  two  crystals  are  reported  to  be  identical  in  each  of 
the  corresponding  phases  (Ref.  4-1).  It  is  the  high-temperature 
transition  (occuring  at  54°C  and  8.5°C  in  DLP  and  DSP,  respective- 
ly) which  is  of  interest  here.  The  remarkable  resemblance, 
visible  in  the  measurements  reported  above  and  elsewhere  (Ref.  4-1), 
between  the  dielectric  behavior  of  the  propionates  at  these 
transitions  and  that  of  ordinary  ferroelectrics  at  second-order 
phase  transitions  has  led  to  the  suggestion  (Ref.  4-5)  that  these 
are  pseudoproper  ferroelectrics,  i.e.,  that  the  order  parameter 
has  the  same  symmetry  as  the  polarization  so  that  a linear 
coupling  between  the  two  is  possible  (Ref.  4-6).  If  in  fact  no 
increase  in  the  size  of  the  unit  cell  occurs  at  the  transition, 
this  view  is  consistent  with  the  observed  reduction  in  point 
group  symmetry  from  422  (D4)  to  4 (C4),  since  it  has  been 
shown  (Ref.  4-7)  that  this  symmetry  change  can  be  induced  by  the 
one-dimensional,  polar  representation  A2  of  422.  If  the  order 
parameter  n is  one-dimensional  and  the  transition  is  second 
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order  (so  that  terms  of  order  ^ and  higher  can  be  neglected), 
the  thermodynamic  potential  can  be  written  in  the  form: 

a (T-T  ) -1 

*(n,P)  = * (T)  + — ° n2  + A n4  + p2 

2 2 (4; 

+ ar|P  - EP 

where  a , g,  x ^ and  a are  taken,  as  usual,  to  be  temper- 
ature independent.  In  the  discussion  of  this  thermodynamic 
potential  it  is  necessary  to  avoid  the  restriction  E = 0,  since 
the  experimental  results  for  unannealed  DLP  show  strong  effects 
due  to  internal  bias  fields. 

The  equilibrium  values  of  n(E,T)  and  P(E,T)  are  determined 
from  the  equations: 

8$/3n  = ao(T-TQ)n  + 6n3  + aP  = 0 

(5] 

3$/3P  = X P + an  - E = 0 

Thus,  the  temperature  and  field  dependence  of  the  polarization 
is  readily  given  in  terms  of  the  temperature  and  field  depend- 
ence of  the  order  parameter, 

P = XE  - axn (E,T)  (6 ; 

while  n(E,T)  is  in  turn  obtained  from  the  solution  of  the  cubic 
equation: 


{fj  »«■ 


T)  + (T-T  )n(E,T)  + (C/ax)E  = 0 


where  we  have  introduced  the  notation, 

T = Tq  + a2x/aQ,  K = ax(ao/6)1//2  and  C = a2x2/aQ 

In  terms  of  the  function  *;(E,T)  determined  by  this  equation, 
one  can  obtain  the  dielectric  constant. 


e(E,T)  = 1 + x + 


[T-Tj.  + 3(ax/K)  n (E,T)] 
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the  pyroelectric  coefficient, 


p(E,T) 


axn (E,T, ) 

[T-T1  + 3(aX/K)2n2(E,T)] 


(9) 


and  the  specific  heat. 


T(ax)2/C  n2 (E,T) 

[T-Tx  + 3(ax/K)2n2(E,T) ] 


(10) 


For  the  special  case  E = 0 when  all  internal  bias  fields  have 
been  annealed  out.  Equation  (7)  can  be  easily  solved  for  n(0,T); 
and  we  obtain  in  that  case  for  the  spontaneous  polarization. 


p 

s 


1K(T1-T)1/2 
0 


T £ 

T > T . 


(ID 


By  least-square  fitting  of  the  P data  of  Takashige  et  al. 

s 

(Ref.  4-1)  between  37  and  49°C,  we  found  for  K the  value  K = 
—8  2 1 /2 

3.1  x 10  C/cm  K ' . Similiarly  where  E = 0 the  dielectric 
constant  takes  the  form. 


e = 


1+X 

1+X 


+ C/CT-T^ 

+ C/2 (T^-T) 


T < T1 
T > T1  . 


(12) 


The  dielectric  measurements  of  Takashige  et  al.  (Ref.  4-1) 
confirm  that  the  dielectric  constant  of  unbiased  DLP  crystals 
has  this  form  with  1 + x = 4 and  C - 69  K. 


With  C and  K thus  determined,  one  has  enough  information  to 
theoretically  determine  the  pyroelectric  coefficient. 


K/2  (T1~T) 


-1/2 


P = 


and  the  specific  heat. 


T < Tx 
T > T . 


(13) 
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c + K2T/2C 
o 

c 

o 

As  expected  (Ref.  4-6),  the  results  obtained  here  are  identical 
to  those  for  a proper  ferroelectric  with  a second-order  transi- 
tion (Ref.  4-3).  Only  the  anomalously  small  Curie  constant, 

2 2 

C = f xQ/3f  produced  by  the  weak  coupling  between  n and  P, 
serves  to  demonstrate  the  "almost  non-polar"  character  of  the 
structural  instability  and  to  provide  a quantitative  meaning  to 
the  distinction  between  proper  and  pseudoproper  transitions.  No 
qualitative  distinction  based  on  symmetry  is  possible  (Ref.  4-6). 

4 . 6 Projected  Figure  of  Merit 

While  the  comparison  of  our  experimental  results  with  the 
predictions  given  above  is  still  in  progress,  one  can  examine 
the  consequences  of  this  analysis  for  the  f igure-of-merit  M(l)  = 
p/ec  in  unbiased  crystals.  Using  the  expressions  in  Equations 
(12),  (13)  and  (14)  above,  we  can  express  the  f igure-of-merit 
as , 

M (1 ) = (K/2)  (T^-T)"172  [(1+X  + (C/2)  (^-T)"1]"1  (cQ  + K2T/2C)  1 (15) 

The  maximum  of  M(l)  as  a function  of  temperature  occurs  at  a 
temperature  T , such  that, 

(TrTm)  - c/e  - (K2/cC)  (T^TJ2  = 0.  (16) 

As  we  shall  show  shortly,  the  third  term  in  this  equation  is 
negligible  in  the  propionates,  so  Tm  is  simply  given  by, 

Tm  = T0  - C/[2  (1  + X))  (17) 

In  this  approximation  we  can  also  take  c = cq  when  evaluating 
M ( 1 ) at  T = T . The  resulting  maximum  value  of  M(l)  is. 


T < T. 


T > T. 


(14) 
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M 

max 


(18) 


= (K/2cq)  (2(1+X)C>  7 

Note  that  M has  no  natural  optimum  with  respect  to  the 
max  c 

parameters  which  characterize  the  phase  transition:  increasing 
X or  decreasing  C will  always  increase  M 

max 

Table  11  summaries  the  experimental  values  (Refs.  4-1,  4-8,  4-9) 

of  the  parameters  required  to  evaluate  M and  to  justify  the 

in  ci  x 

approximation  which  led  to  Equation  (18).  The  results  of  the 

2 

calculation  are  given  in  Table  12.  The  ratio  K /2C  is  a 


TABLE  11.  EXPERIMENTAL  VALUES  USED  IN 
CALCULATION  OF  M 

max 


K(10~8  C/cm2K1/2) 
COO 


T1 

(°C)  a 

1 + 

■ X 

c 

o 

(J/cm3K) 

(a) 

Ref.  4-1 

(b) 

Ref.  4-9 

(c) 

Ref.  4-8 

DLP 

DSP 

3.1 

5.6 

69  a 

73.2 

54 

8.5 

4 3 

4 

1.8 

1.8 

TABLE  12.  PROJECTED  PARAMETERS  FOR 
PROPIONATES  . 


DLP 

K2/2C  (10-5  J/cm3K2)  5.5 

c-c  (10~2  J/cm3K,  at  T=T  ) 1.8 

o u 

T.-T  (K)  8-6 

l m 

M (103  cm2/C)  4.2 

max 


DSP 

24 

6.8 


9.1 


7.3 
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measure  of  the  strength  of  the  specific  heat  anomaly  c-cq  and 
we  find  it  to  be  very  small.  In  fact  the  maximum  value  of 
c-co  found  at  Tc  is  just  1%  of  cQ  in  DLP  and  less  than  4% 
of  cQ  in  DSP.  This  is  consistent  with  the  results  of  Nakamura 
et  al.  (Ref.  4-8)  for  DSP  and  with  our  own  preliminary  measure- 
ments. Clearly,  we  are  justified  in  neglecting  the  third  term 
on  the  right-hand  side  of  Equation  (16)  and  in  using  Equation 

(18)  for  M . The  next  entries  in  Table  12  show  that  T 

max  m 

occurs  close  enough  to  Tq  (roughly  9 K below  in  both  materials) 
so  that  having  neglected  higher-order  terms  in  the  thermodynamic 
potential  probably  has  not  seriously  affected  our  results. 


Finally,  the  values  obtained  for  = are  quite  substantial, 

3 3 7ax 

namely,  4.2  x 10  and  7.3  x 10  cm  /C  for  DLP  and  DSP, 
respectively.  Note  for  comparison  the  value  reported  earlier 
for  triglycine  sulfate  (TGS)  at  40°C  of  slightly  less  than 
4 x 103  cm2/C . 
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5. 


BORACITES 


Initial  dielectric  and  pyroelectric  measurements  were  performed 
on  samples  of  boracites  obtained  from  Plessey  Company.  Major 
difficulties  in  fully  poling  these  samples  were  encounter- 
ed. The  measurements  are  consistent  with  the  data  reported  by 
Plessey  Company  and  with  the  resulting  expectations  of  a high 
f igure-of-merit  for  Ni-Br  boracite. 

Specific  heat  measurements  were  made  on  samples  of  Zn-Br,  Cr-Cl, 
Fe-I,  Ni-Br  and  Cu-Cl  boracite.  In  addition,  infrared  vibration- 
al spectra,  electronic  spectra,  and  second  harmonic  measurements 
were  performed  on  some  of  the  samples  in  order  to  better  char- 
acterize the  material.  A number  of  initial  observations  can  be 
made : 


•The  cubic-orthorhombic  phase  transition  appears  to  be 
strongly  first  order.  This  is  indicated  by  the  large 
latent  heat  of  the  transition,  a detectable  thermal  hyster- 
esis, and  the  discontinuous  jump  in  the  second  harmonic 
intensity  at  the  transition. 

• The  different  members  of  the  boracite  family  may  be  struc- 
turally dissimilar.  This  is  indicated  by  large  differences 
in  the  molar  heat  capacity  of  Fe-I  and  Cr-Cl  boracite  in 
the  cubic  phase.  Large  differences  in  the  molar  latent 
heats  also  support  this  conclusion. 


5. 1 Initial  Dielectric  and  Pyroelectric  Measurements 
on  Boracites 

5.1.1  Poling  of  Samples 

figure  40  shows  an  SEM  image  of  a sample  of  Ni-Br  boracite 
fabricated  by  the  Plessey  Company.  The  electrodes  are  silver 
epoxy.  Measurements  made  on  such  samples  show  them  to  be  very 
difficult  to  fully  pole,  as  has  been  reported.  On  cooling 


through  the  transition  temperature  with  a bias  field  as  large  as 
4.3  kV/cm,  the  pyroelectric  current  contains  many  large  spikes 
characteristic  of  domain  switching  even  a few  degrees  below  the 
transition  temperature.  When  such  a bias  field  is  removed  well 
below  the  transition,  the  sample  would  on  occassion  switch  into 
a polarization  state  opposite  to  that  impressed  by  the  bias 
field.  These  difficulties  in  poling  confirm  the  experiment  at 


Figure  40.  An  electroded  sample  of  Ni-Br  boracite 
fabricated  by  the  Plessey  Co. 

Plessey  where  50  Hz  hysteresis  loops  could  not  be  obtained  with 
fields  as  large  as  30  kV/cm  in  Ni-Br  boracite.  Indeed,  even 
after  resorting  to  a high  voltage  pulse  poling  technique,  they 
observe  coercive  fields  between  50  kV/cm  and  80  kV/cm  only  5 °C 
below  the  transition  temperature. 

In  view  of  the  previous  experience  with  other  improper  ferro- 
electrics  which  are  also  ferroelastic,  (see  Section  2)  it  seems 
probable  these  poling  difficulties  are  associated  with  the 
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mechanical  boundary  conditions.  Thus,  both  growth-induced 
stresses  and  the  effect  of  the  electrodes  play  major  roles  in 
determining  the  orientation  of  polarization. 


The  application  of  a bias  field  should,  nevertheless,  improve 
the  state  of  poling  by  properly  orienting  those  portions  of  the 
sample  not  dominated  by  mechanical  effects.  Such  bias  fields 
were  accordingly  used  in  the  measurements  of  the  dielectric 
permittivity.  Unfortunately,  in  the  measurements  of  pyro- 
electric current,  the  leakage  currents  produced  by  the  bias 
fields  made  reliable  measurements  on  biased  samples  unattainable. 

5.1.2  Ni-Br  Boracite;  Initial  Pyroelectric  Measurements 

A typical  measurement  of  the  pyroelectric  coefficient  of  an 
unbiased  sample  of  Ni-Br  boracite  is  shown  in  Figure  41. 

' The  measurements  were  made  during  cooling  with  no  external  bias 

' 

field  applied.  The  poling  state  is  that  "spontaneously"  assumed 
by  the  sample  and  is  probably  determined  by  the  mechanical 
boundary  conditions.  The  temperature  region  from  123°C  to  the 
transition  temperature  at  126°C  contains  many  large  spikes 
indicative  of  an  incompletely  poled  sample.  The  significance  of 
this  graph  is  not  the  absolute  magnitude  of  the  pyroelectric 
coefficient  but  rather  the  substantial  increase  in  pyroelectric 
coefficient  that  occurs  as  the  transition  temperature  is  ap- 
proached. The  value  at  123°C  (3°  below  the  transition  temper- 
ature) is  about  six  times  the  value  at  90°C. 

An  estimate  of  the  degree  of  poling  can  be  obtained  by  comparing 

_o  2 

the  approximately  constant  value  of  0.26  x 10  C/cm  K (the 

pyroelectric  coefficient  in  the  temperature  range  of  88°C-100°C) 

— 8 2 

to  the  constant  value  of  1 x 10  C/cm  K reported  by  Plessey 
for  Ni-Br  boracite  in  the  10°C-40°C  temperature  range.  If  the 
enhancement  by  a factor  of  six  which  we  observe  near  the 
transition  is  superimposed  on  the  low  temperature  value  for  the 
pyroelectric  coefficient  measured  at  Plessey,  a f igure-of-merit 
substantially  better  than  that  of  DTGFB  would  be  attained  by 
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Figure  41.  Pyroelectric  coefficient  vs.  temperature  for 
Ni-Br  boracite. 
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Ni-Br  boracite  as  much  as  3 K below  its  transition  temperature. 
It  is  to  be  emphasized  that  this  is  projection,  not  a measure- 
ment. 


5.1.3  Ni-Br  Boracite;  Initial  Dielectric  Measurements 

The  dielectric  permittivity  of  Ni-Br  boracite  as  a function  of 

temperature  is  given  in  Figure  42.  The  measurement  was  made  at 

1 MHz  using  a very  low  test  voltage  (10  mV)  and  an  applied  bias 

field  of  20  kV/cm.  The  permittivity  varies  only  slightly  with 

temperature  except  for  the  small  upwards  step  that  occurs  at  the 

transition  temperature.  This  temperature  dependence  is  at 

variance  with  the  data  of  Plessey  which  shows  a smooth  increase 

up  to  the  peak  at  the  transition  temperature.  Differences  in 

the  absolute  magnitude  of  the  permittivity  as  well  as  the 

temperature  of  the  peak  also  occur.  These  discrepancies  can 

have  several  sources.  The  small  size  of  the  sample  (thickness  = 

-3  2 

260  ym;  area  = 3.6  x 10  cm  ; see  Fig.  40)  gives  a capaci- 
tance of  only  0.3  pF  at  its  maximum;  such  small  values  are 
difficult  to  measure.  Also,  the  measurements  were  made  under 
different  conditions.  First,  the  Plessey  data  was  taken  at 
10  kHz,  while  ours  was  taken  at  1 MHz.  Second,  the  character- 
istic test  voltage  of  Plessey's  General  Radio  1615-A  bridge  is 
about  1 V;  while  our  test  voltage  is  0.01  V.  Third,  the  Plessey 
data  were  taken  on  a presumably  unbiased  sample,  while  we  used  a 
20  kV/cm  bias  field.  Thus,  these  discrepanices  are  not  regarded 
as  significant. 

5.1.4  Cu-Cl  Boracite:  Initial  Dielectric  Measurements 

An  initial  measurement  of  the  temperature  dependence  of  the 
dielectric  permittivity  for  an  unbiased  sample  of  Cu-Cl  boracite 
is  shown  in  Figure  43.  The  measurements  were  made  at  1 MHz 
using  a probe  field  of  10  mV  with  no  bias  applied  to  the  sample. 
The  measured  values  agree  with  the  Plessey  report  within  the 
estimated  experimental  accuracy. 


CuCl  -BOR AC  I TE 


TEMPERATURE 


5 . 2 Specific  Heat  Measurements 


Calorimetric  data  were  obtained  with  a Perkin-Elmer  differ- 
ential scanning  calorimeter  Model  DSC-2  interfaced  to  a Tektronix 
programmable  calculator  Model  31.  The  system  is  equipped  with 
an  Autoscan  zero  mini-computer  for  automatic  baseline  corrections. 

Single  crystal  samples,  typically  between  5 and  50  mg,  were 

- 4 

weighed  with  a precision  of  1 (10  ) and  run  under  a nitrogen 

atmosphere.  In  all  cases,  open  aluminum  sample  pans  were 
employed  with  the  largest  crystal  facet  parallel  to  the  bottom 
of  the  pan  to  assure  good  thermal  contact.  The  instrument 
was  calibrated  relative  to  the  melting  point  and  enthalpy  of 

fusion  of  5N  indium  metal;  Tu  = 429.87  K and  AH*  = 3266.8  J mol  1 

M f -1 
and  of  6N  lead  metal;  TM  = 600.65  K and  AH^  = 4768.1  J mol 

Both  the  melting  point  and  enthalpy  of  fusion  were  found  to  be 

independent  of  heating  and  cooling  rates  from  0.625  to  20  K min-'*' 

and  reproducible  to  + 0.02%  and  + 0.75%,  respectively.  Initial 

specific  heat  determinations  were  done  at  5 K min  * over  a 

large  temperature  range  in  each  boracite  (typically  200  K) 

to  obtain  the  temperature  dependence  of  the  heat  capacity.  The 

integration  intervals  was  0.4  K.  Detailed  specific  heat  data 

were  obtained  at  integration  intervals  of  0.078  K at  heating 

rates  of  1.25  K min  * in  the  vicinity  of  the  phase  transition. 

The  raw  data  was  further  refined  by  normalization  to  the  NBS 

values  of  sapphire  using  a sapphire  standard  run  under  the  same 

operating  conditions. 

The  temperature  dependence  of  the  heat  capacity  near  the  43m  - 
mm2  phase  transition  of  some  of  the  boracites  is  shown  in 
Figures  44  - 48.  These  curves  were  obtained  in  a heating  mode 
only  and  are  the  result  of  averaging  at  least  two  runs  on  at 
least  two  different  samples.  The  standard  deviation  of  the 
specific  heat  for  each  data  point  is  on  the  order  of  4%.  The 
transition  temperature  Tc  is  taken  to  be  the  peak  in  the  heat 
capacity.  The  peak  in  the  heat  capacity  at  the  transition  was 
found  to  be  asymptotic  as  a function  of  decreasing  the  heating 
rate  at  a constant  integration  interval.  For  example,  in  Zn-Br 
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Figure  44.  Isobaric  molar  heat  capacity  of  Cr-Cl  boracite 
in  the  vicinity  of  the  43m-mm2  phase  transition. 


* 


Figure  47.  Isobaric  molar  heat  capacity  of  Ni-Br  boracite 
in  the  vicinity  of  the  43m-mm2  phase  transition. 
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boracite  it  is  1660  J mol-1  K-1  at  5 K min-1  and  2406  + 

8 J mol  1 K * at  both  1.25  K min  * and  0.62  K min-'*'.  Similar  re- 
sults were  obtained  for  the  Fe-I  boracite,  and  this  behavior  is 
assumed  to  occur  in  the  other  boracites.  Consequently,  specific 
heat  curves  were  recorded  at  1.25  K min  ^ in  order  to  minimize 
any  instrumental  broadening  of  the  specific  heat  peak  in  the 
vicinity  of  the  transition. 

Table  13  is  a summary  of  the  thermodynamic  values  associated  with 
the  43m  - mm2  phase  transitions  in  the  boracites.  The 
enthalpies  of  the  phase  transitions  were  obtained  by  direct 
computer  integration  on  cooling  and  heating  and  by  graphical 
integration  on  heating  only  of  the  heat  capacity  curves.  All 
three  AH  values  were  in  agreement  to  within  + 3%. 

The  integration  is  performed  by  extrapolating  a line  from  the 

heat  capacity  curve  in  the  low  temperature  phase  to  a line  in 

the  high  temperature  phase.  This  is  done  over  a temperature 

interval  as  small  as  6 K in  Fe-I  and  as  large  as  17  K in  Ni-Br. 

The  entropy  of  the  transition  is  calculated  by  dividing  AH  by 

T . This  assumes  that  the  transition  occurs  over  a small 
c 

enough  temperature  interval  so  that  the  error  in  this  approxi- 
mation is  minimal.  Thus  the  values  of  AH  and  AS  are  upper 
bounds  on  the  true  latent  heat  and  entropy  discontinuity  at  Tc, 
as  they  contain  contributions  also  from  the  rise  in  the  heat 
capacity  just  below  the  transition  temperature. 

Several  conclusions  can  be  made  from  these  specific  heat  curves 
and  computed  thermodynamic  values.  The  sharp  peak  in  the  specific 
heat,  the  large  latent  heat  of  transition,  and  the  thermal 
hysteresis,  although  somewhat  small,  clearly  show  that  the  43m  - 
mm2  phase  transition  is  first  order.  The  anomalously  large  AH 
values  of  the  Zn-Br  boracite  relative  to  the  other  boracites  may 
reflect  the  existence  of  boracite  polytypes. 
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TABLE  13:  THERMODYNAMIC  VALUES  FOR  THE  43m-mm2  PHASE 

TRANSITION  IN  BORACITE  CRYSTALS. 


Boracite 

TC(K) 

AT  (K) 

AH(J  mol.  ) 

A S ( J mol . K 

Cr-Cl 

265.8+0. 2 

1.6+0. 4 

1966+14 

7.40+0.09 

Fe-I 

346.2+0.2 

3. 4+0. 4 

962  +7 

2.78+0.04 

Cu-Cl 

369.7+0.2 

3. 7+0. 4 

1211  +7 

3 . 28+0 . 05 

Ni-Br 

398.4+0.2 

2. 1+0. 4 

2950+20 

7.41+0.09 

Zn-Br 

584.6+0.2 

3. 2+0. 4 

4050+20 

6.95+0.04 

Schmid  et  al.  (Ref.  5-1)  have  previously  reported  the  heat  capacity 
of  Fe-I  boracite  measured  also  by  differential  scanning  calori- 
metry. Their  results  are  generally  in  good  agreement  with 
ours  except  for  a serious  discrepancy  concerning  the  latent  heat 
of  the  transition.  For  Fe-I  boracite  Schmid  et  al.  (Ref.  5-1) 
measure  AH  = 3100  + 300  J mol  1 whereas  we  obtain  AH  = 

962  + 7 J mol"*'*'.  Aside  from  errors  in  measurement  or  calcula- 
tion, this  discrepancy  may  be  explained  by  assuming  either 
different  materials  were  measured  or  that  polytypes  of  Fe-I 
boracite  exist.  This  is  in  agreement  with  the  suggestion 
of  Nelmes  (Ref.  5-2)  that  at  least  three  distinguishable  ortho- 
rhombic polytypes  exist. 

During  the  course  of  this  investigation  some  boracite  samples 
were  observed  to  exhibit  multiple  peaking  at  the  phase  transition. 
An  example  of  multiple  peaking  in  the  heat  capacity  of  Cr-Cl  and 
Ni-Br  boracites  are  shown  in  Figures  49  and  50.  The  general  be- 
havior was  reproducible  in  any  given  sample  with  small  fluctuations 
occurring  in  the  varying  specific  heat  peaks.  More  importantly 
the  AH  values  of  each  successive  determination  were  in  agreement 
to  with  + 3%.  In  the  case  of  Cr-Cl  and  Ni-Br  boracite  no  single 
peak  samples  were  observed  whereas  Zn-Br  and  Fe-I  boracite  were 


144 


more  often  single  peaked  than  multiple  peaked.  No  multiple 
peaked  sample  of  Cu-Cl  boracite  was  observed. 

Tables  14  to  18  show  the  measured  isobaric  (1  atm)  heat  capacit- 
ies in  5 K intervals  for  the  boracites  in  the  mm2  and  T3m 
phases.  Using  a least-squares  curve  fit  we  derive  the  temper- 
ature dependence  of  the  specific  heat  for  boracites  in  their 
appropriate  crystal  phase  over  a given  temperature  range.  The 
resulting  functions  are  shown  in  Tables  19  and  20. 

It  appears  from  this  data  that  either  the  divalent  metal  and/or 
the  halide  contribute  significantly  to  lattice  energies  of 
boracites . 


TABLE  14.  MOLAR  HEAT  CAPACITY  AT  CON- 
STANT PRESSURE  OF  Cr-Cl  BORACITE  IN  THE 
mm2  AND  43m  PHASES  (MW  = 475.10). 
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TABLE  18.  MOLAR  HEAT  CAPACITY  AT  CON- 
STANT PRESSURE  OF  Zn-Br  BORACITE  IN  THE 
mm2  AND  43m  PHASES  (MW  = 559.70). 


T(K) 

C (J  mol"1K'1) 

P 

T (K) 

C„(J  nol_1K_1) 

P 

455.0 

401.7  ± 1.4 

565.0 

566.8 

460.0 

409.2 

570.0 

574.2 

465.0 

416.7 

575.0 

581.7 

470.0 

424.5 

584.6 

T3m-mm2  transition 

475.0 

431.8 

590.0 

585.2 

480.0 

439.  2 

595.0 

588.5 

485.0 

436.  5 

600.0 

592.1 

490.0 

454.2 

605.0 

595.3 

495.0 

461.8 

610.0 

599.0 

500.0 

469.2 

615.0 

602.1 

505.0 

476.  8 

620.0 

605.6 

510.0 

484.) 

625.0 

609.0 

515.0 

491.5 

630.0 

612.5 

520.0 

499.0 

635.0 

615.5 

525.0 

506.7 

640.0 

619.0 

530.0 

514.5 

645.0 

622.5 

535.0 

521.7 

650.0 

625.9 

540.0 

529.0 

545.0 

536.7 

1 

544.0 

551.7 

560.0 

559.2 

8 


TABLE  19.  LEAST  SQUARES  FIT  MOLAR  HEAT  CAPACITY  FUNCTIONS 
FOR  BORACITES  IN  THE  mm2  PHASE. 


Boracite 


Cr-Cl 

Fe-I 

Cu-Cl 

Ni-Br 

Zn-Br 


Temperature 

Interval 

Cp(J 

mol 

-1) 

240  < 

T (K) 

< 

260 

cp  = -3.097 

(102) 

+ 

24.0 

(10_1) 

T 

270  < 

T (K) 

< 

340 

Cp  = -0-^7 

(102) 

+ 

13.0 

(10-1) 

T 

270  < 

T (K) 

< 

355 

Cp  = -3-776 

(102) 

+ 

21.6 

dO"1) 

T 

270  < 

T (K) 

< 

385 

Cp  = -1.986 

(102) 

+ 

16.2 

do"1) 

T 

455  < 

T(K) 

< 

575 



Cp  = -l.R«9 

(102) 

+ 

15.0 

no'1) 

T 

TABLE  20.  LEAST  SQUARES  FIT  MOLAR  HEAT  CAPACITY  FUNCTIONS 
FOR  BORACITES  IN  THE  43m  PHASE. 


Boracite 


Temperature 

Interval 


Cr-Cl 

270  < T (K)  £ 400 

Fe-I 

350  < T (K)  < 440 

Cu-Cl 

385  < T (K)  < 450 

Ni-Br 

410  < T ( K ) < 470 

Zn-Be 

590  _ T (K)  _ 650 

1 

1 

(J  mol 

-1-1 

C 

P 

K i) 

Cp  = 1-374 

(102) 

+ 4.588 

dO"1) 

T 

Cp  = 1-1«4 

(102) 

+ 7.353 

(10"1) 

T 

Cp  = 1-741 

(102) 

+ 5.654 

(10"1) 

T 

cp  = 2-*ii 

(102) 

+ 4.279 

(10'1) 

T 

C = 1.849 

p 

(102) 

+ 6.784 

(10"1) 

T 

5.3 


Infrared  Vibrational  Spectra 


Solid-state  infrared  vibrational  spectra  were  recorded  on  a 
Perkin-Elmer  infrared  spectrophotometer,  Model  283,  with  AgCl 
optics.  Scanning  time  was  150  cm  * min  ^ at  a response  time 
of  2 sec  using  slit  program  7.  Single  crystal  samples  ground 
with  an  agate  mortar  and  pestle  were  run  as  mulls  containing 
approximately  12  mg  of  material.  Fluorolube  mulls  were  examined 
from  4000  to  13000  cm  while  nujol  mulls  were  employed  from 
1300  to  400  cm  The  sample  temperature  in  front  of  the  beam 
was  found  to  vary  from  36  to  38°C.  Therefore,  all  samples  are 
noted  as  being  recorded  at  37  + 1 °C , neglecting  temperature 
fluctuations  due  to  black  body  radiation. 

The  infrared  vibrational  spectra  in  the  solid-state  from  1600  to 
400  cm  1 are  shown  in  Figures  51  - 56  for  the  boracites.  The 
main  absorption  maxima  are  given  in  Table  21.  From  5000  to  1600 
cm  no  absorptions  were  observed  except  for  probable  combina- 
tions and  overtones.  More  importantly  no  absorption  which  could 
be  attributed  to  OH  was  found.  It  is  concluded  that  (with  the 
possible  exception  of  CuCl  which  has  not  been  studied  above 
1300  cm  ^ ) no  boracite  contained  hydrated  water  or  hydroxyl- 
substituted  halogens  in  amounts  greater  than  0.5%  by  weight. 

General  inspection  of  the  spectra  reveal  some  interesting  facts. 
The  spectrum  of  Cr-Cl  is  taken  in  the  ?3m  while  the  other 
spectra  are  taken  in  the  orthorhombic  phase.  It  appears 
that  the  spectra  of  Cu-Cl  and  Ni-Br  are  the  same,  forming  one 
type  while  Fe-I  and  Zn-Br  exhibit  two  other  types  of  spectra. 
Interestingly,  the  spectra  of  "dark  green"  and  "light  green" 
Cu-Cl  boracite  were  identical  within  our  resolution. 

The  absorption  at  1401  and  1348  cm  ^ in  Zn-Br  and  at  1385  and 
1338  cm  in  Ni-Br  are  attributed  to  the  degenerate  stretch 
for  BO-  with  D3h  site  symmetry.  The  latter  absorption  is 
for  B*  while  the  former  is  for  B*1  on  the  approximate  ratio 
of  1:4  corresponding  to  the  natural  isotope  ratio  of  18:82. 
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Figure  51.  Solid-state  infrared  vibrational 
spectrum  of  43m  Cr^B^^Cl. 


Figure  52.  Solid-state  infrared  vibrational 
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TABLE  21.  INFRARED_|BSORPTIONS  OF  BORACITES 
(FUNDAMENTALS  IN  cm"1). 


Fe-I 

Zn-Br 

Ni-Br 

Cu-Cl 

Cr-Cl 

1401 

1385 

— 

1348 

1338 

— 

— 

1261 

1262 

1261 

— 

— 

1169 

1172 

1170 

1175 

1149 

1150 

1158 

— 

— 

— 

1047 

1058 

1062 

1055 

970 

995 

999 

1000 

976 

— 

930 

941 

943 

— 

862 

877 

882 

877 

874 

723 

704 

722 

708 

718 

— 

659 

651 

654 

681 

607 

612 

619 

617 

619 

423 

443 

458 

462 

457 

5 . 4 Electronic  Spectra 

Solid-state  electronic  spectra  were  recorded  on  a Beckman  UV-Vis 
spectrophotometer  model  DK-2  scanned  at  2000  cm-1  min-^. 

Single  crystal  samples  on  the  order  of  10  mg  were  ground  to  a 
fine  powder  with  an  agate  mortar  and  pestle  and  run  as  nujol 
mulls  on  a filter  paper  support. 

The  divalent  metal  ion  of  boracites  in  both  the  43m  and  mm2 
crystal  phases  has  true  D4h  symmetry  but  will  be  approximated 
by  a trans  pseudo-octahedral  environment.  This  is  due  to  a 
broadening  of  the  widths  of  bands  involving  transitions  to  T and 
E terms . 

For  cubic  Cr-Cl  (Fig.  56)  we  observe  an  asymmetric  broad  absorp- 
tion centered  at  14000  cm-1  which  can  be  attributed  to  the 

5E„  — > '’t,.  transition  in  the  Cr++  ion. 
g 2g 
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Figure  57  shows  the  electronic  absorption  spectra  from  two 

different  samples  of  orthorhombic  Cu-Cl  boracite  denoted  "light 

green"  and  "dark  green".  Both  of  these  samples  exhibit  a broad 

absorption  band  centered  in  the  12000  cm-'''  to  13000  cm*'  region 

2 2 

which  can  be  attributed  to  the  — > T„  transition  in  the 

+ + g 2g 

Cu  ion.  This  broad  absorption  band  in  the  "light  green"  sample 

is  shifted  about  650  cm  ^ towards  higher  energies  than  the 

corresponding  absoprtion  in  the  "dark  green"  sample.  This  shift 

is  outside  the  range  of  experimental  uncertainties  and  is  much 

more  than  one  would  expect  from  impurities  or  mild  distortions. 

We  take  this  shift  to  be  evidence  for  the  existance  of  polytypes 

or  substantial  variations  in  stoichiometry  in  Cu-Cl  boracite. 

Finally,  it  should  be  noted  that  no  evidence  for  multiple 
oxidation  states  in  either  crystal  was  observed. 

5 . 5 Second  Harmonic  Generation  in  Ni-Br  Boracite 

Optical  second  harmonic  generation  (SHG ) was  performed  with  a 
second  harmonic  analyzer  (SHA)  (Ref.  5-3)  at  1.06  ym.  Powders  of 
random  size,  212  ym,  ground  from  single  crystals  in  an  agate 
mortar  and  pestle  were  examined  as  monolayer  samples  between  two 
coverglasses . A chromel-alumel  thermocouple  was  placed  near  the 
sample  as  a means  of  monitoring  the  temperature.  Temperature 
fluctuations  were  less  than  + 3°C. 

The  results  of  the  SHG  measurements  on  43m  - mm2  phase  transi- 
tion on  Ni-Br  boracite  are  shown  in  Figure  58.  The  normalized 
second  harmonic  intensity  is  on  the  order  of  0.1  a-SiC^  at  room 
temperature.  The  phase  transition  is  observed  to  occur  at  400  + 

4 K which  is  in  good  agreement  with  the  specific  heat  determina- 
tion of  398.4  + 0.4  K.  No  thermal  hysteresis  was  measured  due 
to  temperature  fluctuations;  however,  the  discontinuous  jump  in 
the  second  harmonic  intensity,  I^W,  supports  the  claim  of  a 
first-order  transition.  Both  phases  are  seen  to  be  noncentro- 
symmetric  in  accordance  with  crystal  point  group  assignments  in 
other  boracites.  The  lack  of  structural  data  on  Ni-Br  boracite 
precludes  a definitive  assignment  of  the  crystal  phases. 
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LITERATURE  SURVEY 


6.1  Basic  References 

There  are  three  recent  review  papers  on  the  subject  of  improper 
f erroelectr ics . The  most  comprehensive  review  of  the  theory 
appears  in: 

A.  P.  Levanyuk  and  D.  G.  Sannikov  1974 

Improper  Fer roelectr ics 

Soviet  Physics:  Uspekhi  199-214. 

Particular  attention  is  given  to  contributions  made  in  the 
Soviet  literature,  some  of  which  are  not  available  in  English 
translation.  The  authors  performed  the  first  theoretical 
analysis  of  the  dielectric  properties  of  improper  fer roelectr ics 
in  1968  and  continue  to  actively  contribute  to  the  field. 
Significantly  more  attention  is  paid  to  experimental  observa- 
tions in: 

V.  Dvorak  1974 
Improper  Ferrcelectr ics 
Ferroelectric  1_,  1-9 

This  author  is  responsible  for  suggesting  the  name  "improper 
fer roelectr ics"  for  this  class  of  materials.  The  name  seems  to 
have  gained  acceptance,  although  papers  continue  to  appear 
occassionally  referring  to  a weak,  atypical,  complex,  faint, 
partial,  or  extrinsic  fer roelectr ics . The  most  recent  review  is 
presented  by: 

Y.  Ishibashi  and  Y.  Takagi  1976 
Improper  Ferroelectric  Phase  Transitions 
Japan.  J.  Appl.  Phys . 15 , 1621-1636. 

In  this  review,  attention  is  focused  on  a number  of  unresolved 
questions,  both  experimental  and  theoretical. 

Much  of  the  theory  of  improper  f erroelctr ics  rests  on  Landau's 
theory  of  continuous  phase  transitions.  The  original  papers, 

L.  Landau  1937A 

Zur  Theorie  der  Phasenumwandlung  I 
Phys.  Z.  Sowjet  11 , 26. 
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and 


L.  Landau  1937B 

Zur  Theocie  der  Phasenumwandlung  II 
Phys.  Z.  Sowjet  11 , 545 

here  appeared  in  English  translation  in: 

L.  D.  Landau,  Collected  Works 
D.  ter  Haar,  Editor 
Pergamon,  1965. 

These  are  still  the  best  expositon  of  this  theory.  The  first 
theoretical  suggestion  that  an  improper  ferroelectric  might 
occur  appears  in, 

V.  L.  Indenbom  1960 

Phase  Transitions  Without  Change  in  the  Number  of 
Atoms  in  the  Unit  Cell  of  the  Crystal 
Soviet  Physics:  Crystallograph  _5,  106-115. 

This  paper  basically  analyzes  the  crystallographic  point  groups 
to  determine  all  possible  phase  transitions  consistent  with  the 
Landau  theory  in  which  the  unit  cell  does  not  change  during  the 
transition.  An  even  more  detailed  analysis  of  precisely  the 
same  probem  is  given  by. 


V.  Janovec,  D.  Dvorak  and  J.  Petzelt  1975 
Symmetry  Classification  and  Properties  of 
Equi-Translational  Structural  Phase  Transitions 
Czech.  J.  Phys.  B25 , 1362-1396. 

This  paper  tabulates  a wealth  of  information.  An  extension  of 
these  two  analyses  to  include  transitions  in  which  the  unit  cell 
changes  has  been  performed  for  the  special  case  of  improper 
fer roelectr ics  which  however  are  not  simultaneously  ferro- 
elastic.  This  partial  extension  is  reported  in: 

P.  Toledano  and  J.  C.  Toledano  1976 
Order-Parameter  Symmetries  for  Improper 
Ferroelectric  Nonferr oelastic  Transitions 
Phys.  Rev.  B14,  3097-3109. 
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6 . 2 Experimental  Literature  On  Boracites 

' » The  boracite  family  of  compounds  (X3B7013Y,  where  X 

is  a divalent -metal  and  Y a halogen)  is  known  to  contain  a large 
number  of  ferroelectric  materials,  some  of  which  have  been 
identified  as  improper  f erroelectr ics . The  large  number  of 

* compounds  and  the  many  interesting  properties  of  these  materials 
has  led  to  a large  experimental  literature.  In  Table  22  this 
literature  is  classified  in  terms  of  physical  properties  and  the 
particular  material.  In  Table  23  the  properties  of  the  observ- 

* ed  ferroelectric  transitions  are  summarized.  All  references  are 
to  papers  listed  in  Appendix  B. 

'•I,'  ' ' 
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TABLE  22.  EXPERIMENTAL  LITERATURE  ON  BORACITES 


-J] 


Dielectric  Constant  ^ 

Mg-Cl:  Jona  1959;  Sonin  & Zheludev  1963;  Bochkov, 

Gavrilova  et  al . 1975;  Albers  et  al . 1976. 

Cr-Cl:  Bochkov  & Drozhdin  1975 

J 

Fe-Cl:  Bochkov  et  al.  1973 

Fe-Br:  Bochkov  et  al . 1973 

Fe-I  : Schmid  et  al . 1976 

Co-Cl:  Bochkov  et  al . 1973 

Co-Br:  Ascher  et  al . 1964 

Co-I  : Ascher  et  al . 1964;  Smutny  & Fousek  1970; 

Smutny  et  al.  1974  B;  Smutny  & Konak  1975; 

Drozhdin  et  al  1975. 

Ni-Cl:  Jona  1959;  Ascher  et  al . 1964. 

Ni-Br:  Ascher  et  al . 1964 

Ni-I  ; Ascher  et  al . 1966;  Smutny  et  al.  If74  A; 

Smutny  et  al.  1974  B. 

Cu-Br ; Drozhdin  et  al.  1975. 

Zn-Cl:  Jona  1959;  Bochknov  et  al . 1973. 

Cd-Cl : Jona  1959. 

Spontaneous  Polarization 


Mg-Cl: 

Bochkov,  Gavrilova  et  al . 1975;  Fridkin 
1974;  Torre  et  al . 1972. 

et  al . 

Cr-Cl: 

Bochkov  & Drozhdin  1975;  Fridkin  et  al . 

1974. 

Fe-Br : 

Fridkin  et  al . 1974 

Fe-I  : 

Schmid  et  al . 1976 

Co-I  : 

Smutny  & Fousek  1970; 
Drozhdin  et  al . 1975. 

Smutny  & Konak  1975; 

Ni-I  : 

Miyashita  & Murakami 

1970 

Cu-Br : 

Drozhdin  et  al . 1975 
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TABLE  22.  EXPERIMENTAL  LITERATURE  ON  BORACITES  (CONT'D) 

* Crystal  Structure 

Review:  Nelmes  1974  - a review  of  the  properties  of 

boracites  with  an  emphasis  on  the 
structual  studies  of  boracites. 

X-Ray  Studies: 

Ito  et  al . 1951  - structure  of  Mg-Cl  (43m-mm2). 

Kobayashi  and  Mizutani  1970  - lattice  strains  of  Fe-I 
strains  and  parameters  of  Fe-I. 

Schmid  1970  - trigonal  boracites. 

Dowty  and  Clark  1972  - atomic  displacements. 

Kobayashi  et  al.  1972  - Fe-I  at  low  temperatures. 

Torre  et  al . 1972  - structural  basis  of  fer roelectr icity 

Felix  et  al . 1974  - x-ray  diffuse  scattering  for  boracite 

Nelmes  and  Thornley  1974  - Cr-Cl,  full  3-diam.  structure 

Thornley  and  Nelmes  1974  - Cr-Cl 

Nelmes  and  Thornley  1976  - Ni-I 

* Thornley  et  al . 1974  - Cu-Cl 

Mossbauer  Investigations: 

Schmid  and  Trooster  1967  - Fe-Cl,  Fe-Br  and  Fe-I. 
f Trooster  1969 

EPR  Spectra: 

Rivera  et  al  1976A  - Mg-Cl(Mn2+) 

, Rivera  et  al  1976B  - Zn-Cl(Mn2+) 

Neutron  Scattering: 

Thornley  et  al.  1976 

Domain  Studies 


Ascher  et  al . 1964  - Ni-Cl  domain  wall  motion 
Kobayashi  et  al.  1968  - Fe-I 
Schmid  1969  - Ni-Cl,  Ni-Br 
Bochkov,  Gavrilova  et  al . 1975  - Mg-Cl 
Drozhdin  et  al.  1975  - Cu-Br 


163 


TABLE  22.  EXPERIMENTAL  LITERATURE  ON  BORACITES  (CONT’D) 


Thermal  Measurements 


The  Latent  Heat  of  the  Transition: 

Bochkov  et  al . 1973  - Fe-Cl,  Fe-Br,  Co-Cl,  Zn-Cl. 
Bochkov,  Gavriola  et  al . 1975  - Mg-Cl 

Specific  Heat: 

Smutny  1972  - Co-I 
Schmid  et  al . 1976  - Fe-I 

Pyroelectricity 

Bochkov  and  Drozhdin  1975  - Cr-Cl 
Torre  et  al.  1972  - Mg-Cl 

Optical  and  Electroopt ical  Properties 

Schmid  and  Trooster  1967  - optical  study  of  phase  transitions 

Kobayashi,  Mizutani,  Hara  et  al.  1970  - optical  refractive 
index  of  Fe-I 

Pisarev  et  al.  1970  - crystal  field  theory  and  optical 
absorpt ion 

Fridkin  et  al . 1974  - temperature  dependence  of  the  optical 
absorption  edge 

Smutny  and  Konak  1975  - Co-I 

Raman  Spectroscopy 

Lockwood  1976  A - Cr-Cl 
Lockwood  1976  3 - Cr-Cl,  Mn-Cl 

Piezoelectricity 


Mg-Cl: 

Sonin  and 
Albers  et 

Zheludev  1963;  Torre  et  al.  1972; 
al.  1976 

Fe-I: 

Kobayashi 
Kobayashi , 

and  Mizutani  1970; 

Mizutani,  Hara  et  al . 1970 

Co-I: 

Smutny  and 

Albers  1972 
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T^BLE  22.  EXPERIMENTAL  LITERATURE  ON  BORACITES  (CONT'D) 


DC  Resistivity 

Sonin  and  Zheludev  1963  - Mg-Cl 
✓ 

Smutny  et  al.  1974  A - Ni-I 
Schmid  et  al . 1976  - Fe-I 

Crystal  Growth 

Jona  1959;  Schmid  1965;  Schmid  1969; 

Nassau  and  Shiever  1972 

Ferromagnetoelectricity  of  Ni-I  Boracite 

Ascher  et  al.  1966 
Miyashita  and  Murakami  1970 

The  Effect  of  Hydrostatic  Pressure  on  the  Ferroelectric  Phase 
Transition  Temperature 

Fousek  et  al . 1972  - Co-I 
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TABLE  23. 

OBSERVED  FERROELECTRIC  TRANSITIONS  IN  BORACITES 

Boracite 

Phase 

Transition 

Nature  of 

Ref. 

Transition 

Temperatures(  C) 

Transition 

Mg- Cl 

43m-mm2 

265 

1 

2 

Cr-Cl 

T3m-mm2 

-9 

1 

2 

Fe-Cl 

T3m-mm2 

320 

1 

1 

-m 

225 

1 

1 

-3m 

203 

1 

1 

Fe-Br 

4 3m- mm2 

222 

1 

1 

-3m 

130 

1 

1 

Fe-I 

43m-mm2 

72 

-55  - -70 

6,7 

Co-Cl 

43m-mra2 

350 

1 

1 

-m 

265 

2 

2 

-3m 

195 

2 

2 

Co-3r 

43m-mm2 

185 

Co-I 

T3m-rara2 

-76 

1 

3,4, 

-235 

3 

Ni-Cl 

T3m-mn2 

337 

1 

V 

Ni-Br 

43m-mm2 

125 

Ni-I 

4 3m-mm2 

-209 

Gu-Gl 

?3m-mm2 

92 

Cu-3r 

4 3m-ram2 

-53 

1 

4 

-248 

Zn-  01 

43m-mm2 

430 

1 

1 

-m 

291 

1 

2 

-3m 

199 

1 

2 

Cd-Cl 

T3m-mm2 

525 

Mg-Br 

21 

3 

Mn-Cl 

407 

3 

Mri-Br 

291-294 

7 

J 

Mn-I 

139 

3 

Zn-Br 

312 

3 

Zn-I 

413 

3 

Cd-Br 

462 

3 

Cd-I 

342 

3 

References 

: (1)  Bochkov  et  al  1973  (2) 

Friakin 

et  al  1974 

(3)  Ascher 

et  al  1964  (4) 

Droshdin 

et  al  1975 

(5)  Fousek 

et  al  1972  (6) 

Kobayashi  et  al  1972 

(7)  Schmid 

& Troo3ter  1367 
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6.3 


Experimental  Literature  on  Molybdates 


Most  members  of  the  molybdate  family  [chemical  composition: 
I^MoO^)^  where  R is  a rare  earth]  have  been  shown  to  be 
ferroelectric.  Among  these,  gadolinium  molybdate  and  terbium 
molybdate  have  been  clearly  established  as  improper  ferroelec- 
trics;  it  is  expected  that  the  other  molybdates  will  also  be 
identified  as  improper  f erroelectrics  when  studied  in  detail. 

In  addition  to  f erroelectr icity , the  molybdates  also  exhibit 
f erroelastici ty . The  resulting  nonlinear  thermal,  electrical 
and  mechanical  properties  of  the  molybdates  have  provided  a 
fruitful  field  for  experimental  study.  This  rather  extensive 
literature  has  been  classified  in  Table  24.  Papers  cited  in  the 
table  are  listed  in  Appendix  B. 


' 

TABLE  24.  MOLYBDATES  EXPERIMENTS 

BHILLOUII!  SCATTERING  (GKO) 

Itoh  ft  Nakamura  1973  ; Itoh  ft  Nakamura  197-1  ; 
luspir.  ft  Hauret  1976  AftB 

J •' LOSTMSTRIC  MB-.SIJ KKMSKTS 

GMC  : Nassau  et  '-.1  1965  ; Fouskova  1969  ; 

Cheung  ft  Ullman  1974  ; Me^-urni  et  -il  1974 
TMC  : Keve  et  al  1570 

J 

1EYSTAL  r.lQ'-rv 

Borchardt  ft  Bisrstedt  1967  ; Drobyshev  et  al  196?  ; 

Keve  et  al  1970  ; Bohm  ft  Kursten  1 97 1 

ORYST'.I  S?:;U  GTIJRE 


X-Ray  Studies  : Bernstein  1965  ; Nassau  et  el  1966  ; 

Aizu  et  al  1969  ; Drobyshev  et  al  1965  ; 

Jamieson  et  al  1969  ; Nassau  et  al  196?  ; 

Newnham  et  al  1969  J Kvapil  ft  John  13nC  ; 

Prewitt  1970  ; Keve  et  al  1971  ; Brixner  et  al  1972  ; 
Jeitschko  197?  ; Kobayasbi  et  al  :972  ; 

Me^umi  et  al  1974 

Neutron  Scattering  : Axe  et  al  197'.  ; Dorner  et  al  1979 
BIELECTRI 0 CONSTANTS 

GMO  : Borchardt  ft  3ierstedt  1 967  ; Gross  et  al  1968  ; 

Aizu  et  al  1969  ; Smith  ft  Burns  1969; 

Fousek  ft  Konak  1972B 

TKC  : Borchardt  ft  Bierstedt  1967  ; Keve  et  al  1970 

DOUBLE  EYSTE-tESIS  LOO? 

Kumada  1973  ; Kumada  1974  ; 

Sawaguchi  ft  Gross  1972 


i 


TABLE  24.  MOLYBDATES  EXPERIMENTS  (CONT'D) 


DOMAIN  STRUCTURES  A ANTI-PHASE  BOUNDARIES 

Borchardt  A Bierstedt  1966  ; Cummins  1970  ; Kumada  1970  ; 
Kumada  et  al  1970  ; Bohm  A Kursten  1971  ; Barkley  et  al  1972  ; 
Kumada  1972  ; Kursten  A Bohm  1972  ; Shepherd  A Barkley  1972  ; 
Suzuki  1972  ; Barkley  A Jeitschko  1973  ; 

Bohm  A Kursten  1973  ; Flippen  A Haas  1975  ; 

Kursten  et  al  1973  ; Meleshina  et  al  1974  ; Flippen  1975 

EFFECT  OF  HYDROSTATIC  PRESSURE  ON  PHASE  TRANSITION  IN  GMO 
Shirokov  et  al  1972 

ELASTIC  PROPERTIES 

GMO  : Cross  et  al  1969  ; Epstein  et  al  1970  ; 

Chizhikov  et  al  1971  ; Hochli  1972  ; 

Courdille  A Dumas  1973  ; Scheiding  1973  ; 

Courdille  A Dumas  1974  ; Esayan  et  al  1974  ; 

Esayan  et  al  1975 
TMO  : Dorner  et  al  1972 

ELECTROOPTIC  PROPERTIES  (GMO) 

Smith  A Burns  1969  ; Nakamura  et  al  1971  ; 

Fousek  A Konak  1972  AAB  ; Kobayashi  & Ascher  1972  ; 

Kumada  1972  ; Kobayashi  A Miyazaki  1975  ; 

Anistratov  et  al  1976  ; Kojima  et  al  1976  ; 

Theophanous  1976 

OPTICAL  PROPERTIES 

Refractive  Indices  A Birefringence  : 

(GMO)Aizu  et  al  1969  ; Smith  A Burns  1969  ; Cummins  1970  ; 
Nakamura  et  al  1971  ; Fousek  A Konak  1 9 7 2 A ; 

Kobayashi  A Ascher  1972  ; Shaldin  et  al  1973  ; 

Fousek  A Glogarova  1975  ; Anistratov  et  al  1976  ; 

Kojima  et  al  1976 

(TMO)Keve  et  al  1970  ; Shaldin  et  al  1973 
( Nd?(MoO^  ) ^Nassau  et  al  1965 
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TABLE  24.  MOLYBDATES  EXPERIMENTS  (CONT'D) 

Transmission  Spectra  : Nassau  et  al  1965  ; Petzelt  1971  ; 

Petzelt  & Dvorak  1971  ; Zeidler  ft  Ullman  1973 
Ferroelectric  Laser  Host  : Borchardt  1963  ; 

Borchardt  ft  Bierstedt  1966 

Optical  Device  Application  : Kumada  1970  ; Kumada  1972 
Second  Harmonic  Generation  : Miller  et  al  1971 

PHOTOELASTICITY  ( GMO) 

Esayan  et  al  1974  ; Esayan  et  al  1975  ; 

Fousek  ft  Glogarova  1975  ; Martynov  et  al  1975  ; 

Anistratov  et  al  1976 

PIEZOELECTRICITY 

Aizu  et  al  1969  ; Cummins  1970  ; Kobayashi  et  al  1972  ; 
Scheiding  ft  Schmidt  1972  ; Scheiding  1973  ; 

Agishev  et  al  1976 

PYROELECTRICITY 

GMO  : Rabinovich  1969  ; Ullman  et  al  1971  ; 

Ullman  et  al  1972  ; Ullman  et  al  1975 
TMO  : Keve  et  al  1970 

RAMAN  SPECTROSCOPY  (GMO) 

Fleury  1970  ; Shepherd  1971  ; Shepherd  & Barkley  1972  ; 

Ullman  et  al  1973  ; Ganguly  et  al  1975  A&B  ; 

Ganguly  et  al  1976  ; Shigenari  et  al  1976  ; 

Takagi  et  al  1976 

SPONTANEOUS  POLARIZATION 

GMO  : Borchardt  ft  Bierstedt  1966  ; Borchardt  ft  Bierstedt  1967 
Aizu  et  al  1969  ; Cummins  1970  ; 

Kumada  et  al  1970  ; Sawaguchi  & Cross  1973 
TMO  : Borchardt  ft  Bierstedt  1967  ; Keve  et  al  1970 

SURFACE  ACOUSTIC  WAVE  VELOCITY  ON  GMO  : Toda  et  al  1973 

Na?e(Mo04)  - A NEW  FEKROELASTIC  : Dudnik  ft  Stolpakova  1976 
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6.4 


Experimental  Literature  on  Ferroelectric  Sulfates 


Tables  25  and  26  summarize  the  results  of  our  search  of  the 
experimental  literature  on  ammonium  sulfate,  (NH4)2  SC>4, 
and  related  materials.  The  following  classes  of  materials  are 
included : 


(1) 

(NH 

4)2S04  analogs: 

(NH4)2X,  with 

X 

= so4 

' BefV 

Se04 

(2) 

Bisulfates:  NH.  HSC. 

4 4 

, RbHS04,  (NH4 

>3 

h(so4)2 

(3) 

ABX 

compounds  with  A 

= Li , Na ; B = 

= Na,  K, 

Rb,  NH4, 

i Cs , 

N2H5 

; X = BeF 4 , S04, 

Se04 

(4) 

Langbeinites:  (Mi+)2 

(M2+)2  (S04) 3 

with 

Mi+ 

= K,  Rb,  Tl,  NH4; 

M2+  = Zn,  Co, 

Mn 

f Mg, 

Cd 

ong 

the 

langbeinites. 

only  ( N H 4 ) 

2 

Cd  2 

(so4)3 

and 

T £ 2 Cd2  ( S04 ) ^ have  been  identified  as  ferroelectric. 
Table  25  lists  all  known  ferroelectrics  among  these  classes, 
gives  the  space  group  (or,  in  parentheses,  the  point  group)  of 
each  phase  if  known,  and  the  temperatures  of  all  reported  phase 
transitions.  Ferroelectric  phases  have  been  characterized  as 
pseudoproper  when  a Curie-Weiss  law  is  obeyed  in  a small  region 
above  the  transition  temperature  and  as  improper  where  there  is 
no  Curie-Weiss  behavior  and  only  a small  dielectric  anomaly.  In 
Table  26  the  literature  is  classified  according  to  the  material 
and  physical  property  studied,  all  references  are  to  papers 
listed  in  Appendix  B. 
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TABLE  25.  THE  FERROELECTRIC  SULFATE  FAMILY 


CTS  Dlcadmlum  D1 thallium 

Sulfate 


TABLE  26.  THE  FERROELECTRIC  SULFATE  FAMILY. 


t LITERATURE  GUI BE 

(NH4)2S04  ft  (ND4)2S04 

CALORIMETRIC  MEASUREMENTS 
I Hoehino  et  al  1958 

CRYSTAL  STRUCTURE 

Neutron  Diffraction  : Hamilton  1969  ; Schlemper  ft  Hamilton  1966 
NMR  : Blinc  et  al  1972  ; Kydon  et  al  1967  ; Miller  et  al  1962  ; 

. O'Reilly  ft  Tsang  1967  ; O'Reilly  ft  Teang  1969 

X-Ray  Scattering  : Okaya  et  al  1958 

DIELECTRIC  CONSTANTS 

Anistratov  ft  Martynov  1970  ; Hoshino  et  al  1958  ; 

/ Ikeda  et  al  1973  ; Kamiyoshi  1957  ; Ohshima  ft  Nakamura  1966  ; 

Sawada  et  al  1975;  Yoshihara  et  al  1976 

ELASTIC  PROPERTIES 

Hoshino  et  al  1958  (Thermal  Expansion)  ; 

Ikeda  et  al  1973  (Elastic  Compliance  Coefficient, Electrostrictive 
Constant)  ; Makita  et  al  1976  (Mechanical  Twins)  ; 

Nakanishi  et  al  1973  (Elastic  Softening)  ; 

Sawada  et  al  1976  (Mechanical  Twins) 

ELECTROOPTIC  PROPERTIES 

Anistratov  ft  Martynov  1970  ; Pousek  ft  Konak  1976 

INFRARED  SPECTROSCOPY 

Belousov  et  al  1975  ; Jain  et  al  1973  ; 

KumaraAcharya  ft  Narayanan  1973  ; Petzelt  et  al  1974  ; 

Schutte  ft  Heyns  1970  ; Torrie  et  al  1972 

OPTICAL  PROPERTIES 

Anistratov  1972  ; Anistratov  ft  Martynov  1970  ; 

Fousek  ft  Konak  1976  ; Strukov  1962 

PYROELECTRICITY 

Sawada  et  al  1975 

RAMAN  SPECTROSCOPY 

KumaraAcharya  ft  Narayanan  1973  ; Torrie  et  al  1972  ; 

Unruh  ft  Ayere  1976 


TABLE  26.  THE  FERROELECTRIC  SULFATE  FAMILY  (CONT'D) 


1 


J 

SPONTANEOUS  POLARIZATION 

Hoshino  et  al  1958  ; Ikeda  et  al  1973  ; 

Matthias  ft  Remeika  1956  ; Sawada  et  al  1975  ; 

Unruh  1970  ; Unruh  ft  Ayere  1976 

THEORIES 

Ikeda  et  al  1973  ; Jain  ft  Bist  1974  ; Sawada  et  al  1973 


f 


(NH4)2BeP4  ft  (NB4)2BeP4 

CALORIMETRIC  MEASUREMENTS 

Hoshino  et  al  1958  ; Strukov  et  al  1973 

CRYSTAL  STRUCTURE 

NMR  : Aleksandrova  ft  Shcherbakov  1970  ; Burns  1961  ; 
Miller  et  al  1962 

X-Ray  ; Makita'ft  Yamauchi  1974  ; Okaya  et  al  1958 
DIELECTRIC  CONSTANTS 

Hoshino  et  al  1958  ; Ohshima  ft  Nakamura  1966  ; 
Pepinsky  ft  Jona  1957  ; Strukov  et  al  1962  ; 

Strukov  et  al  1973 

OPTICAL  ft  ELECTROOPTIC  PROPERTIES 

Anistratov  1972  ; Konak  ft  Matras  1976 

SPONTANEOUS  POLARIZATION 

Pepinsky  ft  Jona  1957  ; Strukov  et  al  1962 

THEORIES  : Dvorak  et  al  1975 

(NH4)2Se04  ft  (ND4)2Se04 

INFRARED  ft  RAMAN  SPECTROSCOPY 

Sankaranarayanan  ft  Krishnan  1972 


! t 


L 


Alexandrov  et  al  1976  ; Flerov  et  al  1975  ; 

Pepinsky  et  al  1958 

CRYSTAL  STRUCTURE 

Miller  et  al  1962  (NMR)  ; Nelraes  1972  (X-Ray)  ; 

Watton  et  al  1976  (NMR) 

DIELECTRIC  CONSTANTS 

Alexandrov  et  al  1976  ; Flerov  et  al  1975  ; Pepinaky  et  al  1958  ; 
Strukov  et  al  1962  ; Taraaki  et  al  1976 


DOMAIN  STUDIES 


Kondo  et  al  1974 


EFFECT  OF  HYDROSTATIC  PRESSURE  ON  PHASE  TRANSITION 
Polandov  et  al  1969 

OPTICAL  PROPERTIES  : Alexandrov  et  al  1976 

SPONTANEOUS  POLARIZATION 

Alexandrov  et  al  1976  ; Flerov  et  al  1975  ; 
Pepinsky  et  al  1958  ; Strukov  et  al  1962 

THEORIES  : Aizu  1974  ; Kondo  et  al  1975 


RbHSO . 4 RbDSO , 


DIELECTRIC  CONSTANT  4 SPONTANEOUS  POLARIZATION 

Alexandrov  et  al  1976  ; Pepinsky  4 Vedam  I960 

SPECIFIC  HEAT  4 OPTICAL  PROPERTIES 
Alexandrov  et  al  1976 

EFFECT  OF  HYDROSTATIC  PRESSURE  ON  THE  PHASE  TRANSITIONS 
Gesi  4 Ozawa  1975 
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TABLE  26.  THE  FERROELECTRIC  SULFATE  FAMILY  (CONT'D) 


(WHi)3H(S04)2 

CRYSTAL  STRUCTURE  : Watton  et  al  1976  (NMR) 

DIELECTRIC  CONSTANTS  : Geai  1976  ; Pepinsky  et  al  1958 
SPECIFIC  HEAT  : Geai  1976 

LiNH4S04  ft  LiND4S04 
CRYSTAL  STRUCTURE 

Aleksandrov  et  al  1975  ; Dollaae  1969  (X-Ray)  ; 

Yuzvak  et  al  1975  (NMR) 

DIELECTRIC  CONSTANTS  ft  SPONTANEOUS  POLARIZATION 
Mitsui  et  al  1975  ; Yuzvak  et  al  1975 

DOMAIN  STUDIES  : Hildrnann  et  al  1975 

ELASTIC  CONSTANTS  : Aleksandrov  et  al  1975 

INFRARED  ft  RAMAN  SPECTROSCOPY  : KuraaraA.cha.rya  * Narayanan  1973 
OPTICAL  ft  ELECTROOPTIC  PROPERTIES  : Anistratov  ft  Mel'nikova  1975 

Li(N2H5)S04 

CALORIMETRIC  MEASUREMENTS  : Pepinsky  et  al  1958 
CONDUCTIVITY  : Howell  et  al  1967  ; Schmidt  et  al  1971 
CRYSTAL  STRUCTURE 

Neutron  Scattering  : Padmanabhan  ft  Balasubramanian  1967 
X-Ray  : Brown  1964  ; Van  den  Hende  ft  Boutin  1964 
NMR  : Cuthbert  ft  Petch  1963  ; Knispel  ft  Petch  1971 

DIELECTRIC  CONSTANTS 

Pepinsky  et  al  1958  ; Schmidt  et  al  1971 


TABLE  26.  THE  FERROELECTRIC  SULFATE  FAMILY  (CONT'D) 


NaNH4Se04,2H20  ft  NaNH4S04> 2H20 
CRYSTAL  STRUCTURE 

NMR  : Aleksandrova  et  al  1972  ; Shirokov  et  al  1975 
DIELECTRIC  CONSTANTS 

Aleksandrov  et  al  1971  ; Makita  & Sekido  1965  ; 

Shirokov  et  al  1975  ; Zaitseva  et  al  1975 

ELASTIC  ft  PIEZOELECTRIC  PROPERTIES  : Zaitseva  et  al  1975 

INFRARED  ft  RAMAN  SPECTROSCOPY  : KumaraAcharya  ft  Narayanan  1973 

OPTICAL  ft  ELECTROOPTIC  PROPERTIES  : Anistratov  * Mel'nikova  1975 

PYROELECTRICITY  : Makita  ft  Sekido  1965 

MeLiBeF4  , MeLiS04  , MeLiSe04 

Me  - K,  Rb,  Tl,  NH4,  Cs  (except  NH4LiS04) 

CRYSTAL  STRUCTURE 

Chung  ft  Hahn  1972  ; Hahn  et  al  1969  ; 

LeRoy  ft  Aleonard  1970  ; LeRoy  ft  Aleonard  1972 

DIELECTRIC  CONSTANTS  : Shiroishi  et  al  1976  (RbLiS04) 

OPTICAL  ft  ELECTROOPTIC  PROPERTIES 

Anistratov  ft  Mel'nikova  1975  (LiNH4Se04) 

INFRARED  ft  RAMAN  SPECTROSCOPY 
Rao  et  al  1976  (LiKS04)  ; 

Sankaranarayanan  ft  Krishnan  1972  (LiNH4Se04) 


1 
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TABLE  26.  THE  FERROELECTRIC  SULFATE  FAMILY  (CONT'D) 


* 

1 


(NH4)2Cd2(S04):,  4 Tl2Cd2(S04), 

CRYSTAL  STRUCTURE  : Ikida  4 Yasuda  1975  (CTS,  X-Ray) 
DIELECTRIC  CONSTANTS 

CAS  : Glogarova  4 Fousek  1973  ; Glogarova  et  al  1972  ; 

Jona  4 Pepinsky  1956  ; Ohshlma  4 Nakamura  1966 
CTS  : Brezina  4 Glogarova  1972  ; Ikida  4 Yasuda  1975 

DOMAIN  STUDIES 

CAS  : Glogarova  4 Fousek  1973  ; Glogarova  et  al  1972 
CTS  : Brezina  4 Glogarova  1972 

ELASTIC  PROPERTIES 

CAS  : Glogarova  et  al  1972 
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This  Appendix  contains  intuitive  physical  overviews  of  the 
theoretical  concepts  used  in  the  study  of  improper  f erroelectr ics . 
There  are  two  languages  in  wide-spread  use.  One  is  provided  by 
Landau's  macroscopic  theory  of  continuous  phase  transitions  and 
the  other  originates  with  the  microscopic  soft-mode  analysis  of 
structural  phase  transitions.  Using  these  ideas  we  distinguish 
between  three  different  types  of  extrinsic  (or  improper  in  the 
larger  sense)  f erroelectrics . Since  all  known  improper  ferro- 
electrics  also  exhibit  f erroelasticity , we  also  summarize  the 
relationship  between  the  full  free  energy  and  the  reduced  free 
energy  which  is  applicable  to  stress-free  measurements. 

1 • 1 Landau  Theory  of  Continuous  Phase  Transitions 

Landau's  theory  is  basically  a macroscopic  thermodynamic 
analysis  appplicable  to  any  continuous,  so-called  second-order 
phase  transition.  In  a second-order  phase  transition,  the 
transition  temperature  marks  the  point  where  the  system's 
properties  begin  to  change;  there  is  no  sudden  jump  in  the 
system's  properties  at  the  transition  temperature,  as  occurs  in 
the  so-called  first-order  transitions.*  While  the  Landau  analy- 
sis is  strictly  applicable  only  to  such  gradual  or  continuous 
phase  transitions,  it  can  be  extended  to  the  first-order  tran- 
sitions in  which  the  discontinuities  exhibited  at  the  transition 
temperature  are  small. 

1.1.1  Free  Energy  Considerations 

The  first  ingredient  of  the  Landau  theory  is  to  exploit  the 
hypothesized  continuity  of  the  phase  transition  by  expanding  the 
relevant  thermodynamic  free  energy  $ as  a power  series  in  the 
so-called  order  parameter,  n . "Order  parameter"  is  the  generic 
name  for  the  property  which  distinguishes  the  two  phases  of  the 
system.  For  example,  in  an  ordinary  ferroelectric,  the  order 
parameter  is  just  the  electric  polarization.  The  order  param- 
eter is  chosen  so  that  it  vanishes  in  one  phase  (the  high 
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symmetry  or  prototypic  phase)  and  gradually  begins  to  appear  at 
the  transition  temperature,  thus  making  the  transition  to  the 
new  phase  (the  low  symmetry  or  ferroic  phase).  This  expansion 
of  the  free  energy  *,  which  is  assemed  to  be  valid  both  above 
and  below  the  transition  temperature,  has,  in  general,  the 
form: 

*(n,p,T)  = «Q(p,T)  + n^fpjT)  + n2*2(P'T)  + 

+ n,*3(p,T)  + n4$4(p,T)  (!) 

where  p denotes  the  pressure  and  perhaps  other  thermodynamic 
properties  which  are  not  involved  in  the  transition. 

This  free  energy,  of  course,  embodies  a complete  description  of 
the  thermodynamic  behavior  of  the  system;  its  consequences  are 
deduced  in  a number  of  ways.  For  example:  the  order  parameter 
n takes  the  value  that  gives  a minimum  for  «;  the  inverse 
susceptibility  ( in  the  ferroelectric  case,  the  dielectric 

stiffness)  is  just  the  curvature  of  * at  that  minimum  = 
32*/3n2;  the  entropy  S = -34/3T,  and  heat  capacity  C = -T32»/3T2 
are  contained  in  the  temperature  dependence  of  the  free  energy. 

Landau's  analysis  reveals  that  if  the  free  energy  $ , given  in 
Equation  (1),  is,  in  fact,  to  describe  a second-order  tran- 
sition additional  requirements  must  be  met.  First,  the  odd 
degree  terms  must  vanish;  second,  the  coefficient  of  the  second 
degree  term  must  change  sign  at  the  transition  temperature;  and 
third,  the  coefficient  of  the  fourth  degree  term  must  be  pos- 
itive. We  shall  see  below  that  these  requirements  amount  to 
n bstantive  restrictions  on  the  nature  of  the  order  parameter 
• • d -.he  nature  of  the  high  symmetry  and  low  symmetry  phases  of 
* -./stem.  At  this  point,  let  us  just  review  the  general 

• ---  jf  the  phasv.  transition,  assuming  that  these  require- 

--•••  * • * ‘jifilled. 
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The  free  energy  for  a second  order  transition  then,  has  the 
form : 


*(n,T)  = + a ( T-Tc ) n 2 + Bn" 


(2) 


t where  B > 0 , Tc  is  the  transition  temperature  and  a > 0 if  the 

high  symmetry  phase  is  also  the  high  temperature  phase.  (This 
is  usually,  but  now  always,  the  case).  In  Figure  1,  * is 
plotted  as  a function  of  the  order  parameter  n for  a number  of 
different  temperatures.  Figure  2 depicts  the  temperature 
dependence  of  the  order  parameter  and  inverse  susceptibility 
deduced  from  this  free  energy.  (In  the  case  of  a proper  ferro- 
electric, these  correspond  to  the  electric  polarization  and  the 
dielectric  stiffness,  respectively.)  From  Figure  1 we  see 
that  for  T > T , the  minimum  of  $ is  at  n = 0 and  x the 
curvature  of  * at  the  minimum  is  quite  large.  As  the  temper- 
ature is  reduced,  the  minimum  remains  at  n = 0 and  the  curvature 

_ i 

x vanishes.  As  the  temperature  is  reduced  below  Tc , the 
central  minimum  bifurcates  into  two  symmetric  minima  which 
gradually  move  further  from  the  origin,  with  their  curvature 
increasing  concomitantly. 


Figure  I:  Free  energy  *(n,T)  vs.  order  parameter 

n for  various  temperatures  above  and  below  critical 
temperature  T where  a second-order  phase 
transition  ccurs. 
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Figure  2:  (a)  Temperature  dependence  of  the 

magnitude  of  the  order  parameter  n for  a system 
which  undergoes  a second-order  phase  transition  at 
T . The  graph  gives  the  value  of  n at  the 
absolute  minimum  of  the  free  energy  which  is  shown 
i-ni^g.  1.  Below  T , n varies  as  (T 
T)  . (b)  Temperature  dependence  of  the  inverse 

susceptibility  x , for  a system  which  undergoes 
a second-order  phase  transition  at  T . This 
graph  gives  the  value  of  the  curvature  of  the  free 
energy  (shown  in  Fig.  1)  at  its  absolute  minimum. 


Landau’s  analysis  of  second-order  phase  transitions  can  be 
extended  to  describe  first-order  transitions  which  are  nearly  of 
second  order.  Such  weak  first-order  transitions  are  describable 
by  a free  energy  of  the  form: 

4 = tQ  + a(T  - Tc)n!  + Bn*  + Cn6  (3) 

where  C > 0,  and  B is  negative  but  small.  Figure  3 depicts 
the  free  energy  at  various  temperatures,  while  Figure  4 exhibits 
the  temperature  dependence  of  the  order  parameter  and  inverse 
susceptibility.  The  order  parameter  makes  a small  jump  at  the 
transition  temperature  and  the  susceptibility  does  not  diverge 
at  T , although  it  exhibits  the  characteristic  Curie  behavior 
X = C/(T  - T„)  with  (T^  < T„)  in  the  high  symmetry  phase. 


Figure  3:  Free  energy  *(n,T)  vs.  order  parameter 

n for  various  temperatures  above  and  below  critical 
temperature  T where  a first-order  phase  tran- 
sition occurs. 


V 


Figure  4:  (a)  Temperature  dependence  of  the 

magnitude  of  the  order  parameter  n for  a system 
which  undergoes  a first-order  phase  transition  at 
T . The  graph  gives  the  value  of  n at  the 
absolute  minimum  of  the  free  energy  which  is  shown 
in  Fig.  3.  (b)  Temperature  dependence  of  the 

inverse  susceptibility  x . for  a system  which 
undergoes  a first-order  phase  transition  at  T . 
This  graph  gives  the  value  of  the  curvature  of  the 
free  energy  (shown  in  Fig.  3)  at  its  absolute 
minimum. 


The  second  major  ingredient  of  the  Landau  theory  of  second-order 
phase  transitions  is  an  analysis  of  the  symmetry  properties  of 
the  order  parameter  and  how  these  relate  to  the  overall  symmetry 
of  the  system  both  above  and  below  the  transition  temperature. 
Let  us  first  consider  the  high  symmetry  phase  ( n = 0 ) of  the 
system  precisely  at  the  transition  temperature.  This  phase  can 
be  characterized  by  the  electric  charge  density  function  p(r). 
Since  the  charge  density  embodies  an  entire  description  of  the 
material,  one  could  in  principle,  calculate  the  free  energy  « 
from  p(r).  Of  course,  one  cannot  in  practice  perform  such  a 
calculation  but  we  can  symbolize  the  dependence  bv  writing 
$[p(r)].  Moreover,  since  the  charge  density  function  reflects 
the  spatial  distribution  of  all  the  electrons  and  nuclei  which 
make  up  the  material,  it  should  be  left  invarient  by  all  the 
crystallographic  symmetries  (rotations,  reflections,  tran- 
sitions, etc.)  possessed  by  the  high  symmetry  phase  of  the 
system. 

As  the  temperature  changes  and  the  system  undergoes  its  phase 
transition  into  the  new  (n  ^ 0)  phase,  the  electric  charge 
density  function  will  change  continuously  from  p(r)  to  p(r)  + 
6p(r).  For  a second-order  phase  transition,  the  electric  charge 
density  will  not  change  abruptly  at  the  transition  temperature; 
rather  6p(r)  will  start  from  zero  at  the  transition  temperature 
and  gradually  increase  in  magnitude  as  the  temperature  moves 
further  into  the  new  (n  ^ 0)  phase.  What  changes  abruptly  at 
the  transition  temperature  is  not  the  magnitude,  but  rather  the 
symmetry  of  the  charge  density  function.  As  we  shall  see  below, 
for  a second-order  phase  transition,  6p(r)  must  necessarily 
possess  fewer  crystallographic  symmetries  than  p(r).  So, 
as  soon  as  the  temperature  moves  an  infinitesimal  amount 
into  the  new  (n  ^ 0)  phase,  6 p ( r ) becomes  non-zero  and  the 
macroscopic  crystallographic  symmetry  of  the  system  is  reduced. 
(This  provides  the  reason  for  naming  the  n ^ 0 phase  the  low 
symmetry  phase,  and  the  n = 0 phase  the  high  symmetry  phase.) 
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Since  <5  p ( r ) brings  about  this  reduction  in  symmetry,  it  is 
instructive  to  analyze  6p(r)  in  terms  of  the  symmetry  of  the 
high  symmetry  phase.  We  shall  decompose  6p(r)  into  a sum  of 
functions  <^n^(r)  which  possess  special  properties  from 
the  point  of  view  of  the  crystallographic  symmetries  of  the  high 
symmetry  phase,  namely: 

«p(r ) = l ( I c[n)  *[n|r) ) (4) 

n i 

where  the  size  of  the  coefficients  of  the  normalized 

/ \ 1 

functions  determine  the  magnitude  of  5p(r)  at  a given 

temperature  (clearly  all  = 0 at  the  transition 

temperature).  The  real-valued  functions  <t>jn^(r)  are 
chosen  so  as  to  transform  according  to  a physically  irreducible 
representation  of  the  symmetry  group  of  the  high  temperature 
phase.  That  is,  when  acted  upon  by  a symmetry  transformation  of 
the  high  temperature  phase,  a given  function,  $jn^(r) 
transforms  into  some  linear  combination  of  the  other  functions 
4>[nNr),  which  have  the  same  index  n. 

For  example,  $!2^(r)  may  become  equal  to  (2<t>j2^(r) 

1 2)  ( 2 ) 

- <J>2  ( r ) + $2  ' (r)//5“  when  its  argument  (r)  = (x,  y, 

z)  is  replaced  by  (x,  y,  - z),  i.e.,  when  it  is  transformed 
according  to  a mirror  reflection  in  the  x-y  plane.  In  addition, 
this  decomposition  of  6p(r)  in  terms  of  families  of  real-valued 
functions  which  transform  solely  among  themselves  is  assumed  to 
be  complete,  in  the  sense  that  a set  of  functions  with  given 
index  n cannot  be  further  decomposed  into  two  subsets  of  real- 


valued functions  which  transform  solely  among  themselves  when 
acted  upon  by  the  symmetry  operations  of  the  high  temperature 
phase  (Ref.  1). 


Ref.  1:  It  is  sometimes  possible  to  further  decompose  a physical' 

ly  irreducible  representation  by  allowing  the  introduction 
of  complex-valued  functions.  (In  such  a case  we  drop  the 
adjective,  physically,  and  call  them  simply,  irreducible 
representations.)  Since  we  are  decomposing  an  intrinsically 
real  function  6p(r),  we  do  not  carry  the  decomposition  further 
than  the  physically  irreducible  representations. 
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In  summary,  the  index  n labels  the  different  physically  irredu- 
cible representations  while  the  subscript  i distinguishes  among 
the  different  independent  functions  belonging  to  a particular 
physically  irreducible  representation;  the  number  of  linearly 
independent  functions  in  any  one  such  family  is  the  dimension  of 
that  physically  irreducible  representation. 

It  is  possible  to  view  Equation  (4)  from  two  different 
perspectives.  Up  to  this  point,  we  have  regarded  the  coeffic- 
ients C[n)  as  constant  and  the  functions  4>j^(r)  as 
transforming  under  symmetry  operations.  An  equally  valid 
prespective  is  to  regard  the  functions  ${n^(r)  as  fixed 
and  then  allow  the  coef f icients to  transform  among 
themselves  under  the  action  of  the  symmetry  operations.  Hence- 
forth, we  shall  employ  the  second  perspective  exclusively:  The 
coefficients  are  regarded  as  transforming  according 

to  the  irreducible  representations  of  the  symmetry  group  of  the 
(n  = 0)  phase. 

The  reason  so  much  attention  is  focused  on  these  coefficients  is 
simply  that  the  coefficient  are  just  the  components 

of  the  order  parameter. 

Each  separate  set  of  coefficients  with  a given  index  n,  (cjn^, 
C^0^...),  is  taken  to  constitute  a separate  order  para- 
meter, distinguishable  from  all  others  by  the  way  it  transforms 
under  the  symmetry  operations  of  the  high  symmetry  phase. 

Clearly  there  are  as  many  different  kinds  of  order  parameters  as 
there  are  different  physically  irreducible  representations  of 
the  symmetry  group  of  the  high  symmetry  phase. 

While  a plethora  of  order  parameters  is  thus  logically  possible, 
it  would  be  extremely  unlikely  if  more  than  one  order  parameter 
were  to  have  precisely  the  same  transition  temperature.  This 
restriction  to  phase  transitions  describable  by  a single  order 
parameter  amounts  to  a substantive  restriction  on  the  kind  of 
symmetry  change  that  may  take  place  in  a second-order  phase 
transition . 
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For  the  case  of  only  one  order  parameter,  we  may  drop  the  index 
n and  write  Equation  (4)  as. 


» * 


«p(r)  = r C-^U). 


(5) 


It  is  conventional  abuse  of  language  for  "order  parameter"  to 

refer  to  both  the  entire  array  of  coefficient  n — (C, , C,,...) 

? 1/2  ^ 

and  to  the  magnitude  of  these  components  n = (jet)  . 

The  expansion  of  the  free  energy  as  a power  series  in  the  order 
parameter  takes  the  form: 


•lotfCj*;]  = »0  + VllCj)  + ♦2P2(Ci)  + 

...  + Vn(Ci>  + l6) 

L,  U 

where  pn(C^)  is  a polynominal  of  the  n1"  degree  in  the 
components  of  the  order  parameter.  (Hence,  the  magnitude  of 
Pn(Ci)  goes  as  nn)-  Since  the  free  energy  must  be  an 
invariant  quantity  under  the  crystallographic  symmetry  trans- 
formations, each  term  which  appears  in  the  expansion  must 
likewise  be  invariant.  A given  order  parameter  has  specified 
transformation  properties  under  the  crystallographic  symmetries 
of  the  high  symmetry  phase,  so  it  is  either  possible  or  not 
possible  to  construct  an  invariant  polynominal  of  some 
desired  degree  ( pn ( C ^ ) ) . Thus  the  transformation  properties 
of  the  order  parameter  will  determine  whether  a given  term  will 
actually  appear  in  the  expansion  (6).  Moreover,  it  is  unreason- 
able to  expect  that  a term  which  is  allowed  by  the  symmetry  of 
the  order  parameter  would  not  appear.  In  summary,  then,  it  is 
the  symmetry  of  the  order  parameter  which  determines  which  terms 
appear  in  the  free  energy  expansion. 
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We  have  seen  in  Par.  1.1.1,  that  the  odd  degree  terms  must 
not  appear  in  a free  energy  which  is  to  describe  a second-order 


phase  transition.  This  requirement  places  further  substantive 
restrictions  on  the  possible  order  parameters  which  describe  a 
continuous  phase  transition.  In  fact,  the  requirement  that  the 
first  degree  term  riot  appear,  rules  out  continuous  phase  trans- 
itions except  those  which  reduce  the  symmetry  of  the  initial 
phase.  (This  fact  has  already  been  assumed  implicitly  in  the 
preceeding  discussion.) 


In  summary,  then,  the  Landau  analysis  of  continuous  phase 
transitions  severely  limits  the  allowable  transitions,  to  those 
describable  by  a single  order  parameter  which  possesses  no  odd 
degree  invariants.  It  is  surprising  how  far-reaching  these 
relatively  simple  requirements  are. 
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Soft  Modes  and  Structural  Phase  Transitions 


The  soft  mode  perspective  provides  a more  intuitive  microscopic 
picture  of  the  order  parameter  which  describes  a continuous 
phase  transition  in  a crystal.  Here  basically  one  identifies 
the  order  parameter  with  the  amplitude  of  some  optical  phonon. 
Such  a description  is  valid  for  a purely  structural  transition 
involving  the  atomic  positions;  it  is  a less  fruitful  approach 
for  transitions  in  which  the  electrons  play  a significant 
role . 


To  appreciate  the  insights  that  this  perspective  does  provide, 
let  us  consider  an  idealized  crystal  and  focus  our  attention  on 
the  motion  of  some  single  atom  within  each  unit  cell.  The  net 
force  acting  on  the  atom  has  two  major  contributions;  first,  a 
local  potential  provided  by  all  the  other  atoms  within  the  same 
unit  cell,  which  causes  the  atom  to  perform  some  oscillatory 
motion  about  its  equilibrium  position;  and  second,  an  inter- 
action with  the  similar  atoms  in  neighboring  cells  which 
couple  these  oscillatory  motions  together.  There  are  two 
complimentary  prespectives  from  which  to  view  a crystal.  In  the 
first,  one  focuses  on  the  oscillation  of  individual  atoms  and 
the  interaction  among  these  individual  oscillators  - as  describ- 
ed above.  In  the  second,  one  focuses  on  the  traveling-wave 
collective  oscillations  of  all  the  atoms  together  - the  optical 
phonons.  When  all  the  local  and  interaction  potentials  are 
harmonic,  the  running  waves  of  atomic  oscillations  are  the  true 
normal  modes  of  oscillation  for  the  crystal,  i.e.,  each  travel- 
ing wave  has  a characteristic  wavelength  and  frequency  and  does 
not  interact  with  all  the  other  traveling  waves  of  atomic 
oscillation.  When  the  local  or  interaction  potentials  are  not 
harmonic,  the  traveling  wave  oscillations  are  not  true  normal 
modes  of  the  crystal  - they  interact  and  scatter  off  each  other. 
Nevertheless,  these  optical  phonons  provide  a fruitful  way  of 
analyzing  the  collective  motion  of  all  the  atoms  within  the 
crystal . 
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Each  optical  phonon  is  characterized  by  its  wavelength  as  a 
running  wave  x (or  equivalently  its  wave  number,  k = 2ir/x),  the 
frequency  of  the  collective  oscillation,  v,  (or  equivalently,  its 
angular  frequency  id  = 2irv)  and  the  amplitude  a^  of  the  oscil- 
lation (see  Fig.  5a). 

Two  particular  optical  phonons  are  of  interest  for  structural 
phase  transitions.  First  the  infinite  wavelength  (i.e.,  k = 0) 
phonon  which  describes  the  in-phase  oscillations  of  all  the 
atoms  up  and  down  together  (called  the  zone  center  or  r point 
phonon;  see  Fig.  5b).  Second,  the  phonon  with  wavelength  of 
two  lattice  spacings  (i.e.,  k = ir/a)  which  describes  the  out-of- 
phase oscillations  where  the  position  of  neighboring  atoms  are 
always  anti-correlated,  (called  the  zone  edge  or  zone  boundary 
phonon,  see  Fig.  5c). 

In  a structural  phase  transition,  the  equilibrium  positions  of 
certain  atoms  within  each  unit  cell  shift.  In  the  "soft  mode" 
picture  of  such  a transition,  we  focus  attention  on  the  optical 
phonon  whose  spatial  pattern  of  oscillation  corresponds  to  the 
spatial  pattern  of  the  atomic  shifts.  For  example,  in  a struc- 
tural transition  in  which  a central  atom  is  displaced  from  the 
center  of  each  unit  cell  in  one  particular  direction,  one  would 
focus  attention  on  the  zone  center  phonon  which  describes  the 
in-phase  oscillation  of  these  central  atoms  in  that  particular 
direction  (see  Fig.  5b).  In  the  high  symmetry  phase,  these 
central  atoms  are  all  oscillating  together  about  the  centers  of 
their  respective  cells  at  some  finite  frequency. 

As  the  transition  temperature  is  approached,  the  oscillation 
slows  down,  i.e.,  its  frequency  approaches  zero  as  shown  in 
Figure  6 and  7a.  It  is  this  slowing  down  of  the  optical  phonon 
that  is  referred  to  as  the  "softening"  of  this  optical  mode. 

As  the  crystal  passes  through  the  transition  temperature  into 
the  low  symmetry  phase,  this  soft  mode  is  described  as  "freez- 
ing in".  Of  course,  this  does  not  mean  that  the  oscillations 
cease;  rather,  it  means  that  the  center  about  which  each  atom 
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Figure  5:  (a)  Idealized  depiction  of  a typical 

optical  phonon,  i,e.,  collective  running  wave 
oscillation  of  ail  of  the  central  atoms  in  the 
unit  cell.  The  wavelength  X is  six  lattice 
spacings  a,  s-"  the  wavenumber  is  k = 2ir/6a. 

One  should  ima0ine  that  the  solid  atoms  are 
all  oscillating  about  the  center  of  their  cells  at 
same  frequency  w The  magnitude  of  this 
collective  oscillation  is  given  by  a^.  (b) 
Idealized  depiction  of  a zone  center  (X  * °°,  k = 

0)  phonon.  All  of  the  atoms  move  up  and  down 
together.  (c)  Idealized  depiction  of  a zone 
boundary  (X  = 2a,  k = ir/a)  phonon 
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Figure  7:  (a)  Temperature  dependence  of  the 

frequency  of  the  zone  center  phonon  which  softens 
as  the  consequence  of  a structural  phase  transtion 
at  Tc . (b)  Heuristic  depiction  of  the  temper- 

ature dependence  of  the  amplitude  of  the  oscillatory 
motion  of  an  atom  which  is  displaced  from  its 
equilibrium  position  when  the  crystal  undergoes  a 
phase  transition  at  . The  solid  curve  gives 
the  temporally  averaged  position  of  the  atom;  the 
dashed  curves  are  the  envelope  which  marks  the 
extremes  of  its  oscillatory  motion. 
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oscillates  has  shifted  in  a spatial  pattern  which  mirrors  the 
spatial  pattern  of  the  oscillatory  motion  (see  Fig.  7b). 

It  is  possible  to  establish  a relatively  strong  link  between 
this  soft  mode  description  and  the  Landau  order  parameter 
description.  In  the  limit  of  a purely  harmonic  crystal  lattice 
whore  the  optical  phonons  are  true  normal  modes  of  the  system, 
one  can  classify  the  spatial  pattern  of  atomic  displacements  of 
each  normal  mode  according  to  the  crystallographic  symmetries. 

In  fact,  one  can  identify  each  optical  phonon  by  the  physically 
irreducible  representation  under  which  it  transforms.  Thus,  its 
behavior  under  the  crystallographic  symmetry  transformations 
enable  one  to  link  a particular  optical  phonon  to  the  particular 
order  parameter  which  transforms  according  to  the  same  physical- 
ly irreducible  representation. 
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1.3 


Extrinsic  Ferroelectrics 


The  figure-of -merit  for  a pyroelectric  detector  material  is  p/ec , 
where  p = (9P/9T)  is  the  pyroelectric  coefficient,  e = (9P/9E)  + 1 
is  the  dielectric  constant,  and  c is  the  volume  specific  heat.  In 
an  intrinsic  (or  proper)  ferroelectric  for  temperatures  where 
the  pyroelectric  coefficient  is  large,  the  dielectric  constant 
is  also  large.  This  is  readily  understood  since  in  an  intrinsic 
ferroelectric  the  same  structural  deformation  of  the  material 
that  is  temperature  dependent  also  carries  the  ferroelectric 
polarization.  In  a regime  where  this  structural  deformation  is 
easily  changed  by  adding  microscopic  thermal  energy  (i.e.,  p = 

9P/9T  is  large),  this  structural  deformation  is  also  easily 
changed  by  an  external  electric  field  (i.e.,  £ = 9P/9E  is 
large)  which  can  couple  directly  to  the  structural  deformation 
via  its  electric  polarization.  Thus,  in  an  intrinsic  ferroelectric 
the  pyroelectric  and  dielectric  response  are  fundamentally 
related  since  they  owe  their  origin  to  the  same  microscopic 
structural  deformation.  Moreover,  and  unfortunately,  the 
relationship  is  such  that  a large  pyroelectric  coefficient 
implies  a large  dielectric  constant,  thereby  producing  a reduced 
detector  f igure-of -meri t . 

In  an  extrinsic  ferroelectric,  there  are  two  structural  deforma- 
tions taking  place  simultaneously.  The  deformation  which  is 
temperature  sensitive  (called  the  primary  or  soft  mode)  carries 
no  electric  polarization  itself.  It  is,  however,  coupled  to  a 
secondary  deformation  which  carries  the  electric  polarization. 

The  description  of  an  extrinsic  ferroelectric  must  include  both 
the  magnitude  of  the  primary  deformation (the  order  parameter,  n) 
and  the  magnitude  of  the  secondary  ferroelectric  deformation 
(the  polarization,  P).  Since  the  thermal  and  electric  response 
of  the  material  are  embodied  in  separate  microscopic  features  of 
an  extrinsic  ferroelectric,  one  no  longer  has  the  tight  coupling 
between  the  pyroelectric  coefficient  and  dielectric  constant' 
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which  occurs  in  intrinsic  f erroelectr ics . The  precise  relation- 
ship of  the  thermal  and  dielectric  response  depends  upon  the 
coupling  between  the  primary  soft  deformation  and  the  secondary 
ferroelectric  deformation.  While  this  coupling  can  vary  con- 
siderably, an  analysis  of  the  symmetry  of  the  order  parameter 
and  the  polarization  enables  one  to  classify  all  extrinsic 
f erroelectr ics  into  three  distinct  classes. 

(a)  Improper  Ferroelectr ics . The  primary  order  parameter  n 
which  softens  at  the  transition  temperature  completely  describes 
the  reduction  in  crystal  symmetry  at  the  phase  transition.  The 
order  parameter  n transforms  under  the  crystallographic  sym- 
metries differently  than  does  the  electric  polarization  P.  The 

n-P  interaction  contributes  to  the  free  energy  a term  of  the 
2 3 

form  n P or  n p.  Such  materials  have  dielectric  properties 
which  are  dramatically  different  from  ordinary  ferroelectr ics . 
Gadolinium  molybdate  is  the  best  known  example  of  such  a material. 

(b)  Pseudoproper  Ferroelectrics . The  primary  order  parameter  n 
also  completely  describes  the  symmetry  reduction  at  the  transition. 
However,  in  this  case  n transforms  under  the  crystallographic 
symmetries  in  the  same  fashion  as  the  electric  polarization 
(although  it  is  not  the  polarization).  The  n-P  interaction 

will  then  contribute  a term  nP  to  the  free  energy.  The  dielectric 
behavior  of  such  a material  is  qualitatively  indistinguishable 
from  an  ordinary  ferroelectric.  KDP  provides  the  best  studied 
example  of  a ferroelectric  transition  of  this  type. 

(c)  Triggered  Ferroelectrics.  The  primary  order  parameter  n 

does  not  completely  account  for  the  symmetry  reduction  at  the 

transition  (thus  the  transition  is  necessarily  of  first  order). 

The  transformation  properties  of  n under  the  crystallographic 

symmetries  are  different  from  those  of  the  electric  polarization 

and  are  such  that  the  lowest  order  n-P  interaction  contributes 
2 2 

a term  n P to  the  free  energy.  Such  a material  would  have 
dielectric  properties  quite  different  from  ordinary  ferroelec- 
trics, but  no  example  of  such  a transition  has  been  identified. 
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1.4 


Ferroelasticity  and  the  Reduced  Free  Energy 

All'  known  improper  f erroelectrics  also  exhibit  ferroelasticity. 
It  is  possible,  however,  to  formulate  the  theory  and  perform 
measurements  in  such  a fashion  that  one  does  not  have  to  deal 
with  the  ferroelasticity  explicitly.  An  improper  ferroelectric 
requires  a minimum  of  three  thermodynamic  coordinates  to  specify 
its  state,  viz.,  temperature  (T),  macroscopic  polarization  (P), 
and  the  order  parameter  describing  the  soft  mode  which  undergoes 
the  phase  transition  [n~(C^ ,C2»  • • ] • For  an  improper 
ferroelectric  which  is  also  f erroelastic , an  additional  thermo- 
dynamic cordinate,  x,  the  strain  field,  must  also  be  specified 
to  describe  the  thermodynamic  state  of  the  system.  The  free 
energy  describing  such  a system  has  the  form. 


4>  (n,P,x,T) 


<J>o  + T ao  (T~TC*  ri2  + T r'2f 4 + T n6f 6 

12  12 
+ ~2  n Pg2  + ~2  h X h2 

12  12 
+ -jAP*  + BxP+-±-Cx 


(7) 


where,  fn,  qn,  and  h^  are  temperature  independent  poly- 
nomials of  degree  n in  the  direction  cosines,  Y • = C./h 

1 2 1^2 

of  the  order  parameter  whose  magnitude  is  n = (£C^  ) 

On  the  first  line  are  the  usual  free-energy  terms  describing  a 
continuous  phase  transition  for  0.  The  terms  in  the  second 
line  describe  the  interaction  of  n with  P and  x.  The  terms  in 
the  last  line  provide  the  stabilization  of  P and  x and  a possible 
direct  interaction  between  these  two  coordinates. 

In  studying  the  pyroelectric  and  dielectric  properties  of  such  a 
material,  the  ferroelastic  features  are  not  of  direct  interest. 
Since  no  external  stresses  are  applied  to  the  sample,  one  can 
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use  the  condition  84>/8x  = 0 to  eliminate  the  strain  coordinate 
x from  the  free  energy  $ (n # P* x , T)  to  give  a reduced  free 
energy  $'(h,P,T)  which  only  depends  on  the  thermodynamic 
coordinates  of  interest.  This  reduced  free  energy  takes  the 
form , 


(n,P,T) 


ao(T-Tc)n' 


f 


6 


(8) 


where,  f^'  = f4  - (h2)2  /2C  , 
A'  - A - B2/C,  and 

g2"  = g2  " h2  B/C  . 


Thus,  the  reduced  free  energy  in  Equation  (8)  has  the  same  form 
as  the  complete  free  energy  given  in  Equation  (7)  if  all  terms 
containing  the  strain  were  omitted.  While  the  form  of  the 
reduced  free  energy  is  the  same  as  for  an  improper  ferroelectric 
which  is  not  ferroelastic,  the  renormalization  of  the  parameters 
given  by  Equation  (9),  which  is  brought  about  by  the  presence  of 
f erroelasticity , alters  the  physical  meaning  of  these  free 
energy  parameters. 
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